Elementary objects of the ID Hubbard model from its symmetry and exact solution I 
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In this first paper of a series of two, we uniquely define the exotic elementary objects that emerge 
from the interplay of the ID Hubbard model recently found global [5(7(2) <g> 5(7(2) ® U{l)]/Z$ 
symmetry with the exact Bethe-ansatz (BA) solution. Their occupancy configurations generate 
the exact energy eigenstates of the model. The results refer to the limit of a very large system in 
which the thermodynamic BA equations apply. The BA solution choses one of the infinite electron 
- rotated-electron unitary transformations such that rotated-electron single and double occupancy 
■ are a good quantum number for any finite interaction value. This allows a precise operational 

definition of three exotic elementary objects whose occupancy configurations generate all exact 
energy eigenstates: spin-less and 77-spin-less c fermions, spin-1/2 spinons, and r;-spin-l/2 ?7-spinons. 
Their occupancy configurations generate as well representations of the c hidden (7(1) symmetry, 
spin 5(7(2) symmetry, and 77-spin 5(7(2) symmetry algebras, respectively. Furthermore, it is found 
that the inverse scattering method BA L-matrix and thus the corresponding monodromy matrix 
^ ■ are naturally expressed in terms of the seven local generators of the model Hamiltonian electron- 

$-H " interaction term local gauge 5(7(2) (g) 5(7(2) (g> (7(1) symmetry and of the related spin-less and 

^ \ 7/-spin-less fermion creation and annihilation operators. Our symmetry-related formulation reveals 

that the M. Takahasi's bound states of electron n-pairs described by complex numbers A' and his 
magnon bound states of n-pairs described by complex numbers A refer to 77-spin-singlet configurations 
involving n anti-bound 77-spinon pairs and spin-singlet configurations involving n bound spinon pairs, 
respectively. In the second paper the scattering theory associated with the present elementary 
objects is constructed. That study reveals that there is no contradiction whatsoever between the 
present elementary-object description and those in terms of traditional spinons and holons. The 
former and the latter are normal ordered relative to the electron vacuum and an initial ground 
O ■ state, respectively. They are found to refer to distinct elementary objects associated with uniquely 

related scattering-state basis choices. The possibility of such different basis choices stems from the 
degeneracy of the excited energy eigenstates that span the subspaces of the corresponding scattering 
theories. 



PACS numbers: 03.65.Ca, 71.10.Pm, 71.27.+a, 03.65.Nk 



I. INTRODUCTION 



\ The Hubbard model with nearest-neighbor hopping integral t and on-site repulsion U is possibly the most studied 
7— I ■ lattice model of correlated electrons. It features electrons that can hop between lattice sites due to the finite hopping 
integral t. Such sites represent atoms, that are arranged in an ordered, crystalline pattern of well-defined geometry. 
When two electrons are on the same site, they have to pay the energy U due to their mutual repulsion. This introduces 
additional electronic correlations beyond those due to the Pauli principle. 

On the one-dimensional (ID) lattice the model is exactly solvable. When the number of ID lattice sites N a ^> 1 is 
very large the thermodynamic Bethe-ansatz (BA) equations introduced in Ref. [l[ apply. The corresponding general 
BA equations were first derived by the coordinate BAJ2]. The same BA equations were later on obtained by use of the 



alternative algebraic BA inverse-scattering method [3|, |4|. The ID Hubbard model is the simplest condensed- matter 
toy model for the description of the role of correlations in the exotic properties of quasi-lD materials It may be 
realized as well in ultracold atomic ID optical lattices @. 

The low-energy physics of ID correlated systems shows universal behavior associated with the so called Luttinger- 
liquid theory Q. At finite excitation energy, beyond the linear Luttinger-liquid regime, the separation of the electronic 
degrees of freedom into exotic different excitation branches prevails. In the case of the ID Hubbard model, the use 
of the pseudofermion dynamical theory (PDT) [HHI revealed that the metallic-phase exponents characterizing their 
finite-energy correlation function singularities are a function of the excitation momentum and differ signicantly from 
the predictions of the linear Luttinger-liquid theory Q . Similar results have been later obtained by a general method 
applicable to both non-integrable and integrable ID correlated systems [r34l7| . 

This paper is the first of a series of two QQ. We find that by combining the model global [SU{2)®SU{2)®U(l)]/Z% 
symmetry [lj| with the exact BA solution [H-Q, an elementary-object operator formulation emerges where such a 
global symmetry, which can be rewritten as [50(4) ® U(l)]/Z2, is made explicit. This allows the identification of the 
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exotic elementary objects whose occupancy configurations generate the model exact energy eigenstates both inside 
and outside the BA solution and define their operator algebra. The corresponding representation involves seven local 
operators whose site summation gives the seven generators of the model global symmetry. The local operator whose 
site summation gives the generator of the global c hidden U(l) symmetry beyond 50(4), is the local density operator 
of spinless c fermions. The local c fermion operators and the six local operators associated with the spin-1/2 spinon 
and 7y-spin-l/2 r/-spinon SU(2) algebras that emerge from the symmetry-related representation are mapped from 
suitable rotated-electron operators by an exact local transformation that does not introduce constraints. 

On the other hand, the rotated-electron operators are related to those of the original electrons by a unitary 
transformation such that rotated-electron double and single site occupancies are good quantum numbers for all finite 
values of u = U /At. There are infinite such transformations. The interplay of the model global symmetry with its 
exacts solution determines the unique choice of one of such transformations. From it elementary objects emerge whose 
occupancy configurations generate the exact energy eigenstates. For large on-site interaction values the ?7-spin degrees 
of freedom become unimportant for the low-energy physics and one recovers the usual charge-spin separation. The 
corresponding charge degrees of freedom refer to the c hidden U(V) symmetry rather than to the U{\) symmetry 
contained in the 77-spin SU (2) symmetry. 

The ?7-spin (and spin) and ry-spin projection (and spin projection) of the energy eigenstates are denoted by S v 
and S* 3 (and S s and Sf 3 ), respectively. The BA solution subspace is spanned either by the lowest- weight states 
(LWSs) or the highst- weight states (HWSs) of both ?7-spin and spin algebras, such that S^ 3 — —S a or S^ 3 — S a , 
respectively, for a = c,s [2(j] . We call Bethe states, the energy eigenstates inside the BA solution. In this paper we 
use the LWS representation of that solution. We emphasize that there is no contradiction whatsoever between the 
global [SU(2) ® 577(2) ® U(i)]/Z% = [50(4) ® U(l)]/Z 2 symmetry found in Ref. [H for the Hubbard model on any 
bipartite lattice and the results of Ref. [2(| , concerning the counting of the A Na energy eigenstates of the Hubbard 
model in the bipartite ID lattice. The important point is that it is confirmed in this paper that the BA solution 
accounts for the quantum number occupancy configurations that generate the representations of the c hidden U(l) 
symmetry algebra beyond 50(4). Moreover, the energy eigenstates outside the BA solution subspace are found to 
have exactly the same c hidden U(l) symmetry algebra representations quantum number occupancy configurations 
as the Bethe states from which they are generated by the off-diagonal 7/-spin and spin operator algebras. 

Although the studies of Ref. [13] have not explicitly considered the model c hidden U(l) symmetry algebra beyond 
50(4), they have used the BA solution, which accounts for it. Thus such studies have implicitly accounted for that 
symmetry, through their counting of the BA solution quantum number occupancy configurations found in this paper 
to generate the representations of the c hidden U(l) symmetry algebra. Interestingly, that the method of Ref. [2(| 
has counted all A Na energy eigenstates gives an exact proof that there are no extra non-Abelian symmetries of the 
ID Hubbard model beyond 50(4) = \SU(2) <g> SU(2)\/Z\. This is indeed the case, since the exact additional hidden 
symmetry found in Ref. [19] is an Abelian U(l) symmetry. 

The generator of the new global J7(l) symmetry algebra beyond 50(4) is the number of rotated-electron singly 
occupied sites operator [l9|, whose eigenvalue we denote by 25 c . Alternatively, such a generator can be chosen to be 
the number of rotated-electron unoccupied and doubly occupied sites operator, whose eigenvalue is 25^ = [N a — 2S C ]. 
The spinon description introduced in this paper is normal-ordered relative to the S v — N a /2;S S — 0; 25 c = 
electron vacuum, in that such a state is not populated by spinons as defined here yet for it the ?7-spinon number 
value is maximum and given by N a . Complementarity, the 77-spinon description is normal-ordered relative to the 
S v = 0; 5 S = 0; 25 c = N a absolute ground state, since that half-filling and zero-magnetization state is not populated 
by ?7-spinons as defined in this paper. For it the spinon number value is maximum and given by N a . The numbers 
25 c and 25^ values are also those of c fermions and c fermion holes, respectively, which live on a lattice identical to 
the original lattice. (Such number values also equal those of the spinons and 77-spinons, respectively.) The c fermion, 
spinon, and ^-spinon operator formulation introduced in this paper applies to all electronic-density and spin-density 
values. 

Here and mostly in the second paper, Ref. [l8j], the relation of the spin-1/2 spinons, ^-spin-1/2 ?7-spinons, and 
spin-less and 77-spin-less c-fermions of the present formulation to other model elementary-object representations is 
discussed. However, a revision of the vast literature about the model physics description in terms of exotic elementary 
objects is not among our goals. The elementary object representations of the ID Hubbard model may be classified 
into two large groups: 

- Elementary-object descriptions normal ordered relative to the electron vacuum, alike the representation in terms 
of c, sv, and cv pseudoparticles of Refs. [2l|, [22[, which is a generalization of that considered in Refs. [H, [H| for a 
well-defined model subspace. The elementary-object occupancy configurations of this type of descriptions generate 
the energy eigenstates from the electron vacuum. The elementary-object operator formulation defined in terms of the 
original electron operators in this paper belongs to this group of descriptions. 

- Elementary-object descriptions normal ordered relative to a well-defined initial ground state. This includes the 



3 



traditional spinon and holon descriptions [7|, I25l430l |. which in general share the same spinon definition as spin- 1/2 
excitations of a S s — ground state. An early example in the literature of spin-1/2 spinons with no charge whose 
spectrum is associated with a "hole" emerging under a transition from a S s — ground state to an excited state in 
a sequence of BA spin quantum numbers are the spin-1/2 color spinors introduced for the solvable ID Gross-Neveu 
model in Ref. [25]. A similar definition was used in Ref. [26j for the spin-1/2 spin waves of the also BA solvable ID 
isotropic Heisenberg antiferromagnetic model. The spin-1/2 spinon description of Refs. 28-30| is a generalization for 
the ID Hubbard model of such spin-1/2 color spinors and spin-1/2 spin waves. For the present model such traditional 
spinon and holon representations are of two types in what the concept of a holon and an antiholon is concerned. 
For those of which that of Ref. (28| is a representative a holon and an antiholon refers to a "hole" and a "particle" , 
respectively, created in the k'jS distribution of Ref. [l[ under the transition from a S s = ground state of arbitrary 
electronic density n to an excited state. On the other hand, for the descriptions of which that of Refs. [HI, H3| is a 
representative, a ?7-spin-l/2 holon and a 7y-spin-l/2 antiholon is an excitation of the S s = S n = absolute ground 
state with 7/-spin projection +1/2 and —1/2, respectively, associated with a "hole" created in the k'^s distribution 
of Ref. [l| under the transition to an excited state. The elementary-object occupancy configurations of this type of 
descriptions generate the excited energy eigenstates from a well-defined initial ground state. 

Our results confirm the validity of the pseudoparticle descriptions of Refs. plM^ . The relation to the traditional 
spinon and holon representations 0, H3] is addressed in the second paper. The studies of that paper have selected the 
results of Refs. [28-30] as a reference. The distinct elementary-object descriptions are found in that paper to follow 
from different choices of scattering states associated with the dressed S matrices and phase shifts extracted from the 
BA solution. For all the corresponding scattering theories the asymptotic scattering states must have a well-defined 
energy. This condition is fulfilled provided that such asymptotic scattering states refer to energy eigenstates of the ID 
Hubbard model. As discussed in the second paper, the degeneracy of the excited energy eigenstates of S s = (and 
S v = 0) ground states allows different choices of scattering states and thus of scatterers, scattering centers, and dressed 
S matrices. Such scatterers and scattering centers are elementary objects whose occupancy configurations refer to the 
spin (and fy-spin and charge) degrees of freedom of the involved excited states. The unitary transformations relating 
the scattering states of the present elementary-object representation to those of the traditional holons and spinons 
of the descriptions of Refs. [28| and [HI, H3| is studied in the second paper for the ID Hubbard model in subspaces 
spanned by excited states of a S s = ground state and the S a = S v = absolute ground state, respectively. 

The results of the second paper confirm that there is no contradiction whatsoever between the present elementary- 
object description and those in terms of traditional spinons and holons as defined in Refs. [28l - [30| . The corresponding 
distinct elementary objects refer to different choices of scattering states allowed by the above mentioned degeneracy 
of the excited energy eigenstates that span the subspaces where the alternative descriptions apply. For instance, the 
spin-1/2 spinons as defined within the present operator formulation and the descriptions of Refs. [28l - f30| are found 
to be completely distinct elementary objects. The elementary objects of Refs. (28| and [29|, 30j have been constructed 
inherently to be excitations of an initial S s = ground state, their number value vanishing for it. In the u 1 limit 
the spin-1/2 spinons of Refs. (28| and [H, |3(| in the excited states of the zero-magnetization and half-filled absolute 
ground state become the spin-1/2 spin-waves of Faddeev and Takhtajan. Those refer to the related ID spin-1/2 
isotropic Heisenberg antiferromagnetic model [26]. On the other hand, within our operator formulation that ground 
state is populated by N = N a spin-1/2 spinons, as defined in this paper. 

The use of our symmetry-related formulation and corresponding counting of the symmetry [SU(2) (g) SU(2) ® 
U(V)\jZ\ algebra representations in each model subspace with fixed c fermion, spinon, and ?7-spinon number values 
provides valuable information on the physical information stored in the occupancy configurations of the BA quantum 
numbers. For instance, the physical origin of the Takahasi's bound states of electron n-pairs described by complex 
numbers A' and his magnon bound states of n-pairs described by complex numbers A considered in Ref. [l|. We find 
that they refer to 77-spin-singlet configurations involving n anti-bound ?7-spinon pairs and spin-singlet configurations 
involving n bound spinon pairs, respectively. The former refer to the 77-spin degrees of freedom of n rotated-electron 
doubly occupied sites and n rotated-electron unoccupied sites. The latter to the spin degrees of freedom of n spin-up 
rotated-electron singly occupied sites and n spin-down rotated-electron singly occupied sites. The imaginary part of 
the complex rapidity numbers that emerge in the BA thermodynamic equations of Ref. [l| refers within our operator 
formulation to the attraction between pre-formed ?7-spin-neutral 77-spinon pairs or spin-neutral spinon pairs. 

The studies of this paper account for the model global symmetry and reveal that the former above spin-related 
2n-site rotated-electron occupancies have additional c hidden U(l) symmetry degrees of freedom, which are described 
by a number 2ro of c fermions. On the other hand, for those of the latter 77-spin-related 2n-site rotated-electron 
occupancies such additional degrees of freedom refer to a number 2n of c fermion holes. The spin- less c fermions carry 
electronic charge e and occupy the sites singly occupied by the rotated electrons. They are insensitive to whether 
such sites are occupied by up-spin or down-spin rotated electrons. That information is stored in the corresponding 
spin-singlet 2n-spinon configuration. The c fermion holes refer to the sites doubly occupied or unoccupied by rotated 
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electrons yet they do not distinguish rotated-electron doubly occupied sites from unoccupied sites. That information 
is stored in the corresponding ?7-spin-singlet 2ri-7y-spinon configuration. (Since we denote the electron density by n, 
the Takahashi's pair positive integer numbers n are in the remaining og this paper called v.) 

We are able to access the statistics in terms of the generalized principle of Ref. [3l[ of the operators associated with 
the composite 2^-spinon sv objects and 2^-77-spinon r\v objects in the above spin-singlet and ?y-spin-singlet configura- 
tions, respectively. (Within our notation, v = 1, ...,00 is the number of spinon and 7/-spinon pairs, respectively.) For 
each energy eigenstate, the number value of spinons and ?y-spinons bound and anti-bound, respectively, within such 
composite objects is fixed. On the other hand, the spinons and 77-spinons left over are found to be invariant under 
the electron - rotated-electron unitary transformation. Those are called in this paper unbound spinons and unbound 
77-spinons, respectively. 

The object representation of Refs. [U [22| was the starting point of the PDT [§4l]|. The operator formulation 
introduced in this paper is used in the second paper to both confirming the validity of the latter dynamical theory 
and deepening the understanding of the microscopic scattering processes behind it. Specifically, its studies address 
the related problem of finding the model natural scatterers and scattering centers. Their dressed S matrix is found 
to be behind the phase shifts that determine the PDT finite-en ergy co rrelation- function exponents momentum and 
interaction dependence [sj-Tllj] . The more recent results of Refs. fl2Hl4l [l7t find the same type of behavior in a wider 
class of ID correlated systems. 

The scattering theory introduced in the second paper is used in its studies to determine the effective charge and 
spin carried by the elementary objects. The sv fermions occupancy configurations are shown to carry the states 
spin currents whereas the unbounded spinons fully control the effects of the spin transverse fluctuations onto such 
currents. (The same applies to the effects of the 77-spin transverse fluctuations on the charge currents carried by the c 
fermion and r\v fermion occupancy configurations, which are controlled by the unbounded 77-spinon occupancies [32j.) 
Interestingly, due to zero-momentum forward-scattering between the c fermions and si fermions, for spin densities 
m 7^ (and electronic densities n ^ 1) the c fermions (and si fermions) can carry with them spin currents (and 
charge currents). Such currents are found in the second paper to be determined by the pseudofermion scattering 
theory phase shifts. On the other hand, the unbound ?7-spinon and unbound spinon occupancies are found to control 
the effects of the transverse ?7-spin and spin fluctuations onto the non-LWSs charge and spin currents, respectively. 

The paper is organized as follows: In Section [TT] the ID Hubbard model is introduced. In addition, our formulation 
in terms of rotated-electron operators such that single and double occupancy is a good quantum number for 11 ^ 
is defined. That such rotated-electron description naturally leads to the emergence of three elementary objects 
whose occupancy configurations generate the exact energy eigenstates is the subject of Section IIIII In Section IIVI 
the three corresponding effective lattices are introduced. The studies of that section also address the consequences 
of the transformation laws of the 77-spinons and spinons under the electron - rotated-electron transformation. The 
occurrence of hard-core ?7-spin-neutral composite 2i/-7y-spinon-anti-bound bosons and hard-core spin-neutral composite 
2i/-spinon-bound objects is confirmed in SectionfVl The BA quantum numbers are then found to label the sv fermions 
and r]v fermions, respectively. Those are generated from the corresponding hard-core bosons by a Jordan- Wigner 
transformation. The anti-binding and binding character of the anti-bound 77-spinons and bound spinons, respectively, 
is confirmed. The consequences of the transformation laws of the sv fermions and r\v fermions under the electron 
- rotated-electron transformation, the relation of the simplest excitations to the representations of the model global 
symmetry algebra, and the properties of the distinct types of elementary objects are issues also addressed in Section 
IVl Section [VI] is devoted to the study of ID Hubbard model in physically interesting subspaces. Most of the results of 
that section focus on the quantum problem described by that model in the subspace with no rotated-electron doubly 
and unoccupied sites. In the m> 1 limit no contradiction whatsoever is found between the present spin- 1/2 spinon 
description and the ground-state normal-ordered spin-1/2 wave description of Faddeev and Takhtajan [26[. The 
concluding remarks are presented in Section [Villi Complementary useful information is provided in three appendices. 

II. THE ID HUBBARD MODEL AND ITS UNIQUELY DEFINED ROTATED-ELECTRON 

REPRESENTATION 

Here we shortly introduce the ID Hubbard model and its global symmetry found in Ref. 19]. Furthermore, we use 
the interplay of the model global symmetry with its exact BA solution to define the rotated-electron operators of our 
formulation. 
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A. The model and its global symmetry 



In a chemical potential [i and magnetic field H the Hubbard Hamiltonian under periodic boundary conditions on 
a ID lattice with a site number N a 3> 1 very large and even can be written as, 

H H symm / f^a 5 a 3 5 (1) 



where 



L symm / 

a— 77, s 



H symm = t [f + Au V D ] , u = U/At ■ 5* 3 = -~[JV« - N] ; 5* 3 = -i[# t - 7V ; ] 



T = — 



E E [ C L c J+i^ + C J>] ! V° = £ (%'.t - 1/2) (fii.4. - 1/2) ■ (2) 



In this equation and throughout this paper de index a reads a =t, 4- when used as an operator label and a = +1/2, —1/2 
otherwise, respectively. Moreover, T is the kinetic-energy operator in units of t, Vd is the electron on-site repulsion 
operator in units of U, u = U/At is the electron on-site interaction in units of it, which is often used in this paper, 
He = 2/z, [i s = 2/ifiiJ, \ib is the Bohr magneton, and S^ 3 and 5f 3 are the diagonal generators of the 77-spin and spin 
SU(2) symmetry algebras [H, H3|, respectively. The operator ct (and Cj^ a ) that appears in the above equations 
creates (and annihilates) a spin- projection a electron at lattice site j = 1, N a . The operator hj >a = ct a Cj >a counts 
the number of spin-projection a electrons at such a lattice site. The electronic number operators read N = 53<7=t I a 
and N a = J3j=i^',o-- The momentum operator is given by P — 53<T=t 4. Sfc ^cr(k)k, where the spin-projection a 

momentum distribution operator reads n a {k) — ct c^, a and the operator ct (and Cfe i(T ) creates (and annihilates) a 
spin-projection er electron of momentum k. 

In this paper we use in general units of both Planck constant h and lattice constant a one. We denote the electronic 
charge by e and the lattice length by L — N a a — N a . The S a and S^ 3 values of the LWSs and HWSs of the r^-spin 
and spin algebras are such that S a = —S* 3 or S a = S^ 3 , respectively. Here a — r\ for ?7-spin and a — s for spin. The 
LWSs have electronic densities n = N/L and spin densities m = [N^ — N±] / L obeying the inequalities < n < 1 and 
< m < n, respectively. The description of the states corresponding to densities such that 0<n<l ] \ < n <2 and 
— n < m < n ; —(2 — n) < m < (2 — n), respectively, is achieved by application onto the latter states of off-diagonal 
generators of the 77-spin and spin SU(2) symmetry algebras. (The expressions of those are given below in Section 
III Bl ) The BA solution refers either to LWSs or HWSs of such algebras. Here we chose the LWS formulation of that 
solution. Hence for simplicity but without loss in generality, most of our results refer to initial ground states with 
densities in the ranges defined by the inequalities < n < 1 and < m < n. 

An exact result valid for the Hubbard model on any bipartite lattice is that for onsite interaction it ^ it has two 
global SU{2) symmetries (HIH, which refer to a global SO (A) = [SUi2) <8> SU{2)]/Z 2 symmetry [HIH. A recent 
study of the problem by one of us and collaborators reported in Ref. [19j , reveals that an exact extra global c hidden 
U(l) symmetry emerges for u ^ 0, in addition to the 50(4) symmetry. Specifically, the Hubbard model on a bipartite 
lattice, such as the ID lattice considered in this paper, has a global [SU{2)®SU(2)®U{1)]/Z% = [SO{A)®U{l)]/ Z 2 = 
50(3) ® 50(3) <8) U(l) symmetry. The origin of the global symmetry is a local gauge SU (2) <g) SU(2) g3 U (1) symmetry 
of the model Hamiltonian electron- interaction term first identified in Ref. 36]. That local symmetry becomes for 
finite U and t a group of permissible unitary transformations. The corresponding local U{1) canonical transformation 
is not the ordinary gauge U(l) subgroup of electromagnetism. It is rather a "nonlinear" transformation (36j . 

The electron - rotated-electron unitary transformation plays a key role in the model physics. The global symmetry 
is closely related to unitary operators V — V(u) and corresponding rotated-electron operators Cj >a — V* ch^V. For 
those single and double occupancy are good quantum numbers for u > 0. The investigations of Ref. |37| reveal 
that there are infinite unitary transformations and corresponding unitary operators V — V(u) and choices of rotated 
electrons. They refer to different choices of the states that span the subspaces with fixed 25 c values. For any of the 
infinite choices of corresponding unitary operators V = e~ s the following expression refers to the same operator, 

25 c = V s Q V . (3) 

Here Q and the related operator D count the number of electron singly-occupied sites and doubly-occupied sites. 
They read, 

^EE^f 1 "^); D = (N-Q)/2, (4) 

j = l <T=f,4. 
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respectively. The operator provided in Eq. © is the above-mentioned global c hidden U(l) symmetry generator 
whose eigenvalue 25 c — 0, 1, ... is the number of rotated-electron singly occupied sites [l9j . 

B. The specific unitary operator associated with the electron - rotated-electron unitary transformation 

performed by the exact BA solution 

A result of major importance for the present operator formulation is that although addition of chemical-potential 
and magnetic- held operator terms to the Hubbard model, as given in Eq. ([T]), lowers its symmetry, such terms 
commute with it. Therefore, the global symmetry being 5*0(3) ® 5*0(3) ® U(l) implies that the set of independent 
rotated-electron occupancy configurations that generate the model energy eigenstates generate representations of that 
global symmetry algebra for all values of the electronic density n and spin density m. Consistent, the number of such 
independent [SU(2) ® 5/7(2) ® J7(l)]/Z| = 50(3) <g> 50(3) ® U(l) symmetry algebra representations equals for the 
present model its Hilbert-space dimension, 4 Na jl9l |. 

The dehnition of the specific electron - rotated-electron unitary transformation performed by the BA solution is 
the first step for the introduction of our operator formulation. 

1. The Bethe states and corresponding n-spin and spin 5(7(2) towers of energy eigenstates 

It is convenient to express the electron numbers as, 

A t = N^ + n v -n s ; N± = Nf + n v + n s , 

/\r = A t + A; = N° + 2n v ; n a = S a + S* 3 = 0, 1, 2S a , a = 1,8. (5) 

Here Aj?, Aj 1 , and A denotes their values for the Bethe states, which as mentioned in Section U are the energy 
eigenstates contained in the BA solution subspace, 

N$ = ^-S, l + S.; Nf = ^-S v -S s ; N° = N a -2S n . (6) 

According to our choice of BA representation, such states are LWSs of both the ry-spin and spin algebras, so that for 
them n n — n s = in Eq. ([5]). The numbers n n > and n s > given in that equation may be expressed in terms of 
the 77-spin S v , spin 5 S , and general electron numbers of Eq. (|5|) as, 

nr, = Sr,-^(N a -N)=0,l,...,2S v ; n s = S s - i(JV t - N±) = 0, 1, 25 s . (7) 

The rotated-electron numbers equal those given in Eqs. (jTj) for the electrons. In addition, for them single and 
double occupancy are good quantum numbers, so that the Bethe-state electron numbers of Eq. ([6]) may be separated 
into two parts: Those corresponding to the 25 c rotated electrons that singly occupy 25 c sites and to the [A — 25 c ] 
rotated electrons that doubly occupy [A — 25J/2 sites, respectively. Hence for rotated electrons the numbers given 
in that equation may be rewritten as, 

A t ° = [25 c + 25 s ]/2+[A Q -25 c -25 r) ]/2, 
A° = [25 c -25 s ]/2+[A Q -25 c -25^]/2, 

A = 2S C + [N a - 25 c - 25 n ] . (8) 

Here [25 c — 25 s ]/2 (and [25 c + 25 s ]/2) gives the number of down-spin rotated electrons (and up-spin rotated electrons) 
that singly occupy sites and [N a — 2S C — 25 r) ]/2 is that of down-spin rotated electrons (and up-spin rotated electrons) 
that doubly occupy sites. 

Each site can accommodate two electrons of different spin projection. One can then associate with each site two 
local orbitals. The number of unoccupied local orbitals, unoccupied spin f local orbitals, and unoccupied spin J. local 
orbitals of a Bethe state are conserving numbers given by A' 1 ' = [2A a - A ], Nb° = [N a -N$], and a[ 1 ' = [N a — N°], 
respectively, where the numbers AP, AP, and A are those given in Eq. ©. 
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Generalization to non-LWSs gives, 





- N R,+l/2 


+ K,-l/2> 




Ni 


= ^,-1/2 


+ ^-l/2> 




N$ 


= ^,-1/2 


T R,+l/2 ' 






— /V s 




(9) 



Here JV* = [N a - JVJ, iVf = [7V - JV t ], and, 

Nh,±i/2 = [25 c ±2^ sT 2n s ]/2, 

^,±1/2 = [N a -2S c ±2S v T2n v }/2, 

Nr = [Nk + i/2+N^ 1/2 ]=2S c , 

n r = K,+i/2 + ^ ; -i /2 ] = [Jv- - . ( 10 ) 

are the following rotated-electron site-occupancy numbers, 

Number of spin f rotated — electron singly occupied sites ; 
Number of spin I rotated — electron singly occupied sites ; 
Number of rotated — electron unoccupied sites ; 
Number of rotated — electron doubly occupied sites ; 
Number of rotated — electron singly occupied sites ; 

Number of rotated — electron doubly occupied plus unoccupied sites , (11) 
which are good quantum numbers for u > 0. We denote the numbers ±1 ^ 2 and ±1 ^ 2 by 

N R%/2 = [2S c ±2S s }/2, 

N R°±i/2 = [N a -2S c ±2S v ]/2, (12) 

respectively, for the particular case of the Bethe states for which n s = n n = 0. 

Let {|^i ,i A ,«)} be a complete set of 4^" energy, momentum, ?7-spin, 77-spin projection, spin, and spin-projection 
cigenstates for u > 0. Such states are chosen so that those among them that are LWSs of both the 77-spin and spin 
algebras are the Bethe states. Here I a is a short notation for the set of four quantum numbers [2S C , S v , S s ,n v ,n s ]. 
The triangle index A in I a symbolizes the three symmetries: The c hidden U(l) symmetry associated with the number 
25c and the two SU (2) symmetries associated with the numbers S n , n v and S s ,n Sl respectively, in [2S Cl S ni S Sl n v , n s ]. 
The index l a represents all remaining quantum numbers, other than those, that are needed to fully specify an each 
energy eigenstate |*i ,i a ,«). 

The energy eigenstates of that set that are not contained in the BA solution subspace are generated from the Bethe 
states as follows, 

l*W° A ,«> • (13) 
Here, 

c a = (* Jo , J o k , u |(s-) n -(s+) n l*« ,,o klU ) 

= [n a l] JJ[2S a + l-/], n a = l,...,2S a ,a = r],s, (14) 

are normalization constants, the 77-spin (a = rf) and spin (a = s) off-diagonal generators 5+ and S~ are given below, 
and I a and Z^ stand for [2S C , S v , S s ,n v ,n s ] and [2S C , S v , S s ,0,0], respectively. (Within our notation, Z^ refers to 
values of the general index I a such that = n s = 0.) 

Analysis of the BA quantum-number occupancy configurations that generate the Bethe states, jo u ), reveals 
that for u > the corresponding wave functions smoothly and continuously behave as a function of u. Hence upon 



N 
N 



R, + l/2 
s 

R-l/2 



iv .R, + l/2 



1/2 

N S R 



l*i.,lA,U> 



n 



=(^) n ° 
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adiabatically increasing u from any finite value to the u — > oo limit, each Bethe state, \^i .io^ u ), continuously evolves 

into a corresponding Bethe state, \^i ot i° & .oo}> an d vice versa. We emphasize that at u = oo there are many choices 

of complete sets of energy eigenstates. Our analysis refers to the specific set of states {|^E'j oi i^ i00 )} generated from 

the set of finitc-u Bethe states {| v f , / oj ;^ jU )} upon adiabatically increasing u to the u — > oo limit. There is a unitary 

transformation connecting the Bethe states \^>i ol o^ u ) and I'J^^.oo)- It is found below that the corresponding unitary 

operator commutes with the six generators of the two SU(2) symmetries in the model global [SU(2)<EiSU (2)<g)?7(l)]/Z| 
symmetry. 

2. The electron - rotated- electron unitary operator 

Both the corresponding sets of 4^" states {\^i a ,i & ,u)} for a given fixed u > value and {|^ , ! ,i A , 00 )}, respectively, 
are complete. The exact BA solution quantum-number independence of u is consistent with the model Hilbert space 
being the same for all u finite values. Hence it follows from basic quantum-mechanics Hilbert-space and operator 
properties that for this choice there exists indeed exactly one unitary operator V — V(u). It is such that any u — > oo 
energy eigenstate |*i oj / A>0 o) is transformed onto the corresponding u > energy eigenstate | v E'/ ,/ A , u ) as, 

|*W Al „> = V-t|* WAj00 ) . (15) 

The energy eigenstates \^i ,i & , v ) = V '\^i ,i^,oo) (one for each value of u > 0) that are generated from the same initial 
u — > oo energy eingenstate \^i ,i &}0 c) belong to the same V tower. The transformation given in Eq. (fT5|) also relates 
the Bethe states l^^.i^.u) t° corresponding u — > oo Bethe states {|^z 0) ^ )00 )}- This is why such a specific unitary 
transformation is that performed by the exact BA solution. 
The corresponding rotated-electron operators are given by, 

cl a = VU) <!T V- c i)£r = Ft Cj> F; fi J> = c] > c J> . (16) 

They are such that rotated-electron single and double occupancy are good quantum numbers for u > 0. This follows 
from the symmetries of the Hamiltonian electron-interaction term. Those are such that all u — ► oo energy eigenstates 
of the set {| v E'z oi z Ai oo)} are as well eigenstates of the electron single-occupancy operator and double-occupancy operator 
of Eq. (fJJ. In the u — > oo limit, the Hilbert space is then classified in subspaces with different numbers of doubly- 
occupied sites. Each of the states {\^i ,i & ,oo)} is contained in only one such subspaces. 

The electron and rotated-electron vacuum |0) is invariant under the above unitary transformation. The u > 
energy eigenstates |vE , z ,/ A , u ) = V J{ \^i o j^ :00 ) are generated from it by rotated-electron occupancy configurations that 
are simply obtained from those that generate the corresponding states |\l/i ol i Al oo) = G|0). The generator onto |0) 
of the corresponding u > energy eigenstates |4 , i oj ; A , ii ) = ^l^io.iA.oo) simply reads G — V^GV. Such u > 
energy eigenstates may then be written as \^i .i & . u ) — G\0). The rotated-electron occupancy configurations that 
generate them from the electron and rotated-electron vacuum |0) are the same that generate the corresponding states 
\^i ,i & ,oo) — G\0), provided that electrons are replaced by rotated electrons as defined here. This is why rotated- 
electron single and double occupancy are good quantum numbers for u > 0. 

The operational representation of the ID Hubbard model introduced in this paper refers to the above specific unitary 
operator and corresponding rotated-electron operators of Eq. (|16p . As given in that equation, the unitary operator 
V = V(u) generates rotated-electron operators from electron operators. We thus call it electron - rotated-electron 
unitary operator. It is associated with one of the infinite unitary transformations considered in Refs. [19. |37|. Except 
for u — ^ oo , the rotated electrons emerging from any of such transformations are distinct from the electrons. Given 
the unitarity character of the transformations associated with the operators V — V(U/t), they may be used in two 
different physical problems, which however are mathematically and technically fully equivalent. Those are considered 
in Refs. [37], HH and in this paper, respectively. 

Specifically, in Refs. [37], HH it is considered that the rotated creation and annihilation operators of Eq. (TT5|) 
refer to electrons. Within that choice, except for u — > oo the Hamiltonian V* H S ymm obtained by rotating the 
Hamiltonian H symm of Eq. ^ is not the Hubbard model. Instead, it is a rotated Hubbard Hamiltonian of a model for 
which electron double occupancy and single occupancy are good quantum numbers. On the other hand, within the 
alternative physical problem studied here, the rotated creation and annihilation operators of Eq. (116[) refer to rotated 
electrons, which for finite u values are objects distinct from the electrons. Within it the Hamiltonian V' H symm V' 
is the Hubbard model expressed in terms of rotated-electron creation and annihilation operators. For that problem 
rotated-electron double occupancy and single occupancy are good quantum numbers for u > 0. 
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We associate with any operator O an operator O — V* O V . It has the same expression in terms of rotated-electron 
creation and annihilation operators as O in terms of electron creation and annihilation operators, respectively. Our 
convention is that marks placed over letters being a caret or a tilde denote operators. (A first exception are the 
electron operators of Eq. ([!()]) ■ which we denote by a and Cj i(7 rather than by cj g and c Ji<T , respectively.) 

The unitary operator V of our formulation is uniquely defined by its A Na x 4 Na matrix elements, i ,i^.u\VY& v ,v^,u) ■ 

That, as justified below, for u > the unitary operator V = V{u) commutes with the six generators of the global 
77-spin and spin SU(2) symmetries, implies that the matrix elements between energy eigenstates with different values 
of S v , S s , n^, and n s and thus of I a vanish. Hence the finite matrix elements are between states with the same 
values. We denote them by Vi ot i> , 

Vl ,l' = (^lo,ljx,u\V\^l',l & ,u) = (^l',l & ,oo\V j \^l a ,l & ,oo)* 

= (^l ,ljx,u\^l',l^,oo) = <*J',JA,ao|*i ,J A ,u)*- (17) 

For 11 — y 00 rotated electrons become electrons, so that the matrix representing the unitary operator V becomes the 
4 Na x 4 Na unit matrix whose entries read ( l I / i ( ,,i A ,oo|Vl v I / i',z' A ,oo) = $i ,i' o <>i & .i> & - Here and throughout this paper 8i d' o 
and ^; A ,i' A are the usual Kronecker symbols. 

In the following we define the electron - rotated-electron unitary operator V by expressing all its finite matrix 
elements Vi ol i' in terms of known amplitudes of the exact BA solution. 



3. The unitary-operator matrix elements between Bethe states 

We start by considering the matrix elements of the specific electron - rotated-electron unitary operator considered 
here between Bethe states. Let |aci, x^o; yi, y N a) denote the state in which the N° electrons are located at 

sites of spatial coordinates x%, ...,xj^o and the N? down-spin electrons among those are at sites of spatial coordinates 

yi, y N o. Consistent with the numbers given in Eqs. jg]) and (fT2|), let | 

^ R R,-l/2 N R,-l/2 

denote the state in which a number = 2S C of rotated electrons singly occupy sites of spatial coordinates 

xf , a^ s ,o, a number N^_^ 2 = [N a — 25c — 25 ?? ]/2 of rotated electrons doubly occupy sites of spatial coordi- 
nates xf, ...,x d vSI , and a number N^°_ 1/2 ~ [2S C — 2S s }/2 of spin-down rotated electrons singly occupy sites of 
spatial coordinates x s ^ , x s ^ Si0 . The relation of the rotated-electron operators to the electron operators is a 

"r,-1/2 

function of u. Thus the expression of the states \u; xf, x s N s ; xf , x d Vi0 ; X\ , x s ^ sM ) in terms of the above 

R N R,-l/2 N R,-l/2 

A^°electron states is u dependent. Furthermore, the sets of states {\u;x{, ...,x s N s ;xf, ...,x d iQ ;x^,...,x s l Z et0 )} 

R N R,-l/2 N R,-l/2 

arc different for different u values. Only at a fixed u value are such states orthogonal. (Since as u — > 00 rotated 
electrons become electrons, one has that |oo; x\ , x s N s ; xf, x d r, t0 ; x i , x s ^ s0 ) and \xi, x^a', Ui, •••> Un ) 

R N R,-l/2 R.,-1/2 ^ 

are two alternative notations for the same state.) 

The above energy eigenstates may be written in two alternative representations as, 



l<Xi<N, 



ft a ,l%,u( x i, -, x N°\Vl, ...,y N a)\xi,...,x N o;y 1 , ...,y N o) , 

= /] ft B ,l<L ,rt{ x li •••) x Nt> X Xi x< M r i> ' X l ' •"' X M 3 <° ) 

' R R.,-1/2 JV R -1/2 

\ s s d d s\. s\. \ /"| q\ 

X \U] Xi, X N s] Xi, X N n,o ',Xi ,...,X Sj0 /, 

R n R,-l/2 "r,-1/2 

respectively. In both cases the summation is over all xi,...,xpjo from 1 to N a . For u > 
the quantity /j oj j^ )U is the amplitude of the state \x\, x^a; yi, y^o) and fi ^^ r t that of the state 
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si si \ 

H.-1/2 "R,-l/2 



/io.i^.w^ii -,XN°;yi, -,yN°) = fa, ...,x N o;yx, -,yN°\^i o ,i & ,u) 

floA ,rt( x l J •••) ^JVS ^li-i^/iri.o i ^1 ) •••) 



J, H,-l/2 JV i?,-l/2 



(u;xl, ...,x s N s;xf, ...,x d v , ; , x a , l0 l^z,,^" . u ) - (19) 

S JV ii,-l/2 iV H,-l/2 A 

By using the first of the state expressions provided in Eq. (1181) in the matrix elements general expression of Eq. 
(1171) . those may be rewritten as, 

Vi ,i> = (*i 01 ^,u|V'l*JMi > u) = (*l.,li,ul*lMi,oo) 

= //t,;^,«( a; i'---' :z; A' o ;yi'---'yjvo)/i',i«,,oo(a;i J ---,a;Ar o ;yi,---,yAft')- (20) 

l<Xi<N a 

To arrive to this expression, we have used the orthonormality of the states \x\, Xj^o; yx, Un°)- 

For u > the state l^^/o^u) is an energy eigenstate and the electron and the S n — N a /2;S S — 0;2S C = 

rotated-electron vacuum |0) is invariant under V. It follows that V\0) — |0) and thus (0|V^ = (0|. That com- 
bined with the unitarity of V reveals that the amplitude /; ; o rt (xl, x s NS ; xf, x c L,, ', a?*, x** a>0 ) given 

A ' R R,-l/2 fl',-1/2 

by (u;xl,...,x s N s;x(,...,x d m0 ; x\ l , aio |*j„io ) equals the amplitude /, jo (x x , ...,x N o;yx, ...,y N o) = 

« JV H,-l/2 JV K,-l/2 A A + 

(ail, Xjyo; t/i, J/i\r° l^/o.i^.oo)) Y e t refers to a different notation. Indeed, 

fi o ,i ( ',rt(xi,---,x s N s;xf,...,x d ^, ;x x ,...,x , i0 ) = (u; x^, x^-s ; Xj, x^^, ; 2^ , 2; „ l^z^.z" ,u) 

A B JV K,-l/2 JV i?,-l/2 R JV fl,-l/2 JV R,-l/2 A 

= (0|[c 1)<71 Q 2 , ffa ...C J v0 (O . Ar0 ]|*j oi j0 v)U ) 

= (0|[t/t CljCT1 fft C2ct2 f...yt Cjy0 ^ o F]Ft|* io)ii)00 ) 

= (0|ci, CTi C2, - 2 ...C W o . n0 |*; o ^ jQO ) 

= (xi, ...,a:jvo|* ;oi ^ i00 ) 

= fl o ,l A ,oc( x i,-i x N°;yi,-,yNf)- (21) 

Here both the relations \^i ot i° tV ) = V - l^j^o^oo) of Eq. (fT5|) . which is also valid for Bethe states, and Cj tCr = c^ a V 
of Eq. (fT6|) were used. 

From the above results one finds that for u > and the rotated electrons generated from the elec- 
trons by the electron - rotated-electron unitary operator V as defined in this paper the amplitude fa ot p jrt 

is independent of u. This is due to an exact compensation of the u dependences of the states \ty la u ) 
and \u;xl,...,x s NS ;xf,...,x d r,M ', x* , £3,,o ), respectively. Specifically, by combining the following re- 

R "r',-1/2 "li, -1/2 

lation between the iV°electron amplitude fa o p u (xi, Xjya; yx, y N o) and A^°-rotated-electron amplitude 



fi ,i ( ).,rt( x i, x s N s ;xf, x d ; 2^ , 



si \ 

> i " ■ • I ■ X A[ s ' > 

R,— 1/2 \R,-l/2 



fi ,i .,u(xi,-.-,X]s[o;yi,...,y N o)— f ll o t (xl,...,x s N s;x d ,...,x d ^, ;Xx^,..., 

£.> 4. ^ — / a' H R.-1/2 



Ki'XjV,, 



,-1/2 JV H,-l/2 



x (a;i,...,a; A ,o;yi,...,?/j V o|rt;xf,...,x^a / ;xf ,...,a^L, '; , x s ^ 3>0 ), (22) 

4 JV H,-l/2 JV H,-l/2 

with the u independence of fa ,0 rt (x\ , ... ), one finds, 

' A ' R H.-1/2 N R,-l/2 

S- fi o ,i A ,u( x ^ ■■■i x N°;yi,-~,yN°) = 

E, / s / s I . d' d '. si' s\ '\ 

Jlo.l . ,rt\ x l i"-i x N a ,Xi ,...,X N n,o ,X 1 , X 3>0 ) 

A " R.-1/2 JV H,-l/2 

® I I . si s I d' d ' si' si '\ /f, \ 

x —{xi,...,x N o]y 1 ,...,y N a\ri;;x 1 ,...,x N s]X l ,...,x N r l ,o ; x{* , x+ ). (23) 

OU * R "R.-1/2 "R.-1/2 
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By use of both the second state expression provided in Eq. (fT5|) and the amplitude relation of Eq. (PHI) in the 
matrix elements general expression of Eq. (|1 T[) those may, alternatively to their expression given in Eq. (|20l) . be 
written in terms of rotated-electron amplitudes as follows, 

Vlo,l' B — 5^ 5^ f* a ,l° .rt( X l> X Ng i x li •••■> x n v ' i ' X /v=>o ) 

' A ' H R.-1/2 " n -1/2 

l<i,<Af a l<x' t <N a 1 
x fl'l°,rt( x l >-"i x Ng i x l > •••) ^aj-ij.o j ---j ^.vts.o ) 

A H JV H,-l/2 JV f?,-l/2 

v / . s s . d d si si \ . sf s ' . d' d s-j- \ /O/l^ 

« JV H,-l/2 "r,-1/2 R JV H,-l/2 JV R,-l/2 

Hence the whole it dependence of this V/ oj p expression is in the state-product factor 

/,,. „,s r,,s . „d d . si si I . si s '. d' d 1 . si' si '\ 

\U, J^i, X N s , d>i , A ,,,0 jX]^ . HAJjX]^ , Xjy-i , , j0 , .1-^ s ,o /• 

R JV H,-l/2 JV i?,-l/2 " "H,-l/2 "r,-1/2 

We proceed by expressing the matrix elements of the electron - rotated-electron unitary operator provided in Eqs. 
(|20| and ([24)) . which are a particular case of those given in Eq. (fT7|) . in terms of known quantities of the exact BA 
solution. This requires the expression of the amplitudes of Eq. (I19[) in terms of such quantities. 

The amplitude fi g ^ u (xi,...,x N o;yi,...,y N o) appearing in the matrix elements provided in Eq. ([20)) has been 

calculated by use of the exact BA solution [39j | Here we write it as, 

fi„,i%,u{ x i, -;X N o;yi, ...,y N o) = C u [gp(xi, x N o) x gp(yi, ...,y N o)] , 

g P { Xl ,...,x N .) = 1)*V^ 
p 

1 



,=1 k Pi x Qi 



(25) 



I 1 ' 

\]Y,i< Xi <N a \9p( x Ii -,x N o) g P (y u .:,y N o)\ 2 

where the permutation Q is defined by the inequality, 

1 < xqi < x Q2 < ... < x QN o < N a . (26) 

The Yip summation runs over all permutations P of the BA quantum numbers kpj considered in Refs. [39| and |40j . 
The function gp{y\, ...,y N o) in Eq. (|2"5|) obeys the equality, 

9p(Vi,— ,Vn°) = tpp(vQi,&Q2, —,<tqn°) , (27) 

where <pp{<tqi,(tq2, &qn°) is defined in Eqs. (2.7) - (2.10) of Ref. |39| . Our notation for it, g P (yi, Un°)> follows 
from that function only depending explicitly on the spatial coordinates y\, ...,y N o of the down spins in the series 
VQl,crQ2,—,crQ N o, in increasing order, 

1 < 1ft. < 2/2 < - < y N ? < N° . (28) 



We profit now from the equality /t oj ;o >r4 — fi ,i a ^,oo °f Eq. (|2ip to express the amplitude 
fl i° rt( x i j •••) x n s I x ii x< nr*i.° > x i ' v,o ) m the matrix elements expression provided in Eq. (f24|) in terms 
of known quantities of the exact BA solution. In the it — >• oo limit, the related amplitude /z oj z<^ j00 has a simpler form, 
first given in Eq. (2.23) of Ref. [40(. Combination of the results of that reference with the properties of the electron - 
rotated-electron unitary transformation confirms that fi o: [0 & ]rt is a function of the rotated-electron singly occupied sites 

of spatial coordinates xf,...,x^ s , rotated-electron doubly occupied sites of spatial coordinates xf,...,x t L r , i0 , and 

R ' "r',-1/2 

spin-down rotated-electron singly occupied sites of spatial coordinates x* , ...,x*T ai0 . For u > this it-independent 

B,-l/2 

amplitude reads, 

si. si \ 

' ,0 ) = 
.-1/2 

si si 



fl ,l° ,rt{ x l' •••! x Nt > x li •••) ^m' 1 ' ' ^l^' 

A R ly R,-l/2 N R,-l/2 



C r t [4>U(l){ x \i — > X Ng ) X < ^ > gt/f2')( a: l' —i X MV,o ) x 4 I SU(2)( X 1 ) ■■■iX S 3y0 )] , (29) 

« w 'H.-1/2 "R,-l/2 
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where. 



017(1) (^li ■■■> x Nf l ) 



si 



' ■•■' a; Ar= 



(-1) C 



(-1) Q det 



) 



R,— 1/2 

= 02(^1 , -,^a,0 ), 



(30) 



and 



^Jj2l< Xi <N a I0{7(l)( x l> x Af= ) <^S(7(2)( X 1' ^ATJ' i/2 ) ( ^S'C/(2) (^l^' — ■> X N B R ° 1/2 ^ 



(31) 



is a normalization constant. An important property is that the amplitude of Eq. (|29j) factorizes into the three functions 



^(l)^!'-'^)' M x 



i I •••! ^jyl.O 



), and fai^, 



Nt 



) associated with the c hidden U(l) symmetry, ?7-spin 



ft, -1/2 A 'R y -l/2 

SU(2) symmetry, and spin SU(2) symmetry, respectively. 

Concerning the function 4>u(i)( x i > ■■■j x n s )■> the permutation Q appearing in the expression of Eq. (|30|) arranges the 



spatial coordinates x\ , . 

the rotated electron with down spin must come first. The spatial coordinates x\ , x s N 



r/ N „ into non-decreasing order with the restriction that from two equal coordinates that of 

refer to the rotated-electron 
correspond to the doubly occupied 



singly occupied sites in the original lattice. The coordinates xf, ...,a;f, 



si 



sites in an effective lattice containing only the doubly and unoccupied sites. The spatial coordinates x^, 

refer to the down spins in increasing order in an effective lattice containing only singly occupied sites. The determinant 
in that function expression depends only on the spatial coordinates of the 2S C sites singly occupied by the rotated 
electrons and not on their spins, <jqi, <Jq2, ^Q2S C - 

Concerning the functions <f>i(xf, x c ' ) and (^2(^1 j ■•■> x tts,o )> for simplicity we omit here their expressions. 

H.-1/2 "r',-1/2 

Indeed, they have exactly the same form as those appearing in the u — > 00 limit in the wavefunction for electrons. 
Such functions are given in Ref. (Specifically, see Eq. (2.24) of that reference for u 3> 1 plus the text in the 

paragraph just after that equation.) They are such that the corresponding functions <p r su(2) 
4> s su(2) = 02 nave the following limiting behaviors, 



(-1)^-2^1/2^ and 



J SU(2) 



= 1 for N% = [N a - 2S C ] = 2S V ; 



<>SU(2) = 1 for N R = 2S c = 25, 



(32) 



In spite of the exact mathematical equality fi o j<^ rt = fi g jo^ t00 , the amplitudes fi o jo & rt and fi ot io^ <00 have different 
physical meanings. On the one hand, the amplitude /; ; o rt and the above rotated-electron spatial coordinates 



„d d 

> 1 ) •••) ft[V,0 



and Xi , . 



si 
C . 
Ni 



on which it depends are valid for the whole u > range. On 



"R.-1/2 "R,— 1/2 

the other hand, the value of the amplitude /; o ,;^ j00 is mathematically identical to that provided in Eqs. (|29j) and 

(|30p . yet physically it refers only to the u — ¥ 00 limit. Such an amplitude describes both electron and rotated-electron 
occupancies, which in that limit are similar quantum particles. 

The expressions given in Eqs. (f20|) - ([30|) uniquely define the matrix elements between Bethe states of the electron - 
rotated-electron unitary operator considered here. 



4- The unitary-operator matrix elements between energy eigenstates beyond the BA solution 

Here we derive the expression of the electron - rotated-electron unitary operator matrix elements between energy 
eigenstates beyond the BA solution. For those the numbers n n and n s of Eq. J7]) have finite values. 

We start by showing that the six generators of the global 77-spin and spin SU (2) symmetries commute with the 
electron - rotated-electron unitary operator. To achieve that goal, it is useful to express the kinetic-energy operator 
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T given in Eq. © as f = T + T +1 + T_ x . Here, 

£y = - 7 = 0,±1, 

To-,j,j> = ^[nj ,-<7 c] >CT Cjv )(T n-j/ _ CT + (1 - % _ CT ) c] j(7 Cj> (1 - _ CT ) + c.c] , 

(T 

= c] >CT CjV )(T (1 - hj> - a ) + hj'-e c], a Cj> (1 - Hj.-o-)] , 

a 

= ^[(l-nj-^ct^c^^nj^^ + il-nj^^c^^Cj^hj-^}. (33) 

a 

While the operator To does not change electron double occupancy, the operators T + i and T_i change it by +1 and 
— 1, respectively. 

For u > the expression of the operator 5* such that V — e~ s can be expanded in a series of t/U, 



0{t 2 /U 2 ) . (34) 



Although there are infinite choices for the operators V = e~ s and S, they share two important properties [19l. I37I. |38| : 
(i) To leading order in t/U all read —jj [T + i — T_i], as given in Eq. (J34J) ; (ii) Their operational expressions involve 
only the kinetic operators T , T +1 , and T_i of Eq. (f3"3")h Such exact properties apply to the specific electron - rotated- 

electron unitary operator V considered here. Since it commutes with itself, the equalities V = e~ s = V = e~ s and 
S = S hold. Hence both the operators V and S have the same expression in terms of electron and rotated-electron 
creation and annihilation operators, so that, 

\f + x - T_il + ... . 



3 U 



(35) 

However, although S = S, except for u — > oo one has that To ^ To, T + i 7^ T+i, and T_i 7^ T_i where, 
^7 = - XI ^r;,w'= 7 = 0,±1, 

^0;j,j' = ^2[nj.,-a ct ff Cj/ )CT _ CT + (1 - gt ff Cj V (1 - + C.C.] , 

a 

a 

T-x-,j,j> = ^[{l-fij-^c^c^^hf _ . + (l-fi i /_ <T )ct,, CT Cj>n i ,_ <r ]. (36) 

Furthermore, for finite u values one has that — jj [T +1 — T_x] ^ — jj [T +1 — T_i]. The equality S = S refers to the full 
expansion in terms of the two sets of three operators To, T+i, T_i and To, T+1, T_i, respectively. Only performing 
the sum of the infinite terms of the corresponding two expansions does that equality holds. 

To confirm that the three generators of the spin SU(2) symmetry, three generators of the j^-spin SU(2) symmetry, 
and the momentum operator P commute with the electron - rotated-electron unitary operator V — V, one uses the 
exact result that the unitary operator V can be solely expressed in terms of the three kinetic operators given in Eq. 
([33)) 0,113] • In Ref. the following twenty one commutators were found to vanish, 

[P,f l ] = [S*°,f l } = [S£,f l ]=0; a = r),s, l = 0,±l. (37) 

Although the algebra involved in their derivation is cumbersome, it is straightforward. The vanishing of such com- 
mutators then implies that the momentum operator and the six generators of the ?y-spin and spin algebras commute 
with the unitary operator V, 

[P,V] = [S* 3 ,V] = [S±,V] =0; a = V ,s. (38) 

Expression of a general operator O in terms of rotated-electron creation and annihilation operators includes a set 
of commutators involving the operators O — O V and 5. It reads, 

d = d + [6 ) s} + ~[[d,s} > s} + .... (39) 
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Combination of this general operator expression with the vanishing commutators of Eq. (|38p confirms that the 
momentum operator and the six generators of the spin and 77-spin SU(2) symmetries have the same expression in 
terms of electron and rotatcd-clcctron creation and annihilation operators. Specifically, the momentum operator 
reads, 



p = E J2 kc L c ^ = p= E E fcc ~L^- 

CT=t-i k c=t,4, k 



(40) 



Furthermore, the six generators of the partial 50(4) symmetry contained in the model global symmetry are given by, 



N a 



N„ 



Q?3 



5„ — 



E' 

3=1 
N a 

E- 

3=1 



s x a 3 E^::.- " 
3=1 



E 

3=1 



°3,a ' 



a = 77, s . 



(41) 



Here the local operators read, 



3,n 



(l-n 



3,t 



3:V 



.)•>) 



(n 



(42) 



The corresponding rotated local operators Sj and s where I = £3, ± are given below. They have similar expressions, 
except that the electron creation and annihilation operators are replaced by rotated-electron creation and annihilation 
operators, respectively. 

The seventh generator of the model global [SU(2) (g> SU(2) ® C/(l)]/Z| = 50(3) ® 50(3) <g> U(l) symmetry, which is 
that of the global c hidden U(l) symmetry, Eq. ([3]), can be written in terms of the rotated-electron number operators 



l 3,° 



A/°3 



1<T given in Eq. (fll)|) as follows, 



N„ 



25, 



E ^>( : 

<T=t4 



"3,~" ) 



(43) 



For u > it commutes with the ID Hubbard model [13 |. On the other hand, except in the u — > 00 limit and in 
contrast to the remaining six generators provided in Eq. (|41[) . the generator of the c hidden U(l) symmetry given 
here does not commute with the unitary operator V. This is behind the hidden character of such a symmetry. 

The six local operators §j a , where a = 77, s and / = X3, ±, and the also unrotated local operator Sj. c = 

S<7=t 1 ™3,<* (1 — ^j,— c) commute with the Hamiltonian electron-interaction operator Vd provided in Eq. 



WD,SLj = WD,S j ,c] = 0: 



a = 77, s . 



(44) 



Indeed, such seven local operators are the generators of the local gauge SU(2) ® SU(2) ® ?7(1) symmetry of the 
Hubbard Hamiltonian at t = [36[ . 

That the electron - rotated-electron unitary operator V commutes with the six generators of the global 77-spin and 
spin 577(2) symmetries implies that its following matrix elements between energy eigenstates have the same value, 



1 



n 

a— rj,s 

n 

a— 77, s 



1 



(45) 



This holds for the set of energy eigenstates {| , 3 / i ,2 A ,ti)} and {\^i> ,i &t u)} with the same Za values and thus n v = 1, 2S V 
and n s = 1, ...,2S S values. This matrix-element equality is confirmed on replacing (5~)™=y(5+)™ a by (S~ S^) na V 
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in the matrix elements of Eq. (|4"5j) and systematically using the SU(2) algebra commutator [S a , 5j] = — 2<5 QjCe 'S'^ 3 
until reaching the following expressions, 



(*i e ,j A ,»|V r |*r,i A ,«) = II F-<*i.,i A ,ul(^a^) n, 'V r |*iM2 



a=r),s 

= II c~^io,i%,u\( § n a ,o + [n a ^ l[[2S a + l-j'])V\*t, A , a ) 

a=ri,s a j' = l 

= (* loA jV\%, AtU ). (46) 

Equation (|T4)) was used to arrive to the last expression. Note that this is the same algebra as that behind the 
normalization of an energy eigenstate \^i a ,i^, u ), 



c 

a—r},s 



n c-( $ i„4»i(^ i o+k!]n[2s a +i-i , ])i\i«,«) 

!=r,,s a j'=l 

*I.,ft,«l*Io,J A ,t.) = 1 - ( 4? ) 



The matrix-element equalities given in Eq. (|46[) are consistent with the properties of the general amplitudes 
(xi, xat; yi, yjv^l^io.ZA,")- The latter are closely related to the amplitude (x%, Xpfo; y\, j/^o l^^o^) of the 
initial Bethe state l 1 ^;^ u ) of Eq. (1431) . This follows from the invariance under the electron - rotated-electron 
unitary transformation of the 77-spin and spin degrees of freedom of the rotated-electron occupancy configurations 
that generate the non-LWSs \^i ,i^, u )- Out of the N = N° + 2n n electrons of the latter states, except for an overall 
phase factor (—1)™'', the amplitudes (x±, xn] yi, ■ yN x \^i ,i & ,u) are insensitive to the spatial coordinates of the 
n n on-site electron pairs originated from n n unoccupied sites of the initial Bethe state |^j OJ i^ )U ). Furthermore, out 

of the ATj_ = + n v + n s down-spin electrons of the non-LWSs |^i ,i A ,u), such amplitudes are insensitive to both 
the spatial coordinates of the n v spin-down electrons in the on-site electron pairs originated from n„ unoccupied sites 
of the initial Bethe state \^i Bt io^ lU ) and the n s spin-down electrons generated by n s spin-flip processes of n s spin-up 
electrons of the same initial Bethe state. Hence the amplitudes (x\, ...,XN',yi, ...,2W 4 .|^ , ;„, ; A ,u) only depend on the 
spatial coordinates x\, a;jv -2S„ of the initial N a = [N a — 2S V ] electrons that are not generated from ?7-spin-flip 
processes and on those y\, yN a /2—S n —S 3 °f the initial N® = — S v — S s spin-down electrons that are not generated 
from spin-flip processes. Hence they may be rather rewritten as (sci, xpf a -2S r) ', yi, VN a /i-s -S s \^1 ,Ia,u) an d 



read, 



fl d & ,u(xi, ...,XN a -2S, 1 ] 2/1, — ) VN a /2-S v -S 3 ) = 
(Xi, ...,X Na -2S^;yi, -,yjV /2-S„-S,l*Jo,«A,u) = 

(-1)"" fl ot l° A , u ( x l> ■-,X N o;y 1 , —,y N °)\N°=N a -2S r ,;N°=N a /2-S ri -S s • ( 48 ) 



Here the Bethe-state electron amplitudes fi ot i° AtU (xi, x^o; y 1; y^o) are those of Eqs. (|2~5 
A similar analysis of the rotated-electron amplitudes involving the non-LWSs \^i ,i^, u ) leads to, 

/( ,i A ,rt ( x i ? ■••! x Ni ! x l ' x m v -° I -El" 1 ■■■] x m s .° ) = 

" J, K,-l/2 "r,-1/2 

(u; x 1 , x N g ; x 1 , x MV ,o ; x^ , x s \^i a ,i & .u) = 

R JV i?.,-l/2 iV ii,-l/2 

(— l) ,l >7 f „ t (x\, . .., x s N .; xf, . .., x d N r, ,0 ; x x , a; , (49) 

where the Bethe-state rotated-electron amplitudes f lg t o rt (xf , x s N s ; xf , x d Vt0 ; x\ , x^Z, :0 ) are those of 

A ' R N R,-l/2 N R,-l/2 

Eqs. dMD-dSID- 
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The energy eigenstates \^i a ,i^,u) may now be written in two alternative representations as, 

|*/o,Za,i*) = ^2 flo,l&,u(x%,-, XN a -2S v ',yi, yN a /2-S n -S e )\xi, —, XN a -2S v ]Vl, ■-, yN a /2-S n -S s ) , 

l<Xi<N a 



— 5 1 fl ,l&,rt(Xi, X N s] X 1 , X MV ,o ')Xi , ... 

z ' B Jv i?.,-l/2 



ATI 



l<Xi<N a 



x |it; a;*, x N s ; , x NV ,o 



R,- 1/2 



R.-1/2 



(50) 



respectively. In both cases the summation is over all x\, XN a -2S v from 1 to _/V a . 

The alternative expressions provided in Eq. (|50|) apply to all A Na energy eigenstates, including the Bethe states. 
From the above results one finds that the electron - rotated-electron unitary operator matrix element between two 
arbitrary energy eigenstates out of the A Na such states that span the Hilbert space is given by, 

(*Io,Ia,«I^I*IM{,.,«> = (*Io,Ia,u|*I',I' a ,oo) = <W A Vlo,V 

= ^a,; a X! fL,i & ,u( x i>—i x N a -2s, 1 ;yi,-,yN a /2-s, 1 -s s ) 

l<Xi<N a 

X fl',l & ,oo(xi,...,XN a -2S v ;yi,—,yN a /2-S rl -S s )- (51) 

Here the amplitudes are those given in Eq. (J3HJ)- Note that both amplitudes refer to the same n v value. Hence the 
product of the two phase factors of the form (— l) Ul1 appearing in Eq. ([41?]) gives 1. This together with the amplitude 
expression provided in that equation is fully consistent with the matrix-element equalities of Eq. (|45[) . 

The matrix element provided in Eq. (|51j) can alternatively be expressed in terms of the rotated-electron amplitudes, 
as in Eq. (|24|) for the Bethe states, with the result, 



(*Wa,«I^I*;',^,u) = (*Jo,iA,u|*J',^,oo) = ^ A ,i' A V lo,V B 

l<Xi<N a Kx'.<N a 



d . s\ 

i x mv<° i X 1 ) 

JV H,-l/2 



■ s 'l 

, 
Ni 



r i si s I . d' d '. si' sj. '\ 

X .] V ,1 A ,rt\ X 1 >---y X Nt ' 1 ) •••) •''ivrl.O ) 2<1 i-jl^o J 

R "r.,-1/2 "H.-1/2 

X (u; X±, !Ejys ; fl^, •••) x'L-n.o X-t , X Si0 

H R, -1/2 JV i?,-l/2 

x joo;:^ , ...,Xjys \x^ ,---,x N ^,o )X-^ ,...,:e^. Sj0 ). 



R.-1/2 



(52) 



For all 4 ° energy eigenstates l^i'^u) the rotated-electron amplitudes appearing here are given by, 

flo.l&^tfali x Ni,i X lt X M 7 l-° > X 1 s0 ) 



ii,-l/2 



C r t [4>U(1) i X l' X N°) X ^5(7(2) ( x l' X< 2i ••■i ^aji.o ) x ( t>SU(2) ( x l'i ■■•i x Si0 )] 

K v ' H.-1/2 JV R,-l/2 



(53) 



where, 



*SC/(2) 



Q 5^(-l) p e <E i=° fcpi ' E <>i = (-l) Q det 
p 

(Xi4 >'{„., ) = ( -1)^-^-2^+2^/2 4> ^ ^ }j 

5(7(2 



' -R.-1/2 R.-1/2 



R,— 1/2 



(54) 



Except for the extra phase factor (—1)™'' in the function ^mjW'^i ■••: a; '^»j,o )> the quantities appearing in Eqs. 



(1551) and (|5"4"]) are the same as in Eqs. (|29l) and (|3T)|) for the Bethe states. Those given here are valid for both the 
latter states and the non-LWSs generated from them. 

The alternative general expressions provided in Eqs. ([51]) and (|52p for the 4 Na x 4^° matrix elements between 
arbitrary energy eigenstates along with Eqs. (|2"5)) - (|2"51) . (|4"5)) - (|4T)1) . and ([53"]) - ([54"]) for the corresponding amplitudes 
uniquely define the electron - rotated-electron unitary operator and corresponding unitary transformation considered 
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in this paper. In the u — > oo limit rotated electrons become electrons, so that the electron - rotated-electron unitary 
operator becomes the unity operator. Consistent, analysis of the amplitude behavior confirms that in such a limit the 
matrix elements given in Eqs. (|51[) and (|52j) simplify to, 

lim (^i .i & , u \V\^i'.i',u) = 6i A » lim V lo>v = £, ,/ 6 lo n . (55) 



III. THE THREE ELEMENTARY OBJECTS WHOSE OCCUPANCY CONFIGURATIONS GENERATE 

THE EXACT ENERGY EIGENSTATES 

Here we show that all A N °- energy eigenstates are generated by occupancy configurations of three rotated-electron 
related elementary objects: the 77-spin-less and spin-less c fermions, Ty-spin-1/2 77-spinons, and spin- 1/2 spinons. 

A. Elementary quantum object operators emerging from the model algebraic structure and symmetries 

The local generator Sj. c = 2 CT =ti (1 — %',-er) given in Eq. (|43[) and the alternative local generator s^ c = 
(1 — §j tC ) of the Hamiltonian electron-interaction term gauge U(l) symmetry in SU(2) ® SU(2) ® i7(l) may be 
expressed as, 

Sj, c = h jiC = f] c 4 C ; s'l c = (1 - u iiC ) = /j,c fl c ■ (56) 

Here /J c and fj >c stand for the following creation and annihilation operators, respectively, of suitable spin- less and 
77-spin-less fermions, 

fj,c = C L I 1 ~ %' + Cj.tr %',-<7 , 

flc = c^ a (1 - hj,. a ) + (-If ct >CT _ a , J = (3/2 - a) = 1, 2 . (57) 

This equation refers to two representations for the spin-less and 77-spin-less fermions, which we call J = (3/2— a) = 1,2 
representations. As confirmed below, the monodromy matrix of the BA inverse-scattering method 0, 14ll - l43l | contains 
/t and fj. c operators referring to the two J — 1,2 representations. The 1 representation (and the 2 representation) 
corresponds to the choice tr =t (and a =1) in Eq. (|57| . However, the fermions whose operators are given in Eq. (|57|) 
are spin-less and 77-spin-less objects that do not carry spin a. Both representations are faithful, being alternative. 

The exact solution of ID integrable models can be reached by two different methods: The coordinate BA used 
by Bethe himself [44| and the more algebraic operator formulation, usually called inverse scattering method (45[. 
For the ID Hubbard model such an algebraic formulation of the Bethe states refers to the transfer matrix of the 
classical coupled spin model, which is its "covering" jUH^]. It is not among the goals of this paper to provide full 
details about the algebraic inverse-scattering BA method. One of its properties is that before arriving to the same 
BA equations as the coordinate BA, the former method involves creation and annihilation fields that act onto an 
extended and partially unphysical Hilbert space, larger than that of the ID Hubbard model [3J, |4|. The two types of 
operators given in Eq. (1571) associated with the ^7 = 1, 2 representations, respectively, act onto that extended Hilbert 
space. Within the inverse-scattering BA method, the unwanted and unphysical terms are eliminated by imposing 
suitable restrictions to the rapidities that change the nature of the fields. Such a procedure is equivalent to project 
the problem onto the model smaller physical Hilbert space. 

Concerning the operators fj c and fj tC , such restrictions impose the use of only one of the above ^ = 1,2 represen- 
tations. Therefore and except for the on-site L-matrices general expression, in this paper we consider only the J = 1 
representation such that, 

fj.c = c] tf (l-n jti .) + (-iycj^hj^, 

fj.c = Cj,f{l-n jti ) + {-iyct t h jti . (58) 

Within the inverse scattering method [3, 0, IHj], the central object to be diagonalized is the quantum transfer 
matrix, rather than the underlying ID Hubbard model. The transfer-matrix eigenvalues provide the spectrum of a 
set of generalized conserved charges [4l| - |43l |. From the point of view of the BA algebraic inverse-scattering method, 
a central role is ascribed to the monodromy matrix. (For further details on that matrix, see Refs. flErfii.) The 
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monodromy matrix Tjj'(X) corresponding to site indices j > j' where — 1, ...,N a can be expressed as a product 
of on-site L-matrices £j(\) as follows, 



3-3 



T j , f (X)= [] ^-i"(A), 3>f- (59) 

2"=1 

The spectral parameter A in the argument of £j(\) and Tj.j'(X) parametrizes the Boltzmann weights Q. For the ID 
Hubbard model a pseudo-vacuum state exists, on which the monodromy matrix acts as a lower triangular matrix. 
Such a state is an eigenstate of its diagonal elements. This has enabled the algebraic construction of the Bethe states 
0, H| , similar to those previously obtained by the coordinate BA 0, Q • 

In the algebraic construction of the Bethe states program carried out in Refs. [!, 0], the charge and spin degrees 
of freedom were solved into two separated steps, respectively. That has decoupled the monodromy matrix problem 
into that of two matrices, referring to the charge and spin degrees of freedom, respectively. The studies of Refs. [U, 0] 
used the HWS formulation of the BA. To solve the charge problem a one-site doubly occupied state, which here we 
denote by \ j, was used as local reference state. (Within our LWS formulation, a similar procedure involves instead 
a one-site unoccupied state, which we call \j, ©).) The corresponding global reference state | fi) was then defined 
by the tensor product | f-l) = Oj=i ® \j> t-l) ■ The latter state is an eigenstate of the transfer matrix. Indeed the 
triangular property is easily extended to a corresponding charge monodromy matrix. Its off-diagonal entries are some 
of the creation and annihilation fields. Their application onto | fl) generates the charge degrees of freedom of the 
energy-eigenstate configurations. On the other hand, to solve the spin problem the authors of Refs. [3|, |4| have used 
another BA. It is equivalent to the solution of the six-vertex model in the presence of inhomogeneities. The solution 
of the spin degrees of freedom involves the diagonalization of the auxiliary transfer matrix associated with a spin 
monodromy matrix 0, 0] ■ Again, the off-diagonal entries of that matrix play the role of the creation and annihilation 
fields that generate the spin degrees of freedom of the energy-eigenstate configurations. 

The local gauge symmetry of the Hamiltonian electron-interaction term and the related model global symmetry 
must be explicit in the algebraic operator formulation of its solution. This is in part why the solution of the ID 
Hubbard model by the algebraic inverse scattering method 0, was achieved only thirty years after that of the 
coordinate BA 0,0]. Indeed, it was expected that on decoupling the monodromy matrix problem into those of 
a charge and spin monodromy matrix, respectively, both such matrices would have the same traditional Faddeev- 
Zamolodchikov ABCD form, found previously for the related ID spin-1/2 isotropic Heinsenberg model (2(|[45|]. Such 
a prediction followed the expected spin SU(2) symmetry and 77-spin charge SU(2) symmetry associated with a model 
global SO(i) = [SU(2) ® SU(2)]/Z 2 symmetry [H|. If that was the global symmetry of the ID Hubbard model, the 
charge and spin monodromy matrices would indeed have the same ABCD form. 

Consistent with the local gauge SU (2) <g> SU (2) ®t7 (1) symmetry of the Hamiltonian electron-interaction term found 
in Ref. 36| and the related model global SO(3) <g> SO(3) <g> U(l) — [5'0(4) ® U(l)]/Z 2 symmetry found more recently in 
Ref. [19 1, all tentative schemes using charge and spin monodromy matrices of the same ABCD form failed to achieve 
the BA equations obtained in Refs. 0, 0| by means of the coordinate BA. The problem was solved in Ref. 0, by 
use of an appropriate representation of the charge and spin monodromy matrices, which allows for possible hidden 
symmetries. Indeed, the structure of the charge and spin monodromy matrices reported in the latter reference is able 
to distinguish creation and annihilation fields as well as possible hidden symmetries. That for u > the ID Hubbard 
model spin and charge degrees of freedom are associated with the SU(2) and U(2) — SU(2) (g) 17(1) symmetries, 
respectively, in [SU(2) <g) U(2)]/Z$ = SO{3) <g> SO(3) (g) C7(l), rather than merely with two SU{2) symmetries, is 
behind the different ABCD and ABCDF forms found in Refs. 0, for the inverse-scattering method BA solution 
spin and charge monodromy matrices, respectively: Due to the extra £/(l) symmetry in U(2) = SU(2) <£> U(l), the 
latter matrix is larger than the former and involves more fields. 

However, the various operators contained in the monodromy matrix enter in the expression for the Bethe states 
found in Ref. 4] in a complicated way. This renders difficult an intuitive physical interpretation. A full algebraic 
construction of the Bethe states as that carried out in that reference is not between the goals of this paper. A less 
ambitious yet related task performed in the following is the expression of the L-matrix Cj in the formal monodromy 
matrix expression provided in Eq. (|59[) in terms of only the seven generators of the local gauge SU{2) <£> SU (2) <g> £/(l) 
symmetry of the Hamiltonian electron-interaction term given in Eqs. (|42l) and (|56[) and related spin-less and 7/-spin- 
less fermion operators of Eq. (|57l) . This confirms the important role played by that local symmetry and the related 
model global [SU(2) ® SU(2) <g> U(l)]/Z 2 symmetry in the exact BA solution. 

By combining that result with the requirement that the Bethe states must correspond to representations of the 
model global symmetry, we find that the generators of such states onto the electronic vacuum may be expressed in 
terms of a set of operators obtained by electron-rotating such seven generators and related spin-less and 77-spin-less 
fermion operators. This confirms that although masked by the very complicated algebraic structure of the inverse- 



19 



scattering method, the BA solution performs the electron - rotated-electron unitary transformation defined in Section 

cm 

In the following we choose the same parameterizing functions a = a(X) = cos A and b = 6(A) = sin A as in Refs. 
0, Ell . Starting from the known expression of the on-site 4x4 L- matrix in terms of electron creation and annihilation 
operators, Eq. (|A1|) of Appendix [Al after some algebra relying on Eqs. (|42l) . (1551) . and (1571) we arrive to the following 
remarkable expression for Cj , 



— iab - 
?t t,h 



2 



-b / JiC h\ - i{-iy a f} c h 2 ib f jtC n\ - (~iy a /] c ht 



iab 



3,s 



-iab + ^-2~Sj 



s X3 



ie h {-iys 



.)•>] 



aflc^i-ii- 1 ) 3 b kc hi ia fj c h\ + (- l) j b f jtC h 2 



ie h {-iys- n 

iaf jtC n l l + (-iybfl c n 1 



iab 



—s h +s X3 

2 6 J,c T b j, V 



smh[2h] = sinh[2/i(A)] = u2ab; a = a(X) = cos A ; b = b(X) = sin A . (60) 

Here s l j V and s l j s with I — x^,± are the local generators given in Eq. (|42"|) of the two gauge SU{2) symmetries, 
respectively. Moreover, Sj :C and Sj c are the two alternative local generators defined in Eq. (|56")l of the gauge U(l) 
symmetry in the t = Hamiltonian gauge SU(2) (g> SU{2) ® U(l) symmetry. The operators fi2 = nj,f, Wy = (1 — ni), 
and nb; = (1 — h^) multiplying the spin-less and ?7-spin-less fermion creation and annihilation operators select one of 
their expression two operator terms given in Eq. (|57j) . To shorten the L- matrix Cj notation of Eq. (|60p . it is implicitly 
assumed that a spinless and ?7-spinless fermion operator multiplying an operator hj or h j with index J = 1 (and 
J = 2) corresponds to the 1 (and 2) representation of Eq. ([57|) . The L-matrix Cj expression of Eq. provides an 
explicit confirmation of the deep relation of the BA solution to the local gauge SU{2) <g> SU{2) ® U(l) symmetry of 
the Hamiltonian electron-interaction term associated with the vanishing commutators of Eq. (|4"4")l . 

We call c fermions the objects whose creation and annihilation operators /j = u c V and fj tC — fj tC V, 
respectively, are generated from those of the spin-less and 77-spin-less fermions of Eq. (|55|) by the electron - rotated- 
electron unitary transformation performed by the BA. Hence in terms of the corresponding rotated-electron creation 
and annihilation operators they read, 

fj,c = ^U^- hj^ + i-iy c^h u , 

fj,c = Cj.t i 1 ~ hj, 4.) + {-I) 3 ct f h u . (61) 
The c fermion local density operator is then given by, 

/;.,./:,,•• (62) 

(Exceptionally, within our notation no upper index as in /J c , /, and n is used onto the (electron-rotated) c fermion 
operators fj_ c and %, c , respectively.) 

One may introduce corresponding c fermion momentum-dependent operators labeled by discrete momentum values 

4,c=(4,c) t = ^E e+ ^"4,ci j = h-,N a . (63) 
Viv 0j ., =1 

For N a 3> 1 and in units of 2-k /N a , the discrete momentum values qj are found below to be the quantum numbers of 
one of the BA excitation branches. 

We consider as well the six rotated local operators s l j a = s l j a V with a = s,T] and I = X3,±. Those are 
generated from the remaining six generators of the local gauge SU(2) ® SU(2) ® U(l) symmetry of the Hamiltonian 
electron- interaction term given in Eqs. (|42[) . 
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The three rotated local spin operators 3j s and the three rotated local 77-spin operators §j _ such that I = ±, a; 3 given 
here are associated with the spin-1/2 spinons and 77-spin-l/2 77-spinons, respectively, as defined in this paper. 

On the use of the general results of Rcf. [19], one straightforwardly confirms that within the present operator 
formulation the site summation Sj=i °f the rotated local operators given in Eqs. (|62[) and (1641) gives the seven 
generators of the model global 50(3) <g> SO(3) ® U(l) = [SU{2) ® SU(2) ® U(l)]/Z$ symmetry. For the six operators 
obtained under that site summation from the six rotated local operators of Eq. (|64l) that is readily confirmed by the 
form of the six generators of the model global symmetry provided in Eq. (|4ip. On the other hand, concerning the 
generator of the global c hidden U(l) symmetry given in Eq. (|43p. a straightforward operator algebra confirms that 
the c fermion local density operator (|62|) is the local operator Sj c appearing in Eq. (|43|) . Hence the global c hidden 
f/(l) symmetry generator may be rewritten as, 

N a N a 

2S c = X> )C = £/j; c / J - c . (65) 

The three local spinon operators s and the three local 77-spinon operators s' are particular cases of general 
SU{2) local 77s quasi-spin operators gj. The latter operators refer to all lattice sites and are associated with a general 
quasi-spin SU(2) symmetry. The rjs quasi-spin operators — q* 1 ±i<jj 2 and gj 3 , where Xi,X2,Xa denotes the 
Cartesian coordinates, have the following expressions in terms of rotated-electron operators, 

Qj = (4,t + £ J'.t) Cj4 > 

g+ = (97 ) f ; 3f = (^-4-1/2). (66) 

Except for unimportant phase factors, the unrotated operators fj c = Vfj c V"^, fj. c — Vfj tC V\ qf = Vqfv\ and 
~ Vqj 3 ^ 1 * corresponding to those defined in Eqs. (|6Tj) and (l66l) are those previously considered in Refs. [46l - |48| . 
The same applies to the related three local spin operators Sj s and three local 77-spin operators si „ given in Eq. (|42|) . 
Hence again except for unimportant phase factors, in the u 3> 1 limit the operators defined in Eqs. (|4"2"]l . (I5T|) . and 
(f66|) become those of Refs. |46l448| . The advantage of our representation is that it is valid for the whole interaction 
range u > rather than only in the «> 1 limit. 

It is convenient to express the local operators Sj :C , Sj c , and Sj a where I = ±, X3 and a = s, 77 as, 

Sj.c — n j,c Sj c = (1 — Tlj iC ) , 

4> = n j>c q\ ; s l j n = (1 - Jij, c ) , ? = ±, x 3 . (67) 
Here n J;C is the c fermion local density operator given in Eq. (|62|) . Hence the 77s quasi-spin operators is given by, 

^=4> + 5 L' z = ±^3- (68) 

The expressions provided in the above equations are consistent with the spin 577(2) and 77-spin SU(2) symmetries 
being particular cases of the quasi-spin SU (2) symmetry. Within them that the three local spinon operators §j s and 
three local 77-spinon operators s^ „ are associated with SU (2) algebra representations defined in two independent sets 
of sites, (i) the 2S C spin-up and spin-down rotated-electron singly occupied sites and (ii) the 2S^ rotated-electron 
doubly-occupied and unoccupied sites, respectively, is made explicit. The c fermion and c fermion hole local density 
operators nj jC and (1 — rij jC ) in the expressions of the operators s l j S and Sj „ provided in Eq. (1671) play the role of 
projectors onto such two sets of lattice-site rotated-electron occupancies, respectively. 

The c fermion operators and 77s quasi-spin operators given in Eqs. (I5T|) and (p7|) and the related spinon and 77- 
spinon operators defined in terms of the local c fermion number operators and 77s quasi-spin operators in Eq. (|66[) 
are mapped from the rotated-electron operators by an exact local transformation that does not introduce constraints. 
Given their direct relation to the generators of the model global symmetry, their occupancy configurations naturally 
generate representations of the corresponding global symmetry algebra. Consistent, inversion of the relations without 
constraints given in Eqs. (|BT|) and (|6"6")l leads to, 

4 = //.cQ-^+c-iy/j-cQ+r) - 

qa = (hc + (-iYfl)qj. (69) 
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We confirm below that the c fermion operators commute with the 77s quasi-spin operators, what justifies the form of 
the expressions given here, whose c fermion operators position is always on the left-hand side. 

In Section III Bl the unitary operator that relates the rotated-electron operators to electron operators was uniquely 
defined. This assures that also the spin-1/2 spinon, ^-spin-1/2 ?7-spinon, and spin-less and 77-spin-less c fermion 
operators are within our formulation defined in terms of those of the original electrons. The M s = 2S C spin-1/2 
spinons describe the spin SU (2) symmetry degrees of freedom of the 2S C rotated electrons that singly occupy sites. 
The c hidden U(l) symmetry degrees of freedom of such 2S C rotated electrons are described by N c = 2S C spin-less 
and 77-spin-less c fermions. Those carry their charge e. On the other hand, the M„ = [N a — 2S C ] Ty-spin-1/2 77-spinons 
describe the 77-spin SU{2) symmetry degrees of freedom of the [N a — 2S C ] sites doubly occupied (77-spin-projection —1/2 
77-spinons) and unoccupied (77-spin-projection +1/2 77-spinons) by rotated electrons. The c hidden U(l) symmetry 
degrees of freedom of such [N a — 2S C ] sites rotated-electron occupancies are described by N r c l = [N a — 2S C ] spin-less 
and //-spin-less c fermion holes. 

Consistent with the form of the rotated-electron amplitudes given in Eqs. (|53| and (|54|) . we find below that the 
M s = 2S C spinons occupy a spin-effective lattice with M s = 2S C sites and the 77-spinons occupy a T^-spin-effective 
lattice with M v — [N a — 2S C ] sites. The c fermions c effective lattice is identical to the original lattice and has N c = 2S C 
occupied sites and N' c l — [N a — 25 c ] unoccupied sites. The important role played by the c hidden U(l) symmetry is 
confirmed by the number values M s = N c = 2S C and M v = — [N a — 2S C ] being determined by the eigenvalue 2S C 
of its generator. 

The rotated-electron interaction operator Vd = V* Vd V of Eq. ([71)1 and three rotated kinetic operators of Eq. 
(1361) have the following expressions in terms of c fermion and 77s quasi-spin operators, 



T 7 — ^7iii>j' ' -^TU'jj' ~~ F T>jtj' Q\f\;j,j' i 7 _ 0,±1 . (70) 

(hi') 



Here, 

^+i;j>J' = (~^y fj-cfj' ,c , 



h\;j,j' 



Q + + qpl + (-i) l7| 2?f , M = o, 1 . (7i) 



For finite u values, the Hamiltonian H symm of Eq. ^ docs not commute with the unitary operator V 



= e- s . 

Hence, alike the general operator of Eq. (f39|) . when expressed in terms of the rotated-electron creation and annihilation 
operators of Eq. (|16p. it has an infinite number of terms, 

•Hsymm — V Hsymm — H syrarri -\- [H syrnrn , S ] + — [[H syrnrni 5], S] ~\~ ... . (*72) 

The commutator [H symm , S] does not vanish except for u — > oo so that H symm ^ H symm for finite values of u. 

The Hubbard Hamiltonian given in Eq. (|72[) may be developed into an expansion whose order refers to the number 
of rotated kinetic operators T 7 , independently of their type, 7 = 0, ±1. To fourth order it reads, 

H 
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This corresponding Hamiltonian expansion can be expressed in terms either of rotated-electron operators or c fermion 
and 77s quasi-spin operators. (The latter expression is obtained by the use of Eqs. (l70l) and ([7T]).) 

First, in the present ID case the BA solution implicitly performs for the whole it > interaction range the sum of all 
Hamiltonian terms on the right-hand-side of both Eqs. (|72|) and (|73]) and beyond. The resulting formal Hamiltonian 
expression is given in terms of c fermion operators and operators of well-defined 77-spin-neutral 77-spinon composite 
objects and spin-neutral spinon composite objects below in Section lVEl The general Hamiltonian expressions provided 
in Eqs. (|72j) and (|73|) are valid for the Hubbard model on other lattices. 

Second, in the u»l limit the rotated-electron creation and annihilation operators in the Hamiltonian expression 
provided in Eq. (|73[) become electron creation and annihilation operators, respectively. In the large-u limit that 
Hamiltonian expression for the model on the ID lattice and other lattices was considered previously by many authors, 
as for instance in Refs. [38l |4^ . l50j . Apparently, such a Hamiltonian t/U expansion is formally similar in terms of 
electron operators, as given in such references, and in terms of rotated-electron operators, as provided in Eq. (|73l) . 
However that is only so for very large u values. For instance, if for intermediate u values one expresses the t 2 /{/, 
t 3 /U 2 , and t 4 /U 3 terms of the Hamiltonian given in Eq. (|73|) in electron creation and annihilation operators, one 
finds many complicated higher order t 3 /U 3 terms where j > 2, j > 3, and j > 4, respectively. Indeed for small 
and intermediate u values the first few terms of the Hamiltonian expression in terms of rotated-electron operators 
describe many higher-order electron processes. 

The use of our operational representation renders the intermediate- u quantum problem in terms of rotated electrons 
similar to the corresponding large-it quantum problem in terms of electrons. As confirmed in later sections, the 
effect of decreasing u involves an increase of the energy bandwidth of an effective spin band associated with neutral 
two-spinon occupancy configurations. The intermediate- it rotated-electron processes may be associated with exchange 
constants describing rotated-electron motion touching progressively larger number of sites. Within the rotated-electron 
representation used here, such Hamiltonian terms emerge naturally upon decreasing the magnitude of u. For moderate 
u the rotated-electron operators generate a much simpler form for the energy eigenstates as well as for complicated 
processes involving a large number of electrons. As mentioned above, in the large-u limit the rotated-electrons become 
electrons, so that the usual picture of Refs. [H, H^, [5(| is recovered. 

We call off-diagonal terms of an operator expansion in terms of rotated-electron creation and annihilation operators 
as that given in Eq. (I39[) . those that do not preserve the numbers of rotated-electron singly and doubly occupied 
sites. An interesting technical detail is that up to third order all the diagonal terms of the expansion provided in Eq. 
([73"]) are generated by the rotated-electron leading-order term, — {t/U)\T+\ — T-i], of the operator S. This is because 
when expressed in terms of electron operators the Hubbard Hamiltonian does not contain any off-diagonal terms with 
more than two electron operators. (In this case the off-diagonal terms refer to electron doubly occupied sites.) 

For the Hamiltonian expansion in terms of rotated-electron operators (and thus of electron operators for «» 1) 
only the Hamiltonian terms H^°\ H^ x \ H^ 2 \ and i?( 3 ' to third order given in Eq. (IT3"|) are universal. The form 
of the terms of fourth and larger order is different for each electron - rotated-electron unitary transformation. For 
the fourth-order term given in that equation only the real-number parameter 9 value is not universal, being 
unitary-transformation dependent (38|. For instance, the methods of Refs. [4!| and [5(| refer to two different large-it 
electron - rotated-electron unitary transformations whose 9 values are 9 = and = 1/4, respectively. Moreover, one 
of the methods of Rcf. [38] refers to an electron - rotated-electron unitary transformation whose 9 value is 9 = 1/2. In 
this paper the value of the parameter 9 associated with the specific electron - rotated-eleectron unitary transformation 
performed by the BA solution is not computed. Indeed, as confirmed below in Section l"V El that solution performs 
the sum of all Hamiltonian terms on the right-hand-side of Eq. (|73[) and beyond, so that such a value is not needed 
for our studies. 

Note that the rotated-electron on-site interaction operator, 

Vd = V* Vd V = ]T V DJ ; V DJ = (n i)t - 1/2) (h u - 1/2) , (74) 
j=i 

which in terms of c fermion operators has a non-interating form given in Eq. (|70[) , appears only once in the Hamiltonian 
expression provided in Eq. (|73l) . Indeed the Hamiltonian terms of order higher than zero involve only the three kinetic 
operators To, T_i, and T+i . This is so because the derivation of such a Hamiltonian expression relies on the systematic 
use of the commutator, 



[Vb,f 7 ]=7T 7 , 7 = 0,±1 



(75) 
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B. The elementary objects operator algebra 

Since the electron - rotated-electron transformation associated with the operator V is unitary, the rotated-electron 
operators a and c Jl(T of Eq. (|16p have the same anticommutation relations as the corresponding electron operators cj a 
and Cj >a , respectively. Similarly, the local c fermion operators of Eq. (|61[) and six local spinon and 77-spinon operators 
of Eq. (|64|) have the same algebra as the corresponding unrotated spin-less and 77-spin-less fermion operators of Eq. 
(|57p and six local operators of Eq. (|42p , respectively. The former operators play a more general role in the physics of 
the model. The latter operators are a particular case of the former operators associated with the «> 1 limit physics. 
Hence for simplicity in the following we provide the algebra of the local c fermion operators of Eq. (|6ip and three 77s 
quasi-spin operators of Eq. (|66|) . The SU(2) algebra of the latter three operators fully determines that of the three 
local spinon operators and three local 77-spinon operators of Eq. (|64p , respectively. Given the unitarity of the the 
operator V, that of the corresponding unrotated local operators is similar. 

Straightforward manipulations based on Eqs. (l6i~j) and (|66|) lead to the following algebra for the c fermion operators, 

{/:, . ff,c} = h f ; {/:.. , /;,, j = {/„ , /,,., j = , (76) 

and the c fermion operators and the local rjs quasi-spin operators, 

If I 4>\ = \fi,o = Vic > S fJ = \ho , 4, J = > 1 = ±^3 , a = 77, s . (77) 

The SU (2) algebra obeyed by rotated local quasi-spin operators q^ where I = X3, ±, such that = q* 1 ± i q* 2 , and 
corresponding 77-spin (a = rf) and spin (a — s) operators Sj a is, 

[qf, qj,] = Sjj, 2 qf ; [qf^f] = T S^, qf , (78) 

and 

[ S 5,a,Sj>, a '} = hi' 5 <x,<*' 2 I I^a.^c'] = , a, a' = ? ? , s , (79) 

respectively. Moreover, one has that [gj, gj,] = and [A t<x , Sj>, a '] — where I = 0, ± and a, a' = 77, s. 

The c fermion and 77s quasi-spin operator algebras refer to the whole Hilbert space. On the other hand, those of the 
?7-spinon and spinon operators correspond to well-defined subspaces spanned by states whose value of the number 2S C 
of rotated-electron singly occupied sites is fixed. This assures that the value of the corresponding 77-spinon number 
= [N a — 2S C ] and spinon number M s = 2S C is fixed as well. 

IV. THREE EFFECTIVE LATTICES AND THE SPINON AND t^-SPINON TRANSFORMATION LAWS 

Here we introduce three effective lattices associated with our operator formulation. Furthermore, we investigate 
the consequences of the spinon and 77-spinon transformation laws under the electron - rotated-electron unitary trans- 
formation. 



A. Three effective lattices 

Within the present operator formulation, the rotated-electron occupancy configurations that generate the 4^" energy 
eigenstates have three degrees of freedom associated with the functions <i>uii\{ x \ > ■■■,x s NB ), 4>su(2\{ x ii x i^ ■■■^trv.o ), 

and 4>su(2)( x i^ ' xS m s -° )' res P ec ti ver yj m the rotated-electron amplitude general expression of Eqs. (|53|) and 

(|54p . Such rotated-electron occupancy configurations three degrees of freedom are within that formulation naturally 
described in terms of independent occupancy configurations of c fermions, 77-spinons, and spinons, respectively. (The 
relationship of the operators of the c fermions to those of the rotated electrons is provided in Eq. (|61l) . That of the 
operators of the 77-spinons and spinons to the rotated-electron operators is given in Eqs. (f66|) -(j67 | .) 

The direct relation in the u — > 00 limit of the rotated-electron amplitudes given in Eqs. (1551) and (1541 to those of 
the electrons derived in Ref. was clarified in Section III Bl We have found that in such rotated-electron amplitudes 
the spatial coordinates x\, ...,x s N „ refer to the rotated-electron singly occupied sites in the original lattice. Within 
our formulation, they correspond to the spatial coordinates Xx, xjsr c of the A^ c = 2^ c fermions in their c effective 
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lattice with N a sites. Moreover, the spatial coordinates xf,..., x d v , in such amplitudes correspond to the rotated- 

R, -1/2 

electron doubly occupied sites that are not generated from 77-spin-flip processes in an effective lattice containing 
only the rotated-electron doubly and unoccupied sites. Within our formulation they are found below to refer to the 
spatial coordinates xi, ...,x M bo of M^ _jy 2 = [Na ~~ 25c ~ 2S', ) ]/2 anti-bound 77-spinons of 77-spin projection —1/2 
in the 77-spin effective lattice with M v = [N a — 2S C ] sites. Indeed, the rotated-electron amplitudes of Eqs. (|53|) and 
(|54[) are independent of the spatial coordinates of the n v rotated-electron doubly occupied sites that are generated 
from 77-spin-flip processes. Those are found below to correspond to the M^ n _ 1 ^ 2 — unbound 77-spinons of 77-spin 
projection —1/2, which are invariant under the electron - rotated-electron unitary transformation. Finally, the spatial 

coordinates x\ , x s ^ Si0 refer to the down spins that are not generated from spin- flip processes in an effective 

s',-1/2 

lattice containing only rotated-electron singly occupied sites. They are found below to correspond to the spatial 
coordinates x\, ...,x M b<, ^ of M%°_y2 = [&c ~ S s ] bound spinons of spin projection —1/2 in the spin effective lattice 
with M s = 2S C sites. Again, the rotated-electron amplitudes are independent of the spatial coordinates of the n s 
down spins of rotated-electron singly occupied sites that are generated by spin-flip processes. Those are found below 
to correspond to the M^j/j = n s unbound spinons of spin projection —1/2, which are invariant under the electron - 
rotated-electron unitary transformation. 

For the Bethe states the numbers given in Eq. ([7]) vanish, n v = n s = 0. Thus all unbound 77-spinons in such states 
have 77-spin projection +1/2 and all unbound spinons have spin projection +1/2. This is equivalent to the absence of 
both rotated-electron doubly occupied sites generated from 77-spin-flip processes and down spins of rotated-electron 
singly occupied sites generated from spin-flip processes. 

The site order invariance occurring in the ID Hubbard model plays an important role in the occupancy configurations 
order of the sites of the spin, 77-spin, and c effective lattices. The occupancies of the sites of the spin effective lattice 
correspond to the spin degrees of freedom of those singly occupied by rotated electrons in the original lattice. A 
given spinon occupancy configuration of that effective lattice refers to the order of the spins of such rotated electrons, 
independently of the positions in the original lattice of the sites unoccupied and doubly occupied by rotated electrons. 
Similarly, the occupancies of the sites of the 77-spin effective lattice correspond to the 77-spin degrees of freedom of those 
unoccupied and doubly occupied by rotated electrons in the original lattice. A given 77-spinon occupancy configuration 
of such an effective lattice corresponds to the order of the rotated-electron unoccupied and doubly occupied sites, 
independently of the positions in the original lattice of the sites singly occupied by rotated electrons. 

Consistent, neither the spinon distribution of the functions <f>su(2){ x i > ■■■! a; tts,o ) n or the 77-spinon distribution 

associated with the relative distribution of the rotated-electron unoccupied and doubly occupied sites of the functions 
(f> 7 l u(2 Jxf,x d ,,...,x d ri:0 ) do change if the chain of rotated-electron singly occupied sites is "diluted" by rotated- 

electron unoccupied and doubly occupied sites, making possible also direct propagation for the rotated electrons. 
Such a propagation is described by independent occupancy configurations associated with the c fermion distribution 
of the function $uii) i x i > •••■> x n 3 )• This important property of the rotated-electron amplitudes of Eqs. (|53p and 
(|54")l is behind the independence of the c effective lattice, 77-spin effective lattice, and spin effective lattice occupancy 
configurations. 

Except for u» 1, such an independence of the rotated-electron occupancy configurations three degrees of freedom 
is lost if one uses descriptions based on electrons rather than rotated electrons, as defined in Section III Bl The c 
effective lattice is identical to the original lattice. Its sites are then labelled by the index j = 1, N a . For subspaces 
with fixed 2S C values, the representations of the c hidden U(l) symmetry algebra in the model global symmetry 
algebra are generated by C N a independent occupancy configurations of the N c = 2S C c effective lattice occupied sites 

relative to its Nf} = [N a — 2S C ] unoccupied sites. Here and throughout this paper the symbol C^ c refers to the usual 
mathematical formula, 



N n \ NJ 



C " : ~Wc)~ N c \(Na-N c )\' Na ~ Nc 

For a given energy eigenstate, the c fermion iV c = 2S C occupied and Nf? = [N a — 2S C ] unoccupied sites have the 
same spatial coordinates as the = 2S C sites singly occupied by rotated electrons and the = [N a — 2S C ] sites 

doubly occupied and unoccupied by rotated electrons, respectively. Hence = C N C = C N C . The c effective lattice 
occupancy configurations and corresponding c hidden U(l) symmetry algebra representations are insensitive to which 
of the original lattice = 2S C sites are singly occupied by spin-up or spin-down rotated electrons and which of that 
lattice N 7 l = [N a — 2S C ) sites are doubly occupied or unoccupied by rotated electrons. As found below, representations 
of the c hidden U(l) symmetry algebra may alternatively be generated by corresponding c band discrete momentum 
value occupancy configurations associated with the c fermion operators defined in Eq. ([55)) . 
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The spinon and 77-spinon distributions over the above-mentioned two sets of N R = 2S C and Nj> — [N a — 2S C ] sites 
of the original lattice, respectively, generate the spin SU(2) algebra and 77-spin SU{2) algebra representations. Those 
correspond to the spin and ?7-spin degrees of freedom of independent rotated-electron occupancy configurations of the 
sets of iVjj = 2S C sites and = [N a — 2S C ] sites, respectively, of the original lattice. Consistent, the M s = = 2S C 
spinons only "see" the = 2S C sites singly occupied by rotated electrons. The M v = ATS = [N a — 2S C ] ?7-spinons only 
"see" the ATS = [N a — 2S C ] sites doubly occupied and unoccupied by rotated electrons. The role of the corresponding 

C^ a — C^ c c fermion occupancy configurations is to store the information on the positions in the original lattice of 
the = 2S C sites singly occupied by rotated electrons relative to the N R = [N a — 2S C ] sites doubly occupied and 
unoccupied by rotated electrons. 

This is consistent with the spinon and 77-spinon distributions being independent. The distributions of M v — N R = 
[N a — 2S C ] 77-spinons and M s = = 2S C spinons use the sets of N R = [N a — 2S C ] rotated-electron doubly occupied 
plus unoccupied sites and = 2S C rotated-electron singly occupied sites as independent 77-spin and spin effective 
lattices, respectively. The number of sites of such effective lattices are in this paper often denoted by N ar) and N aa , 
respectively. They read, 

Na n = M v = Nl = [N a - 2S C ] ; N as = M s = N S R = 2S C , (80) 

and obey the sum-rule iV a = [N av + N as ]. The 77-spinon occupancy configurations refer to the operators Sj _ of Eq. 
(1571) . which act only onto the N a = N R = [N a — 2S C ] sites of the 77-spin effective lattice. As discussed in Section 
1111 A[ the spinon occupancy configurations correspond to the operators 5^ g given in the same equation, which act 
onto the N as = N R = 2S C sites of the spin effective lattice. This is assured by the operators (f — Tij lC ) and rij lC in 
their expressions provided in that equation, which play the role of projectors onto the 77-spin and spin effective lattice, 
respectively. 

Within the present N a 1 limit, the spin effective lattice spacing a s and 77-spin effective lattice spacing a v refers 
to the average spacing between the c effective lattice occupied sites and between such a lattice unoccupied sites, 
respectively, 

L N a 

a a = — = —a; a = 77 , s . (81) 

iv a(3 iy af3 

Within that limit, the spin effective lattice spacing a s (and 77-spin effective lattice spacing a v ) is a well-defined concept 
for finite values of the electronic density n (and hole concentration (1 — n)). Such spin and 77-spin effective lattices obey 
the physical requirement condition that in the (1 — n) — > and (1 — ri) — > ±1 limit, respectively, they equal the original 
lattice. Indeed, in the (1 — n) — ^ (and (1 — n) — > ±1) limit one has that N as — N a and the 77-spin effective lattice 
(and N a = N a and the spin effective lattice) does not exist. Consistent with the expression a a = L/N ap — \N a /N afl \a 
of Eq. (|81[) where a = 77, s, the 77-spin (and spin) effective lattice has the same length as the original lattice. The 
number of sites sum-rule [N a + N ag ] — N a holds. Thus the 77-spin and spin effective lattices have in general a number 
of sites N ari and N ae , respectively, smaller than that of the original lattice, N a . Their lattice spacings provided in Eq. 
([8lj) are in general larger than that of the original lattice. 

We recall that the present spinon description is normal-ordered relative to the S v = N a /2; S s = 0; 2S C = electron 
and rotated-electron vacuum in that such state is not populated by spinons. For it the number of 77-spinons is maximum 
and given by N a . Consistent, for it the spin effective lattice does not exist whereas the 77-spin effective lattice is identical 
to the original lattice. Note that the S v — N a /2; S s — 0; 2S C = electron and rotated-electron vacuum is the ?7-spin 
LWS of a multiplet tower of N a + 1 states. All such states have no spinons, no c fcrmions, N a 77-spinons, and N a 
c fermion holes. For such states electrons equal rotated electrons whose number is A^ = [S v + S^ 3 ] = 0, 2, 4, N a . 
All such electrons refer to A^/2 = [S n + S^ 3 ]/2 = 0, 1,2, ...,N a /2 spin-singlet on-site pairs. On the other hand, the 
77-spinon description is normal-ordered relative to the S v = 0; S s = 0; 2S C = N a absolute ground state, which is not 
populated by 77-spinons, yet the number of spinons is maximum, reading A^ a . Hence for it the 77-spin effective lattice 
does not exist, the spin effective lattice being identical to the original lattice. 

We denote the one-site states that describe the rotated-electron on-site occupancies by \j, lsu(2), h}- Here j = 
1, N a is the site index of the rotated-electron original lattice. The four values of the index lsu{2) — Q> t4-> ti 4- refer 
to the four rotated-electron occupancies associated with the 77s quasi-spin SU (2) symmetry degrees of freedom. Those 
separate into two types of occupancies: (i) The occupancies of the 77-spin effective lattice N ari = N R = [N a — 2S C ] sites 
associated with the 77-spin SU (2) symmetry degrees of freedom; (ii) Those of the spin effective lattice N as — Nf, = 2S C 
sites associated with the spin SU (2) symmetry degrees of freedom. Each on-site occupancy is described by the +1/2 77- 
spinon (l S u(2) = ©) and -1/2 77-spinon (l S u(2) =ti) and the +1/2 spinon (l S u(2) =t) and - 1 / 2 77-spinon (l S u(2) =1), 
respectively. The two index values l c = h c and l c — p c correspond to c effective lattice on-site occupancies. They refer 
to a c fermion unoccupied site and occupied site, respectively. Such c fermion occupancies are associated with the c 
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hidden U(l) symmetry degrees of freedom of the corresponding on-site rotated-electron occupancies. Hence l c = h c 
for l S u(2) = ©,tl and l c = p c for l su{2) =f,4- 

As given in the second expression of Eq. ([50]) . any of the 4 Na energy eigenstates \^i a .i & ,u) may be expressed as 
a superposition of a sub-set of states \u; xf, x s N s ; xf , x d Vt0 ; x\ , je*j, ) describing local rotated-electron 

R "r',-1/2 N R,-l/2 

occupancy configurations. The latter states may be rewritten in terms of spinon, 77-spinon, and c fcrmion occupancy 
configurations as, 

\*loc) = H\j,lsu(2)WcU))- (82) 

3=1 

Each of such subsets contains all states of the form provided in Eq. (|52")l with the same numbers of sites singly occupied 
by spin-up rotated electrons, singly occupied by spin-down rotated electrons, unoccupied by rotated electrons, and 
doubly occupied by rotated electrons. (The site variable j — l,...,N a in the argument of the indices lsu(2){j) — 
©,tj">t>4- an d l c (j) — h c ,p c of the local-rotated-electron states given in Eq. (j52")l accounts of course for their values 
dependence on the on-site index j.) 

Analysis of the expressions provided in Eq. (joTj) reveals that upon acting onto the one-site states \j, lsu(2) > h c ) (and 
\j,lsu(2)>Pc)), in addition to creating (and annihilating) one c fermion, a c fermion operator /J c (and fj >c ) removes 
one site from the 77-spin (and spin) effective lattice and adds one lattice to the spin (and ?7-spin) effective lattice as 
follows, 

/j, c li,0A) = \jA,Pc)\ fl c \j,n,h c ) = (-iy\j,i,p c ), 
4c|j,t,P c ) = li,0A); fj,c\j,UPo) = (-i) j \j,U,h c ) ■ (83) 

The spin and ?7-spin effective lattices are thus exotic. The number of their sites, N as = 2S C and N a — 2S£, 
respectively, varies by ±1 and =f1 upon creation/annihilation of one c fermion. This is why the ?7-spinon and spinon 
operator algebras refer to subspaces spanned by energy eigenstates with fixed 2S C values and thus well-defined ^-spin 
and spin effective lattice site numbers, N av = 25' 1 = [N a — 2S C ] and N as = 2S C , respectively. 

On the other hand, upon acting onto the one-site states \j, lsu{2) > 'c)j the 77s quasi-spin operators of Eq. (|66|) do not 
affect the c fermion occupancy. Thus they preserve the number of sites of both the spin and 77-spin effective lattices 
and transform such states as follows, 

it \j,Q,h c ) = St, (j,©, /i c ) = (-iy \j,U,h c ) , 
qf b', Q,h c ) = s]* n \j, Q,h c ) = -- I j, 0, he) , 

qj \j, U, h c ) = sj v \j, ti, h c ) = {-iy I j, 0, he) , 

qf li, U, K) = sf v I j, U, h c ) = \ \j, U, K) , 

q+ \jA,p c ) = st, \jA,Pc) = \i,UPc) , 
qf \j,t,p c ) = s* 3 s \j,t,Pc) = - g b'.t.Pc) , 
qj \j,i,Pc) = sj s \j,l,p c ) = \j,t,p c ) , 

li, i,p c ) = s% li, l,p c ) = \ |i, ;,p c > . (84) 

For simplicity, in Eqs. (l83l) and (|84|) we have not included the one-site state transformations that give zero. 

The result of the application onto any of the initial one-site states \j, Q,h c ), \j,tl,h c ), \j, t,Pc), and \j,i,p c ) of 
the four rotated-electron operators given in Eq. (|69[) is a trivial problem in terms of the fermionic rotated-electron 
operator algebra. In Appendix [5] it is confirmed that the application onto such states of the operational expressions 
of the four rotated-electron operators in terms of 77s quasi-spin operators and c fermion operators provided in Eq. 
(|69p leads to the same results. This confirms the faithful character of the c fermion representation and 77-spinon and 
spinon representations associated with the 77s quasi-spin algebra. 



B. Consequences of the spinon and ?7-spinon transformation laws 

Here we confirm that the transformation laws under the electron - rotated-electron unitary transformation of the 
?7-spinons (and spinons) as defined in this paper contain important physical information. 
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1. Unbound rj-spinons and unbound spinons 

A well-defined number of 77-spinons (and spinons) remains invariant under that unitary transformation. Since 
analysis of the interplay of the model global symmetry algebra representations with the BA solution reveals that they 
remain unbound, we call them unbound ±1/2 77-spinons (and unbound ±1/2 spinons). The invariance of the unbound 
77-spinons (and unbound spinons) stems from the off diagonal generators of the 77-spin (and spin) algebra, which flip 
their ?7-spin (and spin), commuting with the unitary operator V. This is why such generators have for u > the same 
expressions in terms of electron and rotated-electron operators, as given in Eq. (|41j) . Consistent with that invariance, 
the results of Ref. [22j confirm that creation of one unbound —1/2 7/-spinon (and one unbound —1/2 spinon) onto a 
LWS ground state leads to an increase of one in the number of lattice sites doubly occupied by both electrons and 
rotated electrons (and singly occupied by both spin-down electrons and spin-down rotated electrons). (Within the 
notation of that reference, unbound —1/2 77-spinons and unbound —1/2 spinons are called —1/2 Yang holons and 
— 1/2 HL spinons, respectively.) 

The values of the numbers M'^ r j s _ l , 2 of unbound ±1/2 77-spinons and M"^^ of unbound ±1/2 spinons are fully 
controlled by the 77-spin S v and 77-spin projection S* 3 = — § iV and spin S s and spin projection S* 3 — — N a , 
respectively, of the subspace or state under consideration as follows, 

M™ = [M™_ 1/2 + M™ +1/2 ]=2S a ; M«" ±1/2 = [S a T ; a = r,,s. (85) 

Thus the 77-spin S v , 77-spin projection S* 3 = — | N a , spin S s , and spin projection Sf 3 = — y N a of an energy eigenstate 
are fully determined by the 77-unbound spinons and unbound spinons occupancies. 

For Bethe states with finite spin S s and/or 77-spin S v all unbound spinons and/or unbound 77-spinons have spin up 
and 77-spin up, respectively. Application of the off-diagonal generators of the ?7-spin (a — rj) or spin (a = s) algebra 
provided in Eq. (1411) onto such S a > states, flips the 77-spin (a = 77) or spin (a — s) of the unbound 77-spinons or 
unbound spinons, respectively. The 2S S and 2S n different occupancies of the unbound spinons and unbound 77-spinons, 
respectively, originated by such flip processes generate energy eigenstates outside the BA solution subspace. 

The above invariance means that the 77-spin (a = 77) or spin (a = s) degrees of freedom of the rotated-electron 
occupancy configurations over a set of 2S a sites out of the whole set of N afj sites have for u > the same form in 
terms of electrons and rotated electrons. On the other hand, for finite values of u the corresponding c hidden U(l) 
symmetry degrees of freedom of such rotated-electron occupancy configurations are not in general invariant under 
the electron - rotated-electron unitary transformation. Such configurations are described by the c fermion occupancy 
configurations. The only exception is for spin fully polarized states. For those the on-site interaction does not play 
any role. 

The invariance of the 77-spin degrees of freedom of the above rotated-electron occupancy configurations implies 
that in each symmetry algebra representation there are exactly 2S V — [M^"Li/ 2 + ^^+1/2] sites such that M^_ x , 2 
sites are doubly occupied and M^ n +l ^ 2 sites are unoccupied both by electrons and rotated electrons. Furthermore, the 
invariance of the spin degrees of freedom of the sites singly occupied by rotated electrons implies that there are exactly 
2S S = [M™™ x , 2 + , 2 ] sites of the original lattice such that M™ n _ 1 , 2 sites are singly occupied both for spin-down 

electrons and spin-down rotated electrons and M"" x , 2 sites are singly occupied both for spin-up electrons and spin-up 
rotated electrons. Consistent, it is confirmed below that symmetry implies the number \N as — 2S S ] — [2S C — 2S S ] 
of sites of the spin effective lattice left over to be even and to refer to spin-singlet spinon occupancy configurations. 
Similarly, the number [N a — 2S V ] = [N a — 2S C — 2S V ] of sites of the 77-spin effective lattice left over is even and 
corresponds to 77-spin-singlet 77-spinon occupancy configurations. 

2. Anti-bound rj-spinons and bound spinons 

The rotated-electron occupancy configurations of the [Mr 1 + M^°] = [N a — 2S V — 2S S ] original-lattice sites left over 
are not invariant under the electron - rotated-electron unitary transformation. Indeed except in the u — > 00 limit, 
the occupancy-configuration 77-spin degrees of freedom of = [N a — 2S C — 2S V ] such sites and the spin degrees of 
freedom of the remaining — [2S C — 2S S ] sites are different in terms of rotated electrons and electrons, respectively. 
(The index bo used here in the site numbers and M h ° refers to bound, yet such numbers are shown below in 
Section [V Gl to refer to anti-bound 77-spinons and bound spinons, respectively.) 

Straightforward symmetry arguments given below reveal that both for rotated electrons and electrons these two 
sets of = [N a — 2S C — 2S V ] and = [2S C — 2S S ] sites refer to 77-spin-singlet and spin-singlet occupancy 

configurations, respectively. However, for finite u values such occupancy configurations preserve and do not preserve 
single and double occupancy for rotated electrons and electrons, respectively. The occupancy configurations of the 
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above remaining two sets of M™ = 2S V and M"" = 25*5 sites do preserve it both for rotated electrons and electrons. 
Thus that for finite u values electron double and single occupancy are not good quantum numbers results only from 
the occupancy configurations of [M%° + Mg°] = [N a — 2S V — 2S S ] sites out of N a sites. This is an important physical 
information brought about by the present operator formulation, which makes explicit the model global symmetry. 
The above site numbers, 

M\° = [N av -2S v ] = [N a -2S c -2S ri }, 

M b 8 ° = [N aa - 25 s ] = [2S C - 2S S ] , (86) 

are good quantum numbers. Their values are fully determined by those of the eigenvalue 2S C of the global c hidden 
U(l) symmetry generator and 77-spin S v or spin S s , respectively. Thus M°.° and are not independent quantum 
numbers. 

Although S„ and S s may be integer or half-integer numbers, the factor in the model global [SU(2) €5 SU(2) <g> 
XJ(X)\jZ\ symmetry implies that only 77-spin representations for which [S c + S^] and [S c — S 8 ] are integer numbers 
are allowed. Since the energy eigenstates refer to representations of the model global symmetry, that restriction 
applies as well to such states. Hence those symmetry requirements are behind the numbers = [N a — 2S C — 25^] 
and M'" = [2S C — 2S S ] of Eq. ([56]) being even numbers for all the 4^" energy eingenstates that span the model 
Hilbert space. They then refer to even numbers of ?7-spin-l/2 77-spinons and spin-1/2 spinons, respectively. This is a 
necessary condition for the corresponding rotated-electron and electron occupancy configurations of such two sets of 
= [N a — 2S C — 25,,] and M h ° = [2S C — 2S S ] sites being 77-spin neutral and spin neutral, respectively. 

The 77-spin S v and spin S s values and corresponding projection values S* 3 and 3 of an energy eigenstate are fully 
determined by the occupancy configurations of the unbound 77-spinons and unbound spinons, respectively. This is 
consistent with the rotated-electron and electron occupancy configurations of the two sets of = [N a — 2S C — 2S V ] 
and Mg° = [2S C — 2S S ] remaining sites being 77-spin and spin neutral, respectively. In Section [V] this is found to be 
consistent as well with the BA quantum-numbers occupancies. For rotated electrons M*°/2 = [N a /2 - S c - S v ] of 
such sites are unoccupied and M bo /2 — [N a /2 — S c — S^] sites are doubly occupied. Furthermore, M h ° j2 = [S c — S s ] 
sites are singly occupied by spin-up rotated electrons and M\° j2 = [S c — S s ] sites are singly occupied by spin-down 
rotated electrons. On the other hand, for electrons at finite u values the occupancy configurations of the set of 
[M^° + Mg] = [N a — 2S V — 2S S ] sites are 77-spin-neutral and spin neutral, respectively, yet do not preserve single and 
double occupancy. 

It follows from the known SU (2) operator algebra that application onto S a = states of the 77-spin (a = 77) or spin 
(a = s) off-diagonal generators provided in Eq. (|4Tj) gives zero. For such states, N af) — M a — M^°, so that there 
are no unbound 77-spinons (a = 77) or unbound spinons (a = s). On the other hand, application of these generators 
onto S a > states flips the 77-spin (a — 77) or spin (a = s) of a unbound 77-spinon or unbound spinon, respectively, 
but leaves invariant the rotated-electron occupancy configurations of the set of remaining sites. Such a set 
of (and M^°) sites refers to 77-spin-singlet (and spin-singlet) configurations. The energy eigenstates \$>i ot i &tU ) 
outside the BA solution subspace have exactly the same anti-bound 77-spinon and bound spinons configurations as 
the corresponding LWS |^j 01 j<^ ltt ) of the same 5/7(2) towers. This is consistent with the amplitude equalities given 
in Eqs. (|48|) and (|49|) . Due to the invariance under the electron - rotated-electron unitary transformation of the 
unbound 77-spinon and unbound spinon, such amplitudes are not sensitive to their occupancy configurations, except 
for the phase factor (-1)"" of Eqs. (gSJ) and @SJ where n v = M"™ 1/2 . 



C. Unbound spinon and unbound 77-spinon energies for both states inside and outside the BA solution 

subspace 

In Appendix [Bl the energy of the ±1/2 unbound 77-spinons is derived. This is achieved by combining symmetry 
with the chemical potential curve /1 = /-t(n) in Eq. (|C1[) of Appendix [C] being for n 7^ 1 and in the N a — > 00 limit a 
smooth, monotonous, and continuous decreasing function of the electronic density n. A similar analysis concerning 
the energy of the ±1/2 unbound spinons can be carried out. Again it combines symmetry with the magnetization 
curve 2fis H = 2[Lb H(m) in Eq. (IC6[) of Appendix [Cl being in the iV a — > 00 limit a smooth, monotonous, and 
continuous decreasing function of the spin density m. 

The derivation of the energy for creation onto a ground state of one ±1/2 unbound 77-spinon or one ±1/2 unbound 
spinon requires the clarification of the related problem of the occupancies of that state. On the other hand, the BA 
solution subspace is smaller than and contained in the ID Hubbard model full Hilbert space. Therefore, analysis of 
the energy spectra of the Bethe states alone is not enough to clarify the problem of the ground-state occupancies. In 
Appendix [B] the model global symmetry is combined with the B A solution to address that problem. Specifically, it is 
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found that ground states with electronic density n^l (and spin density m ^ 0) are LWSs and HWSs of the 77-spin 
algebra (and spin algebra) for n < 1 and n > 1 (and for m > and m < 0), respectively. Consistent, a n = 1 (and a 
rn = 0) ground state is both a LWS and HWS of that algebra. In particular, the fi = 0, n = 1, and to = absolute 
ground state of the Hamiltonian H symm of Eq. @ is both a LWS and HWS of the two SU (2) algebras. Hence it is 
both a 77-spin and spin singlet with numbers S v — S s — and 2S C = N a . This last result is fully consistent with two 



exact theorems proved in Ref. [35[, which hold as well for the Hubbard model on lattices other than the ID lattice. 

In Appendix [B] the energy £^±1/2 for creation onto a n ^ 1 ground state of one unbound 77-spinon of 77-spin 
projection ±1/2 is found to be, 

^,±1/2 = 2|/x| ; e^, T i/2 = 0, sgn{(l - n)}l = =Fl . (87) 

The relations given here hold both for finite spin densities in the range to € [0, n] for n £ [0, 1[ and to € [0, 2 — n] for 
rae]l,2]. 

The n = 1 ground state has no 77-spinons. Thus the derivation in Appendix |5] of the energy for creation onto 
that ground state of one unbound ±1/2 77-spinon relative to the zero-energy level of the Hamiltonian H of Eq. (TTJ) 
uses a method other than that for an initial ground state with electronic density n ^ 1. The result obtained in that 
Appendix is, 

e v ,±i/2 = (P° ± (j) at n = 1 and fi e [-/_t°, /_t°] . (88) 

Such a unbound 77-spinon energy is consistent with those provided in Eq. (|87|) for 77 7^ 1. This is confirmed upon 
changing the chemical potential within its n = 1 range, /i G (— fi ,fi ), in the energy expression of Eq. (1881) and 
accounting for it being given by fi — fi° for 77 = 1 + and fi = —fiP for 77 = 1 — . 

A similar analysis to that of Appendix IBI involving the spin SU(2) algebra leads to the following energies for creation 
of one ±1/2 unbound spinon onto a ground state with arbitrary spin density 777 G [0, 77] for 77 € [0, 1] and to € [0, 2 — 77] 
for ?7 € [1, 2]. Relative to the zero-energy level of the Hamiltonian of Eq. (flj in a finite magnetic field H, such a ±1/2 
unbound-spinon energies read, 

£ s ,±l/2 = 2/Xb \H\ ; e s , T i/ 2 =0, sgn{m}l = =pl . (89) 

The energy corresponding to a pair of unbound 77-spinons of opposite 77-spin projection is an important reference 
energy scale for our study of the energies of the anti-bound 77-spinon occupancy configurations. From the energies 
given in Eqs. (|87|) and (|88|) . one finds that relative to the zero-energy level of the Hamiltonian H of Eq. (fTJ), such an 
energy scale is given by, 

£i7.+l/2 + £ 77,-i/2 = 2 M for n + 1 , 

£r,.+i/2 + £17,-1/2 = 2 M° at n = l. (90) 

Also the energy scale associated with the energy of a pair of unbound spinons of opposite spin projection plays an 
important role in our studies. From the energies given in Eq. (|89l) we find that for all spin densities and relative to 
the zero-energy level of the Hamiltonian H of Eq. (TTJ), such an energy scale reads, 

e s ,+i/2 +£« -1/2 = 2/jb \H\ for any to. (91) 



D. The jy-spinon and spinon subspace and corresponding vacuum 

The numbers of ±1/2 77-spinons (a = 77) and ±1/2 spinons (a = s) are fully determined by the rotated-electron 
numbers of Eqs. (TI7j)) and (|lll) as follows, 

M Q , ±1/2 = Nl ±1/2 = [M™ ±1/2 + M b a °/2] = [S a T SS 8 + M h a °/2] = [M a /2 =p 5^] , 

M a = AJj = 7V acv = [M™ + M*°] - [25 a + M b a °] ; a = 77 , s . (92) 

The equalities given here do not imply that the ±1/2 77-spinons (a = 77) and ±1/2 spinons are the rotated-electron 
site occupancies with similar number values given in Eq. (|10p. They rather describe the 77-spin and spin degrees of 
freedom of the rotated-electron site occupancies associated with the numbers ±1 , 2 and ±1 / 2 , respectively. We 
call •q-spinon and spinon subspaces those spanned by states whose value of the number 2S C of rotated-electron singly 
occupied sites is fixed. In it the number values A ajj = [N a — 2S C ] and N aB = 2S C of sites of the 77-spin and spin 
effective lattices and thus M v = [N a — 25 c ] and M s — 2S C of 77-spinons and spinons, respectively, are fixed as well. 
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Within our LWS representation, the theory vacuum of each 77-spinon and spinon subspace is a 77-spin and spin 
LWS. Its maximum spin density m — (1 — x) — n is reached at a critical magnetic held H = —H c , whose absolute 
value is defined in Eq. (|C7|) of Appendix [C] Such a spin fully polarized state has no anti-bound 77-spinons and no 
bound spinons, — = 0. Its electronic density is in the range n £ [0, 1], so that the 2S C value of that state is 
2S C = N. It remains invariant under the electron - rotated-electron unitary transformation and reads, 

|0„ s ) = |0„; N av ) ® \0 s ;N a3 ) <g> \GS C ; 2S C ) . (93) 

Here the 77-spin SU{2) vacuum |0^;JV a ) associated with N a = [N a — N] unbound +1/2 77-spinons, the spin SU(2) 
vacuum |0 S ; N as ) with N as = N unbound +1/2 spinons, and the c band U(l) vacuum \GS C ; 2S C ) with N c = 2S C = N 
c fermions remain invariant under the electron - rotated-electron unitary transformation. Indeed, in this limiting case 
the c fermion operators are also invariant under the electron - rotated-electron unitary transformation. The above 
vacuum \GS C ; 2S C ) may be expressed as \\ q f^ c \GS c ; 0). 

The state \GS c ;0) in ]J q f\ tC \GS c ] 0) is the N c = 2S C = c fermion absolute vacuum. It refers to the S n = 
N a /2; S s = 0; 2S C = electron and rotated-electron vacuum, |0) = |0^; N a ) <S) \0 S ; 0) <g> \GS C ; 0). It is a limiting case of 
that provided in Eq. It corresponds to N ari = [N a — 2S C ] = N a and N as = 2S C — 0. Such a vacuum has no spin 

effective lattice, no spin degrees of freedom, no rotated electrons, no electrons, and no c fermions. Nonetheless within 
our formulation it is populated by M v unbound 77-spinons of 77-spin projection +1/2. It is a 77-spin fully polarized 
state with maximum hole concentration x = (1 — n) = 1 reached at chemical potential fi = — /i 1 = —[U + 4t]/2, Eq. 
(|C5P of Appendix [Cl It has maximum 77-spin value, S v — N a /2. The corresponding 2S C — subspace has dimension 
N a + 1. It is spanned by that vacuum and its 2S C — N a tower states. All such states have no spin effective lattice, 
no spin degrees of freedom, and no c fermions. They are populated by N = 2n v = 2M^ n _ l ^ 2 = 2,4, ...,2N a rotated 
electrons, which are invariant under the electron - rotated-electron unitary transformation. Hence for the 2S C = 
subspace such rotated electrons are electrons. 

Consistent with the c fermions one-site state transformation of Eq. (1831) , creation (and annihilation) of one c fermion 
is a process that changes the N c = 2S C value and thus drives the system into a subspace with a different vacuum, 
Eq. (|93p. As explicitly given in that equation, such a c fermion creation (and annihilation) process involves as well 
addition (and removal) of one site to (and from) the spin effective lattice and removal (and addition) of one site from 
(and to) the 77-spin effective lattice. It follows that from the point of view of the 77-spin and spin degrees of freedom, c 
fermion creation and annihilation processes correspond to a change of quantum system. Indeed, the 77-spin and spin 
lattices and corresponding number of sites change along with the quantum-system vacuum |0^ s ), Eq. (|93p . One can 
then say that for the 77-spinon and spinon representation there is a different quantum system for each eigenvalue 2S C 
of the generator of the global c hidden U(l) symmetry of Eqs. (|4*3"|) and (p5|) . (This is consistent with a quantum 
problem being defined both by a Hamiltonian and the Hilbert space its acts onto.) Specifically, the separate 77-spinon 
and spinon operator algebra representations refer to the quantum problem described by the ID Hubbard model in 
the 77-spinon and spinon subspace defined above. 

Finally, to illustrate the differences between the two types of occupancy configurations of unbound and bound 
spinons, which have finite spin S s and are spin neutral, respectively, we consider a state \ip). (Similar considerations 
apply to unbound and anti-unbound 77-spinons.) For simplicity but without loss in generality, we consider that it 
belongs to a subspace whose spin effective lattice has only four sites. Hence the state has N c = 2S C = 4c fermions, 
= [N a - 2S C ] = [N a - 4] c fermion holes, M s = 2S C = 4 spinons, and M n = [N a - 2S C ] = [N a - 4] 77-spinons, 
where N a 3> 1. Furthermore, we choose it to be a LWS with spin and spin projection S s = —S^ 3 = 1 and 77-spin 
and 77-spin projection S v = —S% 3 = [N a — 4]/2. Thus its numbers of unbound spinons and unbound 77-spinons are 
M "+i/2 = i 2S c ~ 2S s] = 2, M™_ 1/2 = and 1/2 = [N a - 4], M™_ 1/2 = 0, respectively. Furthermore, it follows 
from the above number values that the state has M™ = 2 bound spinons and = anti-bound 77-spinons. 

We denote by — 1,2,3,4) the four-site state that describes the spin-degrees of freedom of the state For 
simplicity, we have chosen the site index values j = 1, 2, 3, 4 to be those of the underlying spin effective lattice. Below 
we also indicate the c fermion occupancies associated with the c hidden U(l) symmetry degrees of freedom of the 
four sites. In terms of the one-site states |j, lsu(2)(j)i lc(j)) on the right-hand side of Eq. (|52"j) . we consider that our 
four-site state is given by, 

\tl>;j = 1,2,3,4) = i[|l,t,Pc)|2,|,p c )|3,t,Pc)|4,t,Pc> - |U,Pc)|2,t,jO|3,t,Pc>|4,t,Pc) 

+ |l,t,Po)|2,t,Pc)|3,t,Pc)|44,Pc) - |l,t,Pc)|2,t,Pc>|34,Pc>|4,t,Pc>] . (94) 
The identification of unbound and anti-bound spinons becomes clear if one writes such a state as a superposition of 
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two states. Each of them is a product of a spin-triplet two-spinon state by a spin-singlet two-spinon state, 

\tf>;j = 1,2,3,4) = -i={[|l,t,Pc)|2,|,p c > - \l,i,Pc)\2A, Pc )} x |3,t,Pc)|4,t,Pc)} 

+ -^={|l,t,Pc>]2,t )Pc > x [|3 > t,Pc)|4,4,Pc) -|3,4.,Pc)|4,t,Pc>]}. (95) 

The two unbound spinons are those of the spin-triplet two-spinon state. The two bound spinons are those of the 
spin-singlet two-spinon state. 

The present Sf 3 = — 1 state \ip) belongs to the subspace spanned by the Bethe states. The S'* 3 = and Sf 3 = 1 
states generated from it are outside the LWS BA solution subspace. Their corresponding four-site states, 

±St |V; j = 1, 2, 3, 4) ® |^) ® |Vc) ; \[St} 2 |V; j = 1, 2, 3, 4) ® |^> (8 |Vc) , (96) 

are given by, 

isj |V;j = 1,2,3,4) = 2^2{[|l J t,Po)|2,4.,Pc) - |l,|,Pc)|2,t,Pc)] x [|3,+,p c )|4,t,p c ) + [|3,t,p c )|4,|,p c )]} 

+ 2^2{[|l,4.,Pc)|2,t,Pc) + |l,t,Pc)|24,Pc)] x [|3,t,p c )|4,|,p c ) - |3,|, Pc )|4,t,p c )]}, (97) 

and 

- A [Stf |V;j = 1,2,3,4) = i={[|l,t,Pc)|24,Pc) - |14,p c )|2,t,p c >] x \3,i, Pc )\A,i,p c )} 

+ -j={\l,l, Pc )\2,l,p c ) x [|3,t,Pc>|44,Pc>-|34,Pc)|4,t,Pc)]}. (98) 

respectively. In the above equations, 5] is the spin off-diagonal generator given in Eq. (|4"Tj). States outside the LWS 
BA solution subspace differ from the corresponding states inside it only in the unbound spinons and/or unbound 
?7-spinon occupancies. Their c fermion, spin-neutral bound spinon, and ^-spin-neutral anti-bound ?7-spinon occupancy 
configurations are exactly the same. 

More generally, the spin effective lattice occupancy configurations of any state with M s = 2S C spinons and M" n = 
2^ unbound spinons such that 2S C > 2S S may be written as a sum of state configurations. Each of such configurations 
is a product of two states: A spin-i's M" ra -spinon state and a spin-singlet M* -spinon state where MJ™ = 2S S and 
M h ° = [2 S c — 2S S ], respectively. Alike for the simple configurations in Eqs. (pj)) , ([57]) . and the overall spin- 

effective lattice state is a superposition of configuration states where the positions of the unbound spinons and 
anti-bound spinons is different. However, their numbers are the same in all such states. 



V. RELATION OF THE c FERMION, SPINON, AND 77-SPINON OCCUPANCY CONFIGURATIONS 
TO THE QUANTUM NUMBERS OF THE BA SOLUTION 

The counting of the independent representations of the model global symmetry algebra in its subspaces with c 
fermion, spinon, and ry-spinon fixed numbers is used here to extract useful information on the relation of Takahasi's 
BA solution numbers to the objects of the present operator formulation. 



A. The thermodynamic BA equations in a notation suitable to the model global symmetry 

The notations used in Ref. [l| do not account for the relation of the BA quantum numbers that label the energy 
eigenstates to the model global SO(3) ® SO(3) <X> U{1) — [50(4) ® U{l)]/Zi symmetry algebra representations. Hence 
here we slightly change the notations used in that paper according to such representations. 

From analysis of the BA results of Ref. [l||, we fi n d l na t the number M^°/2 of anti-bound 77-spinon pairs considered 
in Section HVBI is exactly the number denoted by M' in that reference. Moreover, the quantity J^^Li n M n appearing 
in Eq. (2.13d) of that paper is the number M h s ° j2 of bound spinon pairs. It is confirmed below that the quantities 
denoted in Ref. [lj by M n (and M n ) are the number of specific types of spin-neutral spinon configurations (and 
77-spin- neutral 7/-spinon configurations) involving n bound-spinon pairs (and anti-bound- 7/-spinon pairs). In this paper 
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we denote the numbers M n and M n ' of that reference where n = l,...,oo by N su and N nu , where v = l,...,oo, 
respectively. The sets of numbers {N sl/ } and {N vu } are Bethe-state good quantum numbers, which obey the sum 
rules M b a ° = Y%Li 2v N av for a = r],s. 

The BA quantum numbers whose occupancy configurations generate the Bethe states include as well those associated 
with the degrees of freedom of the c hidden U(l) symmetry beyond 5*0(4), in the model global [SO (A) (g> f/(l)]/i?2 
symmetry. Consistent with the first expression given in Eq. f|B7[) of Appendix[B] the number of c fermions, N c , can be 
written as iV c = 2S C — [N — M*° — 2M"™ x , 2 ]. For the LWS Bethe states that span the BA solution subspace one has 
that 2M™_ 1 y 2 = 0. Hence for the BA solution the N c expression simplifies to Ac = 2S C = [N — M bo ]. The c fermion 
number N c = 2S C and thus the eigenvalue of the c hidden U(l) symmetry generator, 2S C , appears explicitly in the BA 
thermodynamic equations of Ref. where it is denoted by [N — 2M']. (Since M' is the notation of that reference for 
the number M bo /2 of anti-bound 77-spinon pairs, the number [N — 2M'] is indeed the number N c = 2S C = [ N — M bo ] 
of c fermions in the BA solution subspace.) 

The thermodynamic BA equations, Eqs. (2.12a)-(2.12c) of Ref. refer to subspaces spanned by Bethe states 
with fixed values for the set of numbers N c = 2S C and {N al/ } where a = 77, s and v = 1, 00. Such equations involve 
the integers or half-odd integers quantum numbers Ij, J Q ", and J£- (Here the index n = l,...,oo is that denoted in 
this paper by v = 1, ...,00.) We denote the indices a of Ref. 1], appearing in the numbers J a " and J™, by j and 
consider the following related quantum numbers, 

* = J^ I i'> J = h-,N a ; qj = ^-I^- j = l,...,N aa „, (99) 

where Ij = Ij, = 7j I/ , and = I™ , The BA quantum numbers are thus those given here in units of 2ir/N a . The 
latter are confirmed below to be discrete momentum values carried by quantum objects of our operator formulation. 
Although for different Bethe states one has that H, Ij" , and 1^ may be integers or half-odd integers, for each such 
a state the corresponding discrete momentum values qj have the usual momentum spacing, qj+\ — qj = 2ir/N a . 

The number values of the elementary objects whose occupancy configurations are associated with the BA solution 
quantum numbers and its present extension to the full Hilbert space obey the following equalities and sum rules 
controlled by the rotated-electron site occupancy numbers given in Eqs. (fTQ|) and 

N c = M s = N S R ; - M v = N V R , 

OO 

Y,^N av = N%-2S a , a = V ,s, 

v=\ 

00 

M a>±1/2 = M™ ±1/2 + Y J vN a „= JVg, ±1/a , a = 77, s . (100) 

This confirms that the BA solution accounts for the underlying rotated electrons, as defined in Section HT1 
The number of discrete momentum values, N aav , in Eq. f|99[) is for the BA av band given by, 

N aav = [N au + N* v ] , a = 77,s, i/ = l,...,oo. (101) 

Here N au and N^ v are the numbers of the av band occupied and unoccupied momentum values, respectively. From 
analysis of the results of Ref. [J , we find that the latter is given by, 

00 00 
tf*,= [25 a + 2 J2 {v' -v)N avl ] = [M^ n + 2 Y, {^-v)N au ,], a = V ,s. (102) 

Within the notation of that reference, N a /2 and N asv /2 are the numbers on the right-hand side of the two inequalities 
given just above its Eq. (2.13a). The corresponding BA quantum numbers I!j (and If v ) where j — 1, N a (and j = 
1, N aav ) are integers and half-odd integers for J2 av N av (and [N aav — 1]) even and odd, respectively. Furthermore, 
these numbers obey the inequality < [N a - l]/2 for J2 a u N ^ even and ~l N a ~ 2]/2 < IJ < N a /2 for J2au N ^ 
odd (and the inequality [If u \ < [N aau - l]/2 for both [N aau - 1] odd and even). 

In each subspace with fixed values for the set of numbers N c = 2S C and {N av \ where a = 77, s and v = 1, 00, 
the Bethe states are generated by the d c occupancy configurations of the N c c fermions in their c band with N a 
discrete momentum values, d VL , occupancy configurations of N VL/ filled 77^ band discrete momentum values over the 
N a v discrete momentum values for each r\v branch with finite N nll > occupancy, and d sv occupancy configurations 
of N sv filled 77^ band discrete momentum values over the N asv discrete momentum values for each sv branch with 
finite N su > occupancy. Here the dimensions d c , d^, and d su are given by, 

d c = C N ^ ; d^v — Cjy' ; d su — C N ^"^ , v = 1, 00 , (103) 
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respectively. It is confirmed below that the c band occupancy configurations generate representations of the c hidden 
U(l) symmetry algebra, the set of r\v bands occupancy configurations generate 77-spin-singlet representations of the rj- 
spin SU (2) symmetry algebra, and the set of sv bands occupancy configurations generate spin-singlet representations 
of the spin SU (2) symmetry algebra. 

Several physical quantities may be written in functional form in terms of the c, fjv, and sv band distribution 
momentum functions N c (qj), N vu {qj), and N su (qj), respectively. Those are eigenvalues of number operators N c (qj), 
N w (qj), and N s „(qj), respectively, whose expressions are given below in Section IV El The latter operators commute 
with the ID Hubbard model Hamiltonian. The Bcthe states are eigenstates of such operators. They have eigenvalues 
Np(<lj) = 1 f° r occupied momentum values and Np(qj) = for unoccupied momentum values where j3 = c,av, 
a = r),s, and v — 1, 00. 

That the quantum numbers qj are the c band and av band discrete momentum values is fully consistent with the 
following exact functional expression for the energy eigenstates total momentum, 

P = Y.1i Nc(qj) + EE H N ^i) + EE[ T - 1j] N vMi) + * M v,-i/2 > ( 104 ) 

j=l v=l j=l v=\ j = l 

which is additive in such values. From analysis of the results of Ref. one finds that for the Bethe states 

the momentum has the form given in Eq. (|104l) with M Vt _i/2 — M°°/2. By combining the BA solution with 
symmetry, we then straightforwardly confirm that the physical momentum expression given in that equation, where 
^j7,-i/2 = M%°/2 + M™ n _ 1 y 2 , is valid for all energy eigenstates. 

The c band discrete momentum values qj are the conjugate variables of the c effective lattice spatial coordinates j 
where j = 1, N a . Such discrete momentum values label corresponding c fermion operators given by, 

1 N " t 

i = 7fE e+,W 'i; /«.o=(4e) > ( 105 ) 



where we have used units of lattice spacing one. Consistent with the operator algebra of Eq. (|76[) . such c fermion 
operators obey the following anticommutation relations, 

{./;; |C . f qj , a = hr ; {./; , , 4, ,j = {/, , , /,• a = . m 

Within our notation the thermodynamic BA equations, Eqs. (2.12a)-(2.12c) of Ref. may be written in functional 
form as follows, 

kc(qj) = Vj-icrl^l^ ^fe')arctan' 

V=l j' = l 



rEE ^faOarctanf 8111 ^^^^^ ) ; j = l,...,N a , (107) 

a v=\ j' = l 



and 



^fe) = Qj iTr — 2^ ^cfe') arctan J —— J — 

j'=l 

— N a , 



atE E ^ 



Here, 



1/ = 1, ...,oo, j = 1, —,N aat/ , a = r),s. (108) 
k a u{q 3 ) = S a>v 2Re{arcsin(A^(9j) + iuu)} , 2/=l,...,oo, j = 1, ...,JV 0o „ , a = ?7,s, (109) 
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is the otv rapidity-momentum functional. It vanishes for a — s. The function Q v , V '{x) reads, 

u-1 

Qu,w( x ) = arctan(^-) + ^ 4arctan(^) | + (1 — ^,t/')|2 arctan 



1=1 



\v-v'\ 



2 arctan 



T—i)+ 4 arctan X - -)} . (110) 

v + t/J f— ' \\v-v'\+2lJ) y ' 



2=1 

I n 



The quantities A vl/ j = A rjl/ (qj) (and A su j = A su (qj)) are denoted in Ref. [lj by (and A a n ) where a = j = 
1, Na au and n = v = 1, 00. Within the notation of that reference, bound states of electron n-pairs are described 
by complex numbers A' and magnon bound states of n-pairs are described by complex numbers A lj. The quantities 
A^" (and A Q n ) in the above equations are the real part of such complex BA rapidities. 

Within the BA solution literature that followed the results of Ref. [lj], the latter quantities are associated with 
excitations often called strings of length n = v = l,...,oo [30]. In Sections 4.2.1 and 4.2.2 of Ref. [3(| those were 
studied for two very simple limiting states, for which electrons may be approximated by rotated electrons. In the 
thermodynamic limit the two electrons with opposite spin projections in the 77-spin v — 1 string state considered in 
the former section occupy the same site. The v = 1 string involves an unoccupied site, in addition to the spin-singlet 
electron pair. The spin v = 2 string state considered in the latter section involves the spin degrees of freedom of 
the two down-spin electrons and of two out of its N — 2 up-spin electrons, rather than only those of the two former 
electrons. This is in spite of the corresponding wave functions, which are of the general form given in Eq. (|19[) . only 
depending on the spatial coordinates of the spin-down and spin-up electrons and two down-spin electrons, respectively, 
as discussed in Section iBBl 

We denote by A^faj) the complex BA rapidity whose real part A ail {qj) dependence on the av bands momentum 
values qj is obtained by solution of Eqs. (I107p - (|109[) . Its imaginary part is in the N a —> 00 limit independent of 
such momentum values. It depends solely on the on-site repulsion in units of At, u — U/4±, and number of pairs, 
v = 1, 00, 

A L/fe) = A av {qj) + v- 27) u, j=l,...,f, j = 1, N aau , v = l,...,oo, a = rj ) s. (Ill) 

For v > 1 the v imaginary BA rapidity terms i(l + v — 27) u where 7=1, v are within our operator formulation 
related to the attraction between the corresponding v spin-neutral spinon pairs (a = s) or ?7-spin-neutral 77-spinon pairs 
(a = 77). Such v pairs have the same momentum qj. Except for qj = and qj = ±g Q „, where q av = Tr[N aatl — 1]/L, 
they move coherently within the same 2^-spinon (a = s) or 2i/-fy-spinon (a = rj) composite object. (At both qj = 
and qj = ±g Q „ their velocity is found below to vanish.) That attraction is behind the formation of those of such 
composite objects that contain v > 1 pairs. Hence it is to be distinguished from the two-spinon binding and two-77- 
spinon anti-binding of each pair. That for v = 1 and both a — 77, s the BA rapidity A^^qj) = A a \(qj) of Eq. (jllip 
is real follows from it referring to a single pre-formed spin-neutral spinon pair (a = s) or 77-spin-neutral 77-spinon pair 
(a = 77). In the second paper it is confirmed that there are no bound states of such 2i/-spinon (a = s) or 2v- 77-spinon 
(a = 77) composite objects. 

The c band and av bands momentum values qj appearing in Eqs. (|107[) - (|109p and (jllip . which are defined in Eq. 
(|99p , are the actual quantum numbers of the B A solution. The role of the corresponding c band rapidity momentum 
functional k c (qj), r\v band rapidity functional A vu (qj), and sv band rapidity functional A su (qj) is to explicit the 
dependence of the Bethe-states energy spectrum on such quantum numbers occupancy configurations. As confirmed 
below, that energy spectrum depends on the momentum occupancy configurations of the c band and av bands through 
such functionals. 

Combining the BA results of Ref. |l| with the model global symmetry algebra representations, one confirms that in 
the present N a — > 00 limit the BA solution performs a partition of the M^°-7?-spinon 77-spin-singlet (and Mg -spinon 
spin-singlet) configuration into a set of smaller independent 77-spin-singlet (and spin-singlet) configurations. Their 
number is, 

00 

B a = ^N av , a = 77,s. (112) 

v=l 

Each of such 77^ and sv configurations involves 2v anti-bound 77-spinons and 2v bound spinons, respectively. Here 
v = 1, gives the corresponding number of 77-spin-singlet anti-bound- 77-spinon pairs and spin-singlet bound-spinon 
pairs, respectively. For v > 1 the attraction between the v pairs of each configuration is associated with the imaginary 
part of the BA rapidity given in Eq. (jllip . Whether the BA quantum numbers /| of Eq. ([M)) are integer numbers or 
half-odd integer numbers is fully determined by the values of the numbers B v and B s defined here. The BA quantum 
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numbers I? in that equation where f3 = c, av, a = 77, s, and v = 1, oo are determined by the following boundary 
conditions, 

if = 0,±1,±2,... for G p even, 

= ±1/2, ±3/2, ±5/2,... for Gp odd, (113) 

where 

Gp = 6p >c Ba+ 6 H^[Nc + N av - 1] , (114) 

a=r),s 

and B a is defined in Eq. (I112p . Such values are equivalent to those provided in the text below Eq. (|99l) . From 
manipulations of the expression provided in Eq. (|102p . one straightforwardly confirms that N aau = [N av + N£ v ] is an 
odd and even integer number provided that [N c + N au ] is as well an odd and even integer number, respectively.) 

Consistent with both the relation between the c fermion and rotated-electron occupancy configurations and the 
2i/-spinon (a — s) or 2^-77-spinon (a = rj) composite character of the above av objects, the following exact sum-rules 
involving summations over the BA quantum numbers qj in Eqs. (|107[) - (|109j) hold, 

Na 



N c ( qj ) = N C = 2S C , 
[1 - N c ( qj )} = = [N a - 2S C 



N a 



3 = 1 



Y J2 2 " N vAQj) = M b n ° = [N a - 2S C - 2S V ] , 

u=l 3 = 1 

00 N. s „ 

£ 2vN sv ( qj ) - M b ° = [2S C - 2S V ] . (115) 

u=l j=l 

From the use of Eq. (|102j) . one finds that the number N% v provided in that equation simplifies for v = 1. It 
can be expressed only in terms of the number M a of 77-spinons (a — i]) or spinons (a — s) and the corresponding- 
configuration number B a given in Eq. (|112[) as follows, 

N* 1 = [M a -2B a ] ] a = r),s. (116) 

It is confirmed below that the length n = v = 1, ...,00 of the BA string A' n (and A") excitations is indeed the 
number of ?7-spinon pairs (and spinon pairs) in each of the independent ?7-spin-singlet (and spin-singlet) configurations. 
To start with, we call tjv (and sv) composite particle each of such 2^-r/-spinon (and 2^-spinon) configurations. For 
each av momentum band there is one branch of av composite particles. The indices a = T),s and v = 1, 00 label 
the distinct branch types of such av composite particles. 

Within the N a ~ > 00 limit that the thermodynamic BA equations (|107[) - (|109p refer to, it is often convenient to 
replace the discrete momentum values qj of Eq. such that qj+i — qj — 2ir/N a: by corresponding continuous 

momentum variables q. Those belong to domains q € [— Qp,+Qp] where /3 — c,av whose limiting absolute values qp 
read, 

Qc = 7T, 

q a » = ^r{N aav -l) for N aai/ odd, 

= ^rN aav for N aau even. (117) 

For the j3 = av branches the /3 discrete-momentum values distribution is symmetrical. The q av expression provided 
in this equation has no 1/N a corrections. Accounting for such c band corrections leads to q e [q~ ,q^] where the 
limiting momentum values qf are provided in Eqs. (B.15)-(B.17) of Ref. [2l| . 

The thermodynamic BA equations, Eqs. (|107|) - (|109[) . refer only to Bethe states. They can be extended to non-LWSs 
of form given in Eq. (|13l) . This is fulfilled by formally setting some of the rapidities A vl/ and A sl , in such equations 
equal to infinity [39j,|40(. For example, Eqs. (3.23b) and (3.24b) of Ref. [39j describe a 77-spin non-LWS state with 
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numbers S v = 1 and S* 3 = 0. Within our notation, that state rapidity value, A r) i = oo, refers to a 77-spin-triplet 
pair of unbound 77-spinons of opposite 77-spin projection, rather than to one 77-spin-singlet two-spinon composite 77I 
particle. On the other hand, Eqs. (3.23a) and (3.24a) of Ref. [1^] describe a Bethe state with numbers S v = S™ 3 = 0. 
Its rapidity A^i(0) refers to such a 77-spin-singlet two-77-spinon composite 77 1 particle. 

Alternatively, within the present operator formulation the non-LWSs are generated from the corresponding Bethe 
states by application onto those of the off-diagonal generators of the 77-spin or spin SU (2) symmetry algebras, as given 
in Eq. (p]). 



B. Consistency of the representations of the model global symmetry algebra with the BA solution and our 

formulation objects 

Here we confirm that the composite r\v particles and composite sv particles are 77-spin-neutral and spin-neutral 
objects made out of 2v 77-spinons and 2v spinons, respectively, where v = l,...,oo. That such 2v 77-spinons and 2v 
spinons are in anti-bound and bound configuration states, respectively, is confirmed in Section [VGI 



1. Three important subspaces 

The quantum-number occupancy configurations in the thermodynamic BA equations, Eqs. (|107[) - (|109j) . generate a 
well-defined set of Bethe states. Those span the subspace with given fixed values of 5*^, S s , and the sets of numbers 
N c = 2S C and {N au } where a = 77,5 and v = l,...,oo. We call it Takahashi subspace. The extended Takahashi 
subspaces are spanned by all energy eigenstates with given fixed values of S n , S s , and the sets of numbers N c = 2S C 
and {N au } where a = 77, s and v = 1, ...,00. That includes as well the non-LWSs states generated from the Bethe 
states that span the Takahashi subspaces. The dimension of an extended Takahashi subspace then reads, 

d ets = d c x J[ (2S a + l)d av . (118) 

a— r),s 

The dimensions d c and d av appearing here are given in Eq. (|103|) and the S v and S s values are fixed. Those read, 

a = 77, s, 2S C fixed. (119) 



s - 1 - 



OO 



M a -^2vN a „ 

v=l 



We divide the Hilbert space of the ID Hubbard model, Eq. (JTJ, in a set of larger subspaces that contain the 
extended Takahashi subspaces. They are spanned by the states with fixed values of the numbers 2S C , S v , and S s . 
Since the number values of unbound 77-spinons, M™ = 2S V , and unbound spinons, M"" = 2S S , are fixed, we call such 
subspaces unbound rj-spinon and unbound spinon subspaces. Their dimension is determined by the number J\f(S s , M s ) 
of M s -spinon spin SU (2) symmetry algebra state representations, N{S Vl M v ) of M^-77-spinon 77-spin SU(2) symmetry 
algebra state representations, and d c of 2S' c -c-fermion c hidden U{1) symmetry algebra state representations at fixed 
values of the 77-spin S v and spin S s . (We recall that M s = 2S C and M v = 2S'^ = [N a — 2S C ].) It reads, 

d urjS =d c x Af(S v ,M v ) x M(S S , M s ) . (120) 

Here the dimension d c is that provided in Eq. (|103[) and the two a = 77, s dimensions Af(S a , M a ) are given by, 

Af(S a ,M a ) = (2S a + 1) {c^: /2 - S « - C^l' 2 ' 8 "' 1 } , a = r,,s. (121) 

The sets of numbers {N al/ } where a = 77, s and v = 1, 00 of the states that span the smaller extended Takahashi 
subspaces contained in each of the unbound 77-spinon and unbound spinon subspaces have fixed values, which obey 
the following two sum rules, 

00 

M h ° = \M a -Ml n \ = [M a -2S a ] =2j2"N ai/ ; a = rj,s. (122) 



On the other hand, the unbound 77-spinon and unbound spinon subspaces considered here are contained in the 
larger 77-spinon and spinon subspaces introduced in Section IIVDI Those are spanned by the energy eigenstates with 
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fixed value for the c fermion number N c — 2S C , spinon number M s = 2S C , and 77-spinon number M v = [N a — 2S C ]. 
The dimension of the ry-spinon and spinon subspaces is given by, 

I^- 2 ^] 2S < [1 + f_r\l2S a +2S e ]] 

d vss = d nss (2S c ) = J2 J2 d - x II 2 l M(S a ,M a ). (123) 

2S,,=0 2S„=Q a=n,s 

The X^2S anc ^ summations in this equation run over all consecutive integer numbers 2S V =0,1,..., [N a — 2S C ] 

and 2S S = 0, 1, 2S C , respectively. The role of the phase factor (_i)[ 2S c*+ 2S, c] j s ^ se i ec t only the representations of 
the model global symmetry algebra, for which both [2S V + 2S C ] and [2S S + 2S C ] are integer numbers. 

Finally, alike in Ref. [19( for general bipartite lattices, summation over the dimensions of all ?7-spinon and spinon 
subspaces gives indeed the correct dimension of the ID Hubbard model Hilbert space, 

N a 

d H ii= J2 dr, ss (2S c ) = 4 N ° . (124) 

2S C =0 

Here the J^ 2S summation runs over all consecutive integer numbers 2S C = 0, 1, N a . 

A side technical question clarified in the introductory discussions of Section HI is how did the authors of Ref. 
[20j reach the correct number, A Na , of energy eigenstates in spite of assuming that the model global symmetry was 
50(4) = [SU{2)®SU{2)]/Zl rather than SO (3) <8 SO (3) <g> U(l) = [504) ® U{l)]/Z 2 1 As mentioned in SectionH the 
answer is simple. To generate the energy eigenstates outside the BA solution, they used the off-diagonal generators 
of the two non-Abelian SU{2) symmetries. On the other hand, to enumerate the corresponding LWSs l^^o^u) they 
used the BA solution, rather than symmetry. Since that solution accounts for the extra c hidden U(l) symmetry 
algebra representations, they reached the correct number, A Na , of energy eigenstates. 



2. Confirmation of the av particle composition 

The dimension d c given in Eq. (|103p equals the number of occupancy configurations of the N c = 2S C c fermions 
in their c band with N a discrete momentum values. Their allowed occupancies are 1 and 0. On the other hand, 
the oa> composite particles that emerge from the BA solution also carry discrete momentum values qj whose allowed 
occupancies are 1 and 0. Thus the number of occupancy configurations of the set of N au av composite particles in 
their av band with N aav discrete momentum values is that given in Eq. (I103p . d al/ = C^ av ■ It trivially follows from 
the ?7-spin (and spin) SU(2) symmetry algebra that its Af(S v , M v ) (and Af(S s , M s )) representations involve a number 
of 77-spin-singlet (and spin-singlet) representations given by, 

K r f Q at \ N{S ai M a ) n M a /2-S a n M a /2-S a -l / 10 r\ 

N singlet (S a , M a ) = + = C M J - C M J ,a = r),s. (125) 

If the r\v composite particles (and sv composite particles) are 7/-spin neutral (and spin neutral) composite objects 
of 2v 77-spinons (and 2v spinons) and, as given in Eq. (|122l) . all corresponding 77-spinons (a = rj) and spinons 
(a = s) are part of them, the values of their numbers N aav must exactly obey the following equality for all unbound 
?7-spinon and unbound spinon subspaces, 

00 

E IP 

{N m } v=\ 

Here the {N aiy } summation runs over the whole set v — l,...,oo of N au numbers that owing to the conservation 
of the a = T), s numbers and M h s ° exactly obey the subspace sum rules given in Eq. (|122j) . Hence such a 
summation is over all smaller extended Takahashi subspaces contained in the larger unbound 77-spinon and unbound 
spinon subspaces. The number Af S ingiet(S a , M a ) is that given in Eq. (|125[) . 

The general dependence of the number N{£ u of unoccupied av band discrete momentum values on the S a value and on 
those of the set of numbers {N av >} where v 1 = 1, 00 is defined by imposing that Eq. (I126|) is simultaneously obeyed 
in all subspaces. The unique mathematical solution of the problem is indeed that provided by the thermodynamic 
BA solution of Ref. Equation (|126[) is simultaneously obeyed in all subspaces provided that the number N aav — 
[N a u + A^„] of av band discrete momentum values and the corresponding number N^ v of discrete momentum values 
unoccupied by av composite particles are those provided in Eqs. (|101[) and (|102p . respectively. This confirms that 
the BA solution explicitly accounts for the model global symmetry algebra and corresponding representations. 



00 00 

= II C n:1 = E II C nZ„ =Kinglet(S a ,M a ), a^ V ,s. (126) 

{«„„} v=X {N a ,A v=\ 
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Moreover, the equality of Eq. (|126l) confirms that the r\v band (and sv band) occupancy configurations refer to 77- 
spin-neutral (and spin-neutral) 77-spinon (and spinon) occupancy configurations. Hence they generate representations 
of the U{\) symmetry sub-algebra within the 77-spin (and spin) SU(2) symmetry algebra. On the other hand, the 
unbound 77-spinon (and unbound spinon) occupancies generate the remaining representations the 77-spin (and spin) 
SU(2) symmetry algebra. 

Within chromodynamics, the quarks have color but all quark-composite physical particles are color- neutral [5l| . 
Here the M*° 77-spinons (and spinons) that are not invariant under the electron - rotated-electron unitary 

transformation have 77-spin 1/2 (and spin 1/2) but the corresponding 2v- 77-spinon r]v composite particles (and 2v- 
spinon sv composite particles) are 77-spin neutral (and spin neutral). 



C. Subspace dependent av composite particles statistics 

A quantum problem is defined by a Hamiltonian and the subspace it acts onto. Hence the ID Hubbard model in 
each of the three subspaces considered above in Section [VB] corresponds to three distinct quantum problems. Here we 
show that for each of such three quantum problems the av composite particles have a different statistics, as classified 
according to the generalized principle of Ref . [31j . 

In all the three subspaces the number of c fermions is fixed and given by the eigenvalue of the c hidden U{\) 
symmetry generator, iV c = 2S C - Since the 77-spinon and spinon operator algebra refers to the 77-spinon and spinon 
subspace, the dimension relevant for deriving the statistical interaction g a v,a'v' — $a,a' gv.y reads, 

M av = 1 + N* u . (127) 

The statistical interaction is related to the changes SAf al/ of the dimension given here as a result of the deviations 
5N al ,i in the av' composite particle numbers where v' — 1, 00. The statistical interaction does not connect r]v 
and sv' branches, consistently with those being associated with independent 77-spin and spin degrees of freedom, 
respectively. 

The statistical interaction g av , a > v > — <5 a ,a' guy is defined by the equation, 

00 

SATau = ~J2 3v.y SN av , . (128) 



1. The av composite particles statistics in the rj-spinon and spinon subspaces 

In the 77-spinon and spinon subspaces the c fermion number N c = 2S C , spinon number M s = 2S C , and 77-spinon 
number M v = [N a — 2S C ] have fixed values. The conservation of M a for a = 77, s implies that the set of deviations 
and {8N au } where v = 1, 00 must obey the sum rule, 

00 

SM a = 8Ml n + ^2v 6N au = , a = 77, s . (129) 
i/=i 

In order to classify the statistics of the av composite particles according to the generalized principle of Ref. [3l[ , 
it is convenient to rewrite the number of av band holes given in Eq. (|102[) as follows, 

00 

N^^[M a -^{v + v' -\v-v'\)N av ,], M a fixed, a = r],s. (130) 
i/'=i 

From the use of Eqs. (|127|) - (1130[) one readily finds the following av composite particle statistical interaction g v y 
for the quantum problem defined by the ID Hubbard model in the 77-spinon and spinon subspace, 

guy = {v + y' — \y — , v, v' = l, 00 . (131) 

Hence for that quantum problem the av composite particles are neither fermions (g u y — §vy ) nor bosons {g v y = 0). 
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2. The av composite particles statistics in the unbound n-spmon and unbound spmon subspaces 

In the unbound 77-spinon and unbound spinon subspaces the 77-spin S v and spin S s values are fixed, along with 
those of the c fermion number iV c = 2S C , spinon number M s = 2S C , and 77-spinon number = [N a — 2S C ]. The 
conservation of such numbers implies that the set of deviations {8N al/ } where v = 1, 00 must obey the sum rule, 

00 

^2vSN av = 0. (132) 

!/=l 

In order to classify the statistics of the av composite particles according to the generalized principle, the suitable 
expression for the number of av band holes is that given in Eq. (|102[) . By combining such an expression with 
Eqs. (|127p and (|128|) we readily find the following av composite particle statistical interaction g v y for the quantum 
problem defined by the ID Hubbard model in the unbound 77-spinon and unbound spinon subspace, 

gu,u' — 0, v = 1, ...,oo, v' < V , 

= -2(v'-v), z/ = l,...,oo, v'>v. (133) 

Again for that quantum problem the av composite particles are neither fermions nor bosons. 



3. The av composite particles statistics in the SN aL , = subspaces and extended Takahashi subspaces 

A SN^u = subspace is a subspace of a unbound 77-spinon and unbound spinon subspace for which the number 
N^ v vanishes for a specific av branch and 8N av i = for all v' > v branches. The conservation of N^ v and {N al/ >} 
for v' > v, along with that of the unbound 77-spinon and unbound spinon subspace numbers, implies that the set of 
deviations {SN au >} where v' — 1, v must obey the sum rule, 

V 

v 8N av > = . (134) 

v' = \ 

Therefore, for the quantum problem defined by the ID Hubbard model in such a subspace the av composite particle 
statistical interaction g v y is a particular case of that provided in Eq. f|l 33[) given by, 

9»y = • (135) 

It then follows that for such quantum problem the av composite particles are bosons. 

That the concerning the site occupancy configurations of the av effective lattice the av composite particles behave 
actually as hard-core bosons and thus the corresponding av operators generated from theirs by a Jordan- Wigner 
transformation obey an anticommuting algebra is consistent with the following properties: 

- An anticommuting algebra for both the c fermion operators of Eq. (|105[) and the av operators obtained from those 
of the hard-core bosonic av composite particles generates from the electron and rotated-electron vacuum a faithful 
representation of the energy eigenstates that span the 8N^ V = subspaces and extended Takahashi subspaces. This 
applies to all the model 4 Na energy eigenstates, since each of them belongs to a 8N^ V — subspace and extended 
Takahashi subspace associated with a well-defined initial ground state. 

- Such an energy eigenstate representation involves the av band discrete momentum values q,j — [2ir/N a ] of 
Eq. (1991) in the momentum general expression given in Eq. (|104p provided that they are carried by av objects whose 
anticommuting operators are generated from those of av composite particles. This is consistent with within the exact 
BA solution the set of av quantum numbers I" u of Eq. (|113l) for f3 — av and j = 1, N aav with N aau defined by 
Eqs. (|101j) and (|102l) . which are the above discrete momentum values qj = [2n/N a ] in units of [27r/A a ], having 
for all energy eigenstates Pauli-like allowed occupancies and 1, respectively. This is alike the occupancies and 1 of 
the BA c quantum numbers 7| associated with the j3 — c band discrete momentum values qj = I? [27r/iV a ], which are 
carried by the c fermions whose operators are given in Eq. (|105[) and obey the anticommuting algebra of Eq. (|106[) . 

- The discrete momenta qj = {2ir/N a } I? provided in Eq. (|99[) with the quantum numbers I? values as given in Eq. 
(| 1 1 3[) by the exact BA solution are also exactly those that emerge from the boundary conditions brought about by 
the Jordan- Wigner transformation that generates av fermions from the hard-core bosonic av composite particles. 
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- The dynamical and spectral function expressions derived for general u values within the PDT of Refs. (8l-lll| 
agree in the «> 1 limit with those of the method of Ref. [52| provided that alike the c fermion operators of Eqs. 
()105|) and ()106|) . the si operators obtained from those of the si composite particles obey anticommutation relations. 
Within the scattering theory introduced in the second paper, the c and si pseudofermion operators generated from the 
anticommuting c and si fermion operators, respectively, have modified phase-shift dependent exotic anticommutators 
that control the PDT quantum overlaps of the dynamical and spectral function matrix elements between the ground 
state and the excited states. (Within the method of Ref. [EH, which is valid for u ^ 1, the charge degrees of freedom 
are described by spinless fermions that refer to the c fermions, whereas the spin degrees of freedom are described by 
the wave function of the ID spin-1/2 isotropic Heisenberg model.) 

Hence we rename the the av composite particles, av bond particles. The av effective lattices are defined below. 
Their introduction simplifies the description of some particular types of the spin (a = s) or ?7-spin (a = rf) effective 
lattice occupancy configurations. 

For the particular case of the al branches, Eq. (| 134)) is only fulfilled for an extended Takahashi subspace. An 
extended Takahashi subspace is a SN^ V — subspaces for all its av bands with finite occupancy. For it the same 
numbers as for a unbound ry-spinon and unbound spinon subspace are fixed, along with those of the sets of av 
composite bound-particle numbers {N av } where a — rj,s and v = l,...,oo. The Bethe states contained in such 
subspaces span the Takahashi subspaces associated with the thermodynamic BA equations, Eqs. (|107|) - (|109p . It then 
follows from Eq. (|102j) that numbers of av band holes are fixed as well. In general we consider that the the 6N^ V = 
subspaces are extended Takahashi subspaces. 



D. The av fermions emerging from the hard-core av bond particles 

Above we have established that the spinons (and ry-spinons) that are not invariant under the electron rotated- 
electron transformation are part of spin (and 77-spin) neutral 2zy-spinon (and 2z/-?7-spinon) composite sv (and r]v) 
bond particles. 

The concept of a av effective lattice is well-defined in the extended Takahashi subspaces. For the ID Hubbard 
model in such subspaces, the numbers N aav , N av , and N£ u in Eqs. (|101[) and (|102l) give as well the numbers of 
the av effective lattice sites, occupied sites, and unoccupied sites, respectively. An occupied site of the av effective 
lattice corresponds to 2v = 2,4, ...,00 sites of the ?7-spin (a — rf) or spin (a = s) effective lattice. Consistent with 
the expressions provided in Eqs. (jlOll) and (|102l) . the 2S V (and 2S S ) sites of the Ty-spin (and spin) effective lattice 
occupied by unbound ?7-spinons (and unbound spinons) and 2(i/ — v) sites out of the 2v' sites of that lattice occupied 
by each hard-core r/v' (and sv') bond particle such that v' > v play the role of unoccupied sites of the rjv (and sv) 
effective lattice. 

Provided that the ratio N aav /N a remains finite as N a ~ > 00, the av effective lattices can be represented by ID 
lattices with spacing, 

T N a N a0 



a = ^7 a a ; N aau > 1 , (136) 

where v = 1, 00 and a — 77, s. The arguments behind the ID av effective lattice average distance a av between its 
sites playing the role of lattice spacing are similar to those used in Section llV Al for the spacing of the 77-spin and spin 
effective lattices. 

The traditional operational formalism of second quantization can be applied to the av bond particles for transitions 
within the extended Takahashi subspaces. Indeed, for such subspaces the dimension Af au = [1 + N„ v ] of Eq. (|127[) 
remains constant. In terms of the occupancies of the av band discrete momentum values, such transitions involve 
particle-hole processes. As discussed below, the av bond particle creation or annihilation operators also play a role 
in transitions between different extended Takahashi subspaces whose description is beyond the standard operational 
formalism of second quantization. 

Hence one may introduce creation and annihilation operators for the av bond particles provided that they are used 
in the expressions of the generators of the allowed processes and transitions. In addition, such operators must obey 
an algebra compatible with the av bond particles being hard-core bosons. Thus the operators av and Qj^av that 
create and annihilate one av bond particle at a av effective lattice site of index j — 1, N aav must anticommute on 
the same site, 

{f],avi flW> = 1 ; {gl au , g). au } = {gj,au, gj, av ] = , (137) 

and commute on different sites, 

[4,c«" & v >°"] = [4,o*'4',cJ = =0; ■ (138) 
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Such -qv (and sv) bond-particle operators can be expressed in terms of products of 77-spinon (and spinon) operators, 
Eq. (|67|) . Since their expressions involve only such operators, it follows from Eq. (|77|) that they commute with the c 
fcrmion operators, 

ifj,c > 9j>, av ] = [fj,c , 9j>, av ] = [fj, c > 9j',au] = [fj,c , 9r, av ] = , a = T),s, v = l, 00 . (139) 

Each site of the rjv (and sv) effective lattice refers to 2v sites of the ?7-spin (and spin) effective lattice. For each 
occupancy configurations of the latter effective lattice, the 2v sites of a rjv (and sv) bond particle are neither occupied 
by another rjv' (and sv') bond particle of a different v' ^ v branch nor by unbound 77-spinons (and unbound spinons). 
Moreover, the sites of the ?7-spin and spin effective lattice correspond to different sites of the original lattice. Hence 
it follows from the ?7-spinon and spinon operator algebra of Eq. (|79j) that the operators of av and a'v' bond particles 
corresponding to different av ^ a'v' branches commute between themselves, 

[Sj,c«" ft/',*'*/'] = Wj, av l9f >a > vl ] = [Si.m^Sj'.aV'] = J av ^ Oi V . (140) 

When acting onto their extended Takahashi subspace, the av bond-particle creation and annihilation operators 
appear in operators whose expressions contain an equal number of such creation and annihilation operators. However, 
they also appear in expressions of operators containing only av bond particle creation or annihilation operators. 
Specifically, this is the case of the operators that generate transitions between extended Takahashi subspaces belonging 
to the same 77-spinon and spinon subspace. Those can be written as a product of two operators. One of them adds 
or removes lattice sites to or from, respectively, the initial subspace av effective lattices. The other operator is one 
or a product of av bond-particle creation or annihilation operators suitable to the final subspace. We note that 
adding or removing av effective lattice sites also adds and removes, respectively, discrete momentum values from the 
corresponding av momentum band. 

Analysis of the spectral-weight distributions obtained by use of the PDT formalism [8l4li1| reveals that most r\v 
(and sv) bond particle spectral weight corresponds to a 2^-site compact domain of the former lattice. The spatial 
coordinate Xj — a nl/ j (and Xj = a sv j) of the r\v (and sv) bond particle refers to the position of such a domain middle 
point. However, more diluted and thus spread 2^-site 77-spin (and spin) effective lattice occupancy configurations of 
increasing length but centered at that point contribute as well to the r]v (and sv) bond particle. 

Fortunately, within the N a 3> 1 limit associated with our operator description the problem of the internal degrees 
of freedom of the composite av bond particles separates from that of their position in the corresponding av effective 
lattice. The latter corresponds to the spatial coordinate of the above ry-spin lattice (a — rj) or spin lattice (a = s) 
middle point of each av bond particle 2^-sites configuration. Hence here we do not need to address the very complex 
problem of the explicit expression of the r\v (and sv) bond particle operators in terms of products of Ty-spinon (and 
spinon) operators, which is implicitly accounted for by the BA solution. 

Within the Jordan- Wigner transformation method of Ref. (53[, the hard-core av bond particles are mapped into 
av fermions with operators given by, 

fj,au = e 9j,au i fj,av = & Qj av , (141) 

where, 

4>j,av = j "ji t0iv fj' ,otv <$>j' ,j,av ■ (142) 
J'¥=J 

For ID the phase (f>j>,j,au appearing here has for all av branches only two possible values, (j)j'^ tav — and (f>j>,j,au = n - 
Indeed, the relative angle between two sites of the av effective lattice in a ID chain can only be one of the two values. 
Then the av phase factor of Eq. (| 142|) is such that [53| , 

e i(0 i+ i, Q „-0j, Q „) _ e <ir/j )OV /i,af _ (143) 

The c fermion operators have the anticommuting relations given in Eq. (|76[) . On the other hand, the av fermion 
operators that emerge from the Jordan- Wigner transformation associated with Eqs. (|141l) and (|142p have similar 
anticommuting relations given by, 

{fj,av ' fj',av} — ; {fj,av ' fj',ai/} = {fj,<*v > fj',av} = . (144) 

In addition, it follows from Eqs. (|139p ~ (|142l) that the av fermion operators commute with those of the c fermions, 

lfj,c i /j'.oJ = lfj,c . fj>,ad = [fj,c » fj'^u] = [fj,c , fj', av ] = 0, a = r),s, v = 1, oo , (145) 
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and av and aV fermion operators corresponding to different av ^ olv' branches commute between themselves, 

[/U> fjw] = \ftav> f]>,<*u>] = f ?,*>*>] = ; «V . (146) 

One can introduce av fermion operators /J. a „ labeled by discrete momentum values qj such that j = 1, ... : N aau . 
Those are the conjugate variables of the av effective lattice real-space coordinates of site index j = 1, ...,N aav . It is 
confirmed below that their values exactly coincide with those of the BA solution discrete momentum values given in 
Eq. For subspaces for which the ratio N aau /N a involving the number N aav of sites of the av effective lattice 

given in Eqs. (jlOll) and (|102[) is finite, such operators are given by, 

N a 



ft = 1 V e + il }j a ^]' ft - = 1 jy 

1 

- -7vp= E e^ "" j j = l,...,JVa w - (147) 

In ID the phase factor e^ 3 > a " does not have any effect when operating before fj au . It follows that in ID the 
expression of the Hamiltonian does not involve the phases c/>j a v given in Eq. (|142j) . On the other hand, the Jordan- 
Wigner transformation phases 4>3,uii have direct effects on the boundary conditions. Those determine the discrete 
momentum values qj of both the c and av fermion operators of Eqs. (|63|) and (11471) . respectively. In ID the periodic 
boundary conditions of the original electron problem are ensured provided that one accounts for the effects of the 
Jordan- Wigner transformation on the boundary conditions of the c fermions and av fermions upon moving one of 
such objects around the chain of length L once. That in ID only the values 4>y .j,av — an< i <Pj',j,av = n are allowed 
in the cj>j iav expression given in Eq. (|142l) implies that in units of 2n/L the av fermion discrete momentum values 
either are integer numbers or half-odd integer numbers. 

We recall that the rotated-electron occupancies of the sites of the original lattice separate into two degrees of 
freedom only. Those of the = 2S C sites of the original lattice singly occupied by rotated electrons separate into 
(i) N c = Nfi = 2S C sites of the c effective lattice occupied by c fermions and (ii) iV a , = M s = iVjj, = 2S C sites of the 
spin effective lattice occupied by spinons. Those of the Np, = [N a — 2S C ] sites of the original lattice doubly occupied 
and unoccupied by rotated electrons separate into (i) = 25' 1 = N'^ = [N a — 2S C ] sites of the c effective lattice 
unoccupied by c fermions and (ii) N a = M v = 2S^ = = [N a — 2S C ) sites of the 77-spin effective lattice occupied by 
77-spinons. Indeed, the number values equalities N c = M s = and N' c l = M v = given in Eq. (|100[) do not imply 
a similar equality for the objects under consideration. Specifically, the N c c fermions (and c fermion holes) and M s 
spin-1/2 spinons (and M v ry-spin-1/2 ry-spinons) describe different degrees of freedom of the same rotated-electron 
singly occupied sites (and rotated-electron doubly occupied and unoccupied sites) whose rotated-electron site 
occupancy numbers are defined in Eqs. (pTO]) and ([TTj). 

The c fermions do not emerge from a Jordan- Wigner transformation. On the other hand, each av fermion corre- 
sponds to a set of 2v sites of the original lattice different from and independent of those of any other a'v' fermion. It 
then follows that upon moving around its av effective lattice of length L, a av fermion only feels the Jordan- Wigner- 
transformation phases of its own lattice. Those are associated with both the av fermions and av fermion holes. Hence 
its discrete momentum values obey the following periodic or anti-periodic boundary conditions, 



N a 

iqj L 



J^J J jg»(& + l,af-^,auO t e K4>i+-L,<*v-<l>i,C.u) I _ e i7T[JVa Q „-l] _ _ e i^N aav _ (148) 



Here the phase factor reads 1 and —1 for [N aav — 1] even and odd, respectively. The term —1 in [N aatl — 1] can be 
understood as referring to the site occupied by the av fermion moving around its effective lattice, which must be 
excluded. 

For the av fermions the unoccupied sites of their av effective lattice exist in their own right. Indeed, note that 
according to Eq. (|141j) the expression of both the creation and annihilation operators of such objects involve the 
Jordan- Wigner-transformation phase <j)j_ av . As a result, such a phase affects both the av fermions and av fermion 
holes. That justifies why the phase factor e l ^ Naa " of Eq. (|148[) involves all the N aav = [N aw + N^ u ] sites of the av 
effective lattice. The only exception is that occupied by the moving av fermion. Hence it involves both the [N av — 1] 
sites occupied by the remaining fermions of the same av branch and the corresponding N^ v av fermion holes. 

In contrast, the c fermions are only affected by the sites occupied by av fermions. Indeed, only the sets of 2v sites 
of the spin (and 77-spin) effective lattice associated with each occupied site of the sv (and r\v) v — 1, ...,00 effective 
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FIG. 1: The c fermion energy dispersion e c (q) = e°{q) - £° c {q Fc ) defined by Eqs. (fToTj) and ([162)) for /3 = c plotted as a function 
of the momentum in units of t for a set of U values (in units of t), electronic density n — 0.59, and (a) spin density m = and 
(b) limit of maximum spin density m — > n — 0.59. In the spin-density m — > n limit the shape of the energy dispersion e c (q) 
becomes independent of U/t. In both figures the vanishing-energy ground-state energy level is marked by a horizontal line, 
which overlaps the energy dispersion at the c Fermi points q = ±qFc ~ ±2/cf = ±7rn. 



lattices and the sites of the spin (and 77-spin) effective lattice occupied by unbound spinons (and unbound 77-spinons) 
correspond to sites of the original lattice whose degrees of freedom associated with the c hidden U(l) symmetry are 
described by the occupancy configurations of the c effective lattice. 

The unbound spinons and unbound 77-spinons of 77-spin projection +1/2 do not contribute to the states momentum. 
The momentum 7r associated with the unbound 77-spinons of 77-spin projection —1/2 remains unchanged for all Hilbert- 
space energy eigenstatcs. Consistent, unbound spinons and unbound 77-spinons do not undergo any Jordan- Wigner 
transformation. Thus they do not play any role in the boundary condition that determines the c band momentum 
values. On the other hand, due to the Jordan- Wigner-transformation phase <j)j>,a.v of each of the N av av fermions 
at sites j' = 1, ...,N aav of their av effective lattice, the c fermion discrete momentum values are determined by the 
following periodic or anti-periodic boundary condition, 

e iq 3 L = YIY[ e i ^'+^""-^'<^ = e n ^ N ^ . (149) 

QLV j'=l 

The phase factor on the right-hand side of Eq. (1149)) reads 1 and —1 for Ylau ^ av even an d odd, respectively. 

The above results imply that the discrete momentum values qj of both c and av fermions have the usual momentum 
spacing qj+i — qj = 2ir/L. Moreover, they are exactly those of the BA solution given in Eq. (IM1) . Here the same 
values of the quantum numbers I*" and J? given in Eq. (|113[) were imposed by the boundary conditions of Eqs. 
(11481) and (|149[) . This is fully consistent with the discrete momentum values of Eq. being within our operator 
formulation carried by the c fermions and av fermions as defined here. 

E. The general energy spectrum, the /3 fermion energy dispersions, and the Hamiltonian expression 

Let us denote the c and av fermions by /3 fermions, where ft = c, av. The corresponding momentum distributions 
Np(qj) appearing in the functional representations of the total momentum provided in Eq. (| 104|) and thermodynamic 
BA equations given in Eqs. (|107[) - ([109l) are the eigenvalues of the number operators, 

M<li) = fl,pf«.(i, P = c,av. (150) 

As mentioned in Section IV A[ such operators commute with the ID Hubbard model Hamiltonian. Consistent, the 
Bethe states are eigenstates of the number operators provided in Eq. (|150[) . Those have eigenvalues Np(qj) = 1 for 
occupied momentum values and Np(qj) = for unoccupied momentum values. The non-LWSs generated from such 
Bethe states differ from them only in the unbound 77-spinon and unbound spinon occupancies. Hence the c fermion and 
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av fermion operator algebras and corresponding occupancy configurations apply to all energy eigenstates, including 
those outside the BA solution subspace. 

The quantum numbers in the thermodynamic BA equations, Eqs. (|107|) and (1108|) . are the discrete momentum 
values qj of Eq. (|99|) . Those are eigenvalues of corresponding j3 translation generators qp, which commute with 
the Hamiltonian, momentum operator, and electron - rotatcd-clcctron unitary operator. That the rotated-electron 
occupancy configurations that generate the Bethe states, |^{ O) j0 )U ), are for u > independent of u is consistent with 
the corresponding f3 fermion momentum values qj occupancy configurations being also independent of u. On the 
other hand, the corresponding rapidity momentum functions k c (qj) and rapidity functions A ail (qj) appearing in the 
thermodynamic BA equations, Eqs. (|107|) - (|109|) . are the eigenvalues of corresponding operators, 

k c {qj); k av (qj), a = 7], s, v = \,...,oo. (151) 

Such operators also commute with the ID Hubbard model Hamiltonian. Thus the Bethe states are eigenstates of the 
rapidity operators. However, except for u — > oo they do not commute with the electron - rotated-electron unitary 
operator. 

The rapidity operators of Eq. (|151[) have a well-defined expression in terms of the sets of number operators N c (qj) 
and {N av (qj)} of Eq. (|150|) where a — r),s and v = 1, oo. The dependence of such rapidity operators on the latter 
number operators has exactly the same form as that of the corresponding BA rapidities k c (qj) and h av (qj) on the 
set of momentum distributions N c (qj) and {N au (qj)} where again a = to s and v = 1, oo. The latter dependence 
is defined by the thermodynamic BA, Eqs. (|107[) - (ll09p . Hence such equations define as well the rapidity operators of 
Eq. (Tim 

At both finite chemical potential and magnetic field the relation provided in Eq. (|B3|) of Appendix [B] between the 
energy eigenvalues E = {^i d^,u\H\^i j & , u ) and E symm = {^i ,i^, u \H symm \^i ,i^, u ) of the Hamiltonians flTJ and 
pj). respectively, can be rewritten as, 

E = E symm + M [M™ l/2 - M ™ 1/2 ] + H [M™ 1/2 - 1/2 ] . (152) 

As found in Ref. [l[ , here the energy eigenvalues E symm depend on u through the dependence of that parameter on 
the rapidity functionals defined by the thermodynamic BA equations, Eqs. f)10T|) - (|110[) . Within the present operartor 
formulation, the energy eigenvalues obtained in that reference can be expressed as functionals of the occupancy 
configurations of the /? = c, rjv, sv bands discrete momentum values qj as follows, 



E symm = -2t^2 WciUj) (cosk c (qj) +u)- u/2] 



i=i 



r f ^ 



(153) 



From the interplay of symmetry with the BA solution, one can extend the results of Ref. [l| for the energy spectrum 
of the ID Hubbard model in the BA solution subspace to its full Hilbert space. The energy of the unbound 77-spinon 
and unbound spinon occupancy configurations that generate the non-LWSs corresponds to the energy terms on the 
right-hand side of Eq. (1152[) other than E symm . On the other hand, the E symm expression provided in Eq. (1153j) is 
valid for all A Na energy eigenstates that span the ID Hubbard model Hilbert space. 

The limiting momentum values of Eq. (|117|) have simple expressions for the ground states and their excited states 
generated from them by occupancy-configuration deviations involving a finite number of our formulation objects. 
For simplicity and without loss in generality, we consider ground states with densities in the ranges n € [0, 1] and 
to G [0, n] for which, 

q c = Tr; q sl = k F -f ; q sv = [fof - k F {\ =irm, v > 1 ; q w = [it - 2k F ] =it(l-n), (154) 

where we have ignored corrections of order 1/N a . 

Consistent with the results of Appendix [B] concerning /3 fermion occupancies, ground states are described by 
compact c and si fermion finite occupancies for momentum values in the ranges q G [— q F p, +QF/3] where j3 = c, si. 
On the other hand, the av ^ si fermion branches have vanishing ground-state occupancy. For the above ground 
states with densities in the ranges n G [0, 1] and to G [0, n] the c and si Fermi momenta are given by, 

q Fc = 2k F ; q Fs i = k Fi . (155) 



Here we have again ignored corrections of order l/N a . 




FIG. 2: The si fermion energy dispersions e s i(q) = £°i(q) and e° 1 (q) defined by Eqs. (|161[) an d (|162[) for /3 — si plotted as a 
function of the momentum in units of t for a set of U values (in units of t), electronic density n = 0.59 and (a) spin density 
m — and (b) limit of maximum spin density m — > n = 0.59, respectively. The energy dispersions vanishing-energy ground- 
state energy level is marked by a horizontal line, which overlaps them at the si Fermi points q — ±gFsi = ±/cf = ±[7r/2] n at 
m — and q = ±gFsi = ±2&f = inn as rn — s- n, respectively. As discussed in the text, for general values of the spin density 
m 6 [0, n] the energy dispersion £si(q) < can for the absolute- value momentum range \q\ £ [qFsi,q s i] be identified with the 
two-spinon binding energy within the spin- neutral si fermion. The energy dispersion e®i(q) < plotted here in figure (b) for 
m — > n — 0.59 is thus such a binding energy. 



We denote by k®(qj) and ArJqj) the rapidity-function solutions of Eqs. (|107[) - (|109p that refer to a ground state. 
Those play an important role both in the PDT of Refs. (8l-tTo| and pseudofermion scattering theory introduced in the 
second paper. Upon suitable manipulations of Eqs. (|107l) - (|109[) . such ground-state functions k®(qj) and Ajg(<jj) may 
be defined in terms of their inverse functions qj = qj(Ag •) as, 

qj = ^(A°fe)) + (-l)^ f* dkic,J^,$^) , j = l,...,N ag , 

J-Q \ U U J 

F C °(A) = arcsin(A) ; F c °(A°fe)) = k° c ( qj ) , 

F°„(A) = <5 Q ,, ; 2Re [arcsin(A + ivu)] . (156) 

Here N afl = N a for f3 = c, N aav is given in Eqs. (|101[) and (|102l) . the functions $ Ci/ g(r, r') are defined in the second 
paper, and the parameter Q — tc®(2kp) is the value of the ground-state rapidity function fc^fe') a ^ the c Fermi point 
q Fc = 2k F - 

The P fermion momentum-distribution function deviation, 

SNpfa) = Npfa) - N^fa) , j = l,...,N afi , p = c,av, (157) 

also plays an important role both in the PDT of Refs. [8l-fl0| and pseudofermion scattering theory introduced in the 
second paper. In it, Np(qj) is the excited-state /3 fermion momentum-distribution function. Here Np(qj) is the initial 
ground-state j3 momentum distribution function. 

Within our general notation, the momentum deviations 5Np(qj) given in Eq. (|157[) are labelled by the /3 band 
momentum value qj. The corresponding /3 band momentum function distribution Np(qj) can only have the values 
Np(qj) — 1 for occupied momentum values and Np(qj) = for unoccupied momentum values. Therefore, the (3 
band momentum function distribution deviation SNp(qj) can only have the values SNp(qj) — 0, 8Np(qj) = — 1, and 
SNp(qj) = +1. Moreover, its values SNp(qj) = —1 and SNp(qj) = +1 can only occur provided that in the j3 band 
momentum function distribution Na(qj) the corresponding discrete momentum value qj, which may be shifted for 
the excited state, is occupied and unoccupied, respectively. Since for the discrete momentum values qj for which 
8Np(qj) = the /3 momentum distribution function deviation SNp(qj) does not contribute to the excitation spectra, 
we need only to define it by those for which SNp(qj) = —1 and 5Np(qj) = +1. (The set of (3 band momentum values 
qj left over are then those for which SNp(qj) = 0.) 

Within the present N a 3> 1 limit, for simplicity we often approximate the discrete momentum values qj by a 
continuum variable q. Consistent with [qj+i — qj] = 2ir /N a , in this case the above (3 band momentum distribution 
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function deviation values 5Np(qj) = —1 and 5Np(qj) = +1 become 6Np(q) = —2n/N a and 8Np(q) = +2n/N a , 
respectively. The BA quantum numbers Ij given in Eq. (|113|) where (3 = c, av, a = rj, s, and v — 1, oo are the (3 
band discrete momentum values qj — [2Tr/N a ] Ij in units of 2n/N a . Depending on Gp as defined in Eq. (|114[) being an 

even or odd integer number, the quantum numbers I? are integer or half-odd integers, respectively. Hence if under a 
transition from the ground state to an excited state the Gp value changes by an odd integer the corresponding (3 band 
discrete momentum values qj = [27r/iV a ] 1^ undergo a collective shift given by ±n/N a . Within the representation of 
the discrete momentum values qj by a continuum variable q, such an effect is captured by deviations SNp(q) — ±ir/N a 
for q at the /3 Fermi points. Those are additional possible j3 band momentum distribution function deviation values in 
addition to the above 5Np(q) — —2n/N a and 5Np(q) = +2ir/N a values associated with creation of one (3 pseudofermion 
hole and one /3 pseudofermion of momentum q, respectively. 

Consistent, the general form of the corresponding momentum distribution function deviation 5Np(q) is, 

6Np(q) = j^~^2S(q~ q p ) - —^2S(q- q h ); if 8Gp even for (3 = c, si , 

a p=l a h=l 

2tt N$ 2n N " 

= ^^8(q-q p )- —^2S(q-q h )±— ^ j8(q-jq FI3 ); if 8Gp odd for f3 = c, si , 

iVa p=l iVa h=l iVa j=±l 

2n N " 

8N p (q) = —J28(q-q p ); for p = av £ s\ , (158) 
a P =i 

where here and throughout this paper 8{x) denotes the usual Dirac delta- function distribution, qi,...,q N p are the 

excited-state momentum values of the Np added (3 fermions, q\, ...,q N h those of the added f3 fermion holes, and 

8Gp stands for the deviation in the value of the number given in Eq. (|114[) for the (3 = c, si branches. The additional 
term in the expression given in Eq. (j!58[) that emerges when the deviation SGp is for the /3 = c, si branches a odd 
integer number results from the corresponding /3 = c, si band momentum shift ±7r/iV Q . 

For ground-state - excited-state transitions for which the number of discrete momentum values N a/3 changes, it is 
considered that the removal or addition of such values happens after the f3 band momentum shift, if it occurs. The 
removed or added discrete momentum values are always located just before the excited-state limiting momentum 
— qp or just after the excited-state limiting momentum q@. Here qp is given in Eq. (| 1 1 T[) for (3 ^ c. For the c 
momentum band the limiting momentum values q^r may either be symmetrical, q+ +q~ = 0, or slightly asymmetrical, 



q+ + q~ = ±2ir/N a . (We recall that such values are provided in Eqs. (B.15)-(B.17) of Ref. (2l|.) The order of (a) 
the processes involving (3 fermion creation or annihilation and (b) the exotic processes involving (3 band discrete- 
momentum- values addition or removal is always such that the latter processes refer to unoccupied momentum values. 
Hence they are zero-momentum and zero-energy processes. This is why they are not explicitly included in the j3 
momentum distribution function deviations SNpfaj) as defined above. 

The use in the rapidity functional equations provided in Eqs. (I107p - (|108[) and general energy spectra given in 
Eqs. (|152[) and (jl53[) of distribution functions of form Np(qj) = Ng(qj) + 8Np{qj) and their combined and consistent 
solution up to first order in the deviations SNp(qj) given in Eq. (I157p . leads to the energy functionals, 

SE symm =EE4 fa) 5N f> (ij ) . (159) 

3=1 

and 



3 = 1 



N aav 

= E Y, £ ^ SN ^+ E E £ lMSN av (qj) 

+ 2| M | M„,_ 1/a + 2^ B \H\ [*Af a ,_i /2 - SN sl ] , (160) 
associated with such energy spectra, respectively. The energy dispersions £®(qj) appearing here can be expressed in 
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FIG. 3: The s2 fermion energy dispersion £^2(9) defined by Eq. (|161[) for /3 — s2 plotted as a function of the momentum 
in units of t for electronic density n — 0.59, limit of maximum spin density m — >■ n = 0.59, and a set of U values (in units 
of t). The energy dispersion u — > 00 limit corresponds to the vanishing-energy horizontal line. The s2 momentum band and 
corresponding energy dispersion do not exist at m = 0, both their momentum and energy bandwidths vanishing as m — s> 0. 
The energy dispersion £^2(1) < plotted here is the four-spinon binding energy within the spin- neutral s2 fermion. 



terms of the initial ground-state rapidity functions hPJqj) with A°(gj) = sinfc°(gj) for ft = c as follows, 



u 



2t 



dk<$> 



sin k AS (qj 



c,f3 



sin fc , 



E° c (qj) = ---2;cos/c c %), 
Cfe) = Sa, v {-uU + mB 



1 



(A°„(% 



(161) 



Furthermore, the energy-shifted dispersions £p(qj) are given by, 

£c(qj) = £°(g,) - e° c (q Fc ) = e°(g,) + (1 - 5 n>1 ) + <5„ 4 M ° - ^ 1^1 , 

e r ,u(qj) = el u (q 3 ) + (l-5 nA )2v\n\+5 nA 2v[i a 



e a sl (q 1 )-e a sl (qF S i)^e sl (q J ) + 2^ B \H\ 
e° sv {q j ) + 2v l i B \H\, 



(162) 



The functions $ Ci| g(r, r') appearing in Eq. (|161j) are defined in the second paper and the inverse functions of the 
singled value ground-state functions A° (qj ) are defined in Eq. (|156|) . As illustrated below in Section IV Fl the one- 
and two-electron energy excitation spectra can be expressed in terms of simple occupancies of the ft fermion energy 
bands defined in Eqs. ([1611) and (1T62]) . 

The standard operational formalism of second quantization can be applied to the av fermion for transitions within 
the extended Takahashi subspaces, which correspond to ai^-band particle-hole processes. Alike for the av bond-particle 
operators, the av fermion creation or annihilation operators also play a role in transitions between different extended 
Takahashi subspaces, whose description is beyond the traditional operational formalism of second quantization. We 
note that the standard formalism does not apply only in that the transition generators cannot be expressed in terms 
of av fermion creation or annihilation operators alone. However, they can be written as a product of two operators, 
one of which is a av fermion creation or annihilation operator. The other operator raises or lessens the number 
of discrete momentum values N aav of the initial subspace av momentum band. Concerning the excitation energy 
spectrum though, only the processes generated by the av fermion creation or annihilation operators play an active 
role. Indeed in the N a — > 00 limit raising or lessening the number of the av band discrete momentum values is a 
zero-energy and zero-momentum process. This is why the av fermion energy dispersion defined in Eqs. (| 16 1[) and 
(|162l) is the bare energy for addition/removal of one av fermion. 

For electronic densities n ^ 1 the values of the chemical potential /1 = fJ-(n) and magnetic-field energy 2\x B H = 
2/ib H(m) appearing on the right-hand side of Eq. (|160[) , whose ranges are given in Eqs. (|Cip and (|C6[) of Appendix 




FIG. 4: The rjl fermion energy dispersion eSife) denned by Eq. (|161[) for /3 = 77 1 plotted in units of t for a set of U values 
(in units of t), electronic density n = 0.59, (a) spin density m = and (b) limit of maximum spin density m — > n = 0.59. The 
energy dispersion u — > 00 limit corresponds to the vanishing-energy horizontal line. The energy dispersion £°i((?) > plotted 
here is the two-??-spinon anti-binding energy within the 77-spin-neutral rjl fermion. 



[Cl respectively, are fully controlled by the energy dispersions of the (3 — c,sl fermion branches defined in Eq. (11611) . 
They read, 



(j, = sgn{(l - n)} 



2/i B H = sgn{m}e° sl (q Fsl ) , n^l. (163) 



The expressions given here are valid for the whole range of densities n ^ 1 and to. 

The energy functional of Eq. (|160[) includes interaction terms, which are of second order in the deviations SNp(qj), 
beyond those given in that equation. Such terms correspond to 1 /N a corrections to the energy functional provided in 
Eq. (|160p . Due to the occurrence of an infinite set of conservation laws associated with the model integrability in the 
BA solution subspace [41H431 |54| . such interactions are of zero-momentum forward-scattering type. Such second-order 
in the deviations SNp(qj) interaction energy terms and related dressed phase shifts, dressed S matrix, and scattering 
theory are issues studied in the second paper. 

The energy dispersions e c (q), e s i(q) and £®i(q), ^(q), £ni(q), and s^ 2 {q) ar e plotted as a function of the momentum 
q in Figs. (TJ [^1 El SI and[5j respectively, for several U/t values, and electronic density n = 0.59. The two dispersions 
plotted in each of the Figs. [THUSl and [5] refer to the spin-density limits to — > and m — > n = 0.59, respectively. In 
the to — y limit the momentum and energy bandwidths of the s2 fermion energy dispersion vanish, so that in Fig. [3] 
it is plotted only for the spin-density limit to — > n = 0.59. For to = the electronic density n = 0.59 is that used in 
the studies of Ref. [ll[ of the photoemission spectrum of the quasi-lD metal TTF- TCNQ. The figures U/t values are 
U/t = 0,2,4.9,5.61, 10,oo. Here U/t = 4.9 and U/t = 5.6 are the values suitable to the spectral lines of the TCNQ 
and TTF molecular stacks, respectively, in the investigations of Ref. [ill ]. 

Limiting expressions of the energy dispersions s c (q), e s i(q), and £^ v (q) corresponding to u — > and u 3> 1 
are given in Eqs. (|C12I) - (|C14[) of Appendix [C] for electronic density n 6 [0, 1] and spin density m = 0. In the 
maximum spin-density limit, to — > n, all /3 fermion energy dispersions have simple closed-form expressions. Those 
are given Eqs. (|C15[) - flC16|) of Appendix [Cl In that limit the ground-state si band momentum limiting values read 
±q s i = ±kp^ = ±2kp = ±7r n and have the same values as those of the remaining ground-state sv band momentum 
limiting values, ±g s „ = ±[fcj"|- — = ±2kp = ±irn. On the other hand, the ground-state rjv band momentum 
limiting values are in that limit given by their general n dependent expression ±q VL , — ±[tt — 2kp\ — ±tt [1 — n]. 

The corresponding av fermion dispersion energy bandwidths are given by, 

W av = \el v {q av ) - e° av {Q)\ = | £ce!/ (0)| ; W c = At . (164) 

The c fermion energy dispersion bandwidth, W c = 4i, is independent of n, to, and u. On the other hand, from 
manipulations of Eqs. (|16ip and (|T62|) we find that for av ^ si branches the av fermion energy dispersion bandwidths 
can be expressed in terms of the unbound TOspinon (a = 77) and unbound spinon {a — s) energy scale [e a +1/2+^0,-1/2] 
given in Eqs. (|9Tj)) and (|9T1) as follows, 

Wav = i w [e a ,+i/2 + £0,-1/2] ! <£,,„<!, av^sl, a = 77, s . (165) 




FIG. 5: The rfl fermion energy dispersion £® 2 (<lj) defined by Eq. (|161[) for ft = rfl piotted in units of t for a set of U values (in 
units of t), electronic density n = 0.59, and (a) spin density m — and (b) limit of maximum spin density m n = 0.59. The 
energy dispersion u — ¥ oo limit corresponds to the vanishing-energy horizontal line. The energy dispersion £^2(3) > plotted 
here is the four-jy-spinon anti-binding energy within the jy-spin-neutral r)2 fermion. 



For a ground state with densities n £ [0, 1] the number i av appearing here decreases continuously for increasing values 
of u. For a = 77 it has the limiting behaviors i niJ — > 1 for u — > and — > for u — > 00. On the other hand, for a = s 
and v > 1 it has the limiting values i sv — > 1 for u — > and i sv — > 1/v for it — > 00. 

For a ground state with densities n £ [0,1] and to € [0, n] the corresponding si band is unoccupied for \q\ € [fcp-j,, fc^t] 
and occupied for |q| e [0,fci?j_]. Again from the use of Eqs. (|161|) and (|162|) . we find that the si energy bandwidth 
W s i defined in Eq. (|164l) can for av = si be expressed as, 



w sl = wA + wji 

= [e.,+i/2 + e 



-,-1/2] = 



= [e al (*n)-e-i(*n)]; w !i 

2fi B \H\ . 



[£si(kFi) - £ s i(0)] 



(166) 



Here M 7 ^ and Wf 1 are the energy bandwidths of the ground-state si fermion hole and si fermion seas, respectively. 
For the electronic density range n £ [0, 1], their limiting behaviors are Wgi — > for m — > 0, Wj\ — > W s i for m — > n, 
Wfj VF s i for rn 0, and -> for to -> n. 

For n G [0, 1] and m — > n the energy bandwidths VF sl , and W vl , have closed-form expressions, provided in Eqs. 
(IC17I) - (|C20|) of Appendix [C] At n — 1 and to = the W sl = bandwidth u dependence is for the whole u > 
range given in Eq. (|C21[) of that Appendix. The behaviors for both u — > and u — > 00 of the energy bandwidths W su 
and Wrjv are reported in Eqs. (|C22j) and (|C23|) of Appendix [C] for electronic densities n £ [0, 1] and both spin densities 
to = and m — > n. That such energy bandwidths vanish as u — > 00 follows from all spin and ry-spin configurations 
being degenerate in that limit. In the limits n — > 1 and m — > the rjv bands and v > 1 sv bands collapse into the 
point £^(0) = and £^(0) = 0, so that W vv — > and W sv — > 0, respectively. 

The BA solution performs the sum of all Hamiltonian terms on the right-hand side of both Eqs. (|T2")) and ([75)1 and 
beyond. The corresponding expression of the ID Hubbard model Hamiltonian H syrnm appearing in such equations 
and in Eqs. ([TJ and ([2j is simply obtained by replacing the /3 fermion numbers Np(qj), rapidity momentum functions 
k c (qj), and rapidity functions A au (qj) in the general energy expression of Eq. (|153|) by the corresponding ft fermion 
number operators, c rapidity momentum operators, and av rapidity operators, respectively. This gives, 



H 



symm 



N a 

-2t^2 [fic(qj) (cosfccfe) +uj -u/2 



3=1 



u=l j=l 



Re 



1 - (A nu (qj) + ivu) 2 \ - v 



(167) 



As discussed above, the c rapidity momentum operators and av rapidity operators are implicitly defined by the 
thermodynamic BA equations, Eqs. (I107p - (|109[) . 
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F. Simple excitations spectra in terms of the c and si fermion energy dispersions 

The one- and two-electron excited states mostly considered in the vast literature of the ID Hubbard model are 
those of an initial ground state with electronic density in the range n£ [0,1] and vanishing spin density, m — 0. The 
energy spectrum of such excited states has simple expressions in terms of the c and si energy dispersions defined in 
Eq. (|162l) . Here we revisit some of such excited states. One of our aims is to describe the microscopic processes that 
generate such states in terms of the elementary objects of our operational formulation. The corresponding spectra 
are expressed in terms of such objects energy dispersions. 

Consistent with the results of Appendix [B] that m = ground-state occupancies of the objects of our formulation 
are, 

M s = M b s ° = N = nN a ; M s bo = N ; M™ ±lJ2 = , 
M v = M% n = M™ 1/2 =N a -N = (l-n)N a , 

N c = N; N sl = N/2 ; N> s \ = . (168) 

For simplicity, we consider that N/2 is an odd integer number, so that N c — N and N s i = N/2 are even and 
odd integer numbers, respectively. Such a m = ground state has no unbound spinons and no si fermion holes, 
M™ = and Nl\ = 0, respectively. For it all si band momenta q E [—c[Fsi> +QFsi] = [— &f, +&f] where kp = nn/2 
are occupied. Moreover, the c band momenta are occupied for q € [— qFc, QFc] = [—2kp,2kp] and unoccupied for 
\q\ £ [iFctQc] = [2fc_F,7r] where 2kp — ixn. Here we have ignored 1/N a order corrections, which vanish within the 
thermodynamic limit. The fx = 0, n = 1, and m = absolute ground state is a particular case of the m = ground 
state for which the number values of Eq. (|168l) apply with n — 1 and thus N — N a . For the absolute ground state 
the c band is full and there are no unbound 77-spinons. 

The N%i = 2 spin-singlet excited states of such a m = and N^ = ground state have deviations given by, 

SN C = -6M™ 1/2 = Q; 6N sl = -2; SN S \=N^=2; 6N s2 =N s2 + l, 

±ir/L — c — band — momentum — shift . (169) 

The number values of such excited energy eigenstates are the sum of those of the corresponding initial m = 
ground state provided in Eq. (|168[) plus their deviations given here. Nearly all the spectral weight associated with 
the spin-singlet excitations is originated by the processes associated with such deviations. Within the ground-state 
transitions to such excited energy eigenstates, the c band momentum values undergo alternative collective shifts, 
±7t/L. Those lead to a macroscopic momentum contribution given by ±2kp — ±itn. Within such transitions, the 
number of annihilated si fermions exactly equals that of holes emerging in the si band, so that the number of discrete 
momentum values N agl is preserved. 

Under such = 2 spin-singlet excitations the spinon pairs of the two annihilated si fermions are broken. The 
four spinons left over reorganize in the final excited state within a spin-singlet four-spinon composite s2 fermion 
configuration. The emerging composite s2 fermion has vanishing momentum, energy, binding energy, and velocity. Its 
s2 effective lattice is reduced to a single occupied site, which corresponds to four sites of the underlying spin effective 
lattice. 

The spectrum of such spin-singlet excited states is degenerate with that of the spin-triplet excited states, whose 
deviations are given by, 

6N C = -SM™ +1/2 = ; SN sl =-l; $jv£ = = 2 ; 5M" n = M"" = 2 , 
±tt/L — c — band — momentum — shift , 

=F7r/L — si — band — momentum — shift . (170) 

Nearly all the spectral weight associated with the spin-triplet excitations is originated by the processes associated 
with such deviations, which generate exact energy eigenstates. The c and si band momentum alternative collective 
shifts ±tt/L and ^fir/L, respectively, result here from the changes in the boundary conditions that determine the 
integer of half-odd integer values of the BA c and si quantum numbers, Eq. (|113p . Those of the c band lead to a 
macroscopic momentum contribution given by ±2kp = ±7rn. On the other hand, the si band momenta collective 
shifts ±7r/L generate a macroscopic momentum ±kp = ±7rn/2. For the lowest-energy states of such an excitation 
branch that microscopic momentum exactly cancels another momentum contribution, ^fkp = =p7rn/2. The latter 
momentum arises from the emergence of a hole at si band momentum q = ±kp = ±nn/2, due to the si fermion 
removal. Therefore, the si band momentum collective shifts do not lead to any net momentum contribution in the 
general spin spectrum given below. 
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Spin spectrum energy-momentum domain Spin spectrum energy-momentum domain 

n=0.59,U=4.9 n=l, U=4.9 




FIG. 6: The (k, tj) plane domain of the spin spectrum of Eq. (|171[) with the energy in units of t for U/t = 4.9, spin density 
m = 0, and electronic densities (a) n = 0.59 and (b) n = 1. The shape of this spin spectrum is determined by the occupancy 
configurations of two holes in the si fermion energy dispersion plotted in Fig. [2] (a). 



A second si band hole appears at the opposite si band limiting side, due to the emergence of a discrete momentum 
value, which is associated with the deviation SN aBl = [SN s i + SNfjj] = 1. This process does not generate any 
momentum but drives the system into a virtual unphysical state whose si band limits are not symmetrical. For all 
4 Na energy eigenstates the av momentum bands are symmetrical around zero, with the limiting momentum values 
reading ±q aL , where q av is given in Eq. (|117[) . The above si band momentum collective shift ±7r/L restores that 
symmetry and leads to the final physical state. (Actually, all such processes cannot be separated, so that the virtual 
state is an artifact of our separation of the overall process into intermediate virtual processes.) 

In contrast to the iVj\ = 2 spin-singlet excited states, the number of discrete momentum values N asl of the spin- 
triplet states is not preserved. However, since for the spin-singlet excited states there is no si band momentum 
collective shift, in the thermodynamic limit their energy and momentum spectra are exactly the same. Both the 
transitions from the initial m = ground state to the spin-singlet and spin-triplet excited energy eigenstates whose 
number deviations are given in Eqs. (|T69|) and ([T70]) involve the emergence of two holes in the si fermion band. 
(The spin weight originated by four-sl-band-hole excited energy eigenstates is very small for u finite and vanishes as 
u — > co.) 

The degenerate spectrum of the = 2 spin-singlet and spin-triplet excited states can be expressed as, 
w(fc) = -e sl (q) - e sl {q') ; k = [±2k F - q - q'] = [±nn - q - q'] , 

q,q' G [-k F ,k F ] = [-• 7rn/2,7rn/2] . (171) 

Its {k, uj) plane domain is represented in Fig. [B]for U/t — 4.9, spin density m = 0, and electronic densities (a) n = 0.59 
and (b) n = 1. It is generated by all two-hole si band occupancy configurations of the si fermion energy dispersion 
plotted in Fig. [1(a). 

The number and deviation values given in Eq. (|170l) apply to the three branches of spin-triplet excited states. 
For that inside the BA solution subspace the excited-state unbound-spinon number 5M^ n — M"" = 2 refers to 
-^"+1/2 = 2; M"™ 1 , 2 — 0. The remaining two spin-triplet excited state branches are outside the BA solution 
subspace, yet are accounted for by the present operator formulation. Specifically, they have unbound-spinon numbers 
-^T+i/2 = l'-^C-i/2 = 1 an d -^T+i/2 = O'^-C-i/2 = respectively, and are generated from the corresponding 
= 2; ^™-i/2 = spin-triplet excited state inside that subspace by spin-flipping one and two unbound 
spinons, respectively. 

Finally, we consider the simplest one-electron excited states of the \x — 0, n = 1, and m = absolute ground 
state. Within our rotatcd-clcctron operational description, the one-electron lower-Hubbard-band (LHB) and upper- 
Hubbard-band (UHB) are well-defined concepts. Specifically, the one-electron LHB is generated by transitions from a 
n < 1 (and n > 1) ground state to excited energy eigenstates with no rotated-electron doubly occupied sites (and no 
rotated-electron unoccupied sites). For an initial n — 1 ground state such transitions do not exist. The first UHB is 
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generated by transitions from n < 1 (and n > 1) ground states to excited energy eigenstates with one rotated-electron 
doubly occupied site (and one rotated-electron unoccupied site). For one-electron addition excitations, the spectral 
weight generated by transitions to excited energy eigenstates with two rotated-electron doubly occupied sites is very 
small. Those transitions generate the second UHB. 

Here we discuss only the case of the one-electron addition first UHB of the initial fj, = Mott-Hubbard-insulating 
absolute ground state. By simplicity, we call UHB the first UHB. A selection rule prevents the occurrence of one- 
electron addition excited states involving creation of r\v fermions onto that ground state. Indeed, the corresponding 
SN = 1 deviation is in that case given by 5N — [6N C + Y^vLi 2v8N rju ] = 1. On the other hand, the corresponding 
deviation in the 77-spin value reads SS V — ~^[6N C + Y^L^ 2vSN v ^] = — |. Since S v > 0, the latter deviation is only 
allowed for initial ground states whose 77-spin value obeys the inequality S v > 1/2. The one-electron addition UHB of 
the the fi = absolute ground state is generated under creation onto it of one unbound 77-spinon of 77-spin projection 
— 1/2, which is an allowed process. 

The number values of the Mott-Hubbard insulator m = 0, n = 1, and fi — absolute ground state are those 
provided in Eq. (|168j) for N = N a . The UHB one-electron addition spectral weight associated with the creation of 
one rotated-electron doubly occupied site onto that state is generated by processes whose number deviations are given 
by, 

SN C = -1 ; 5N sl = -1 ; 5N& = 1 ; 8M%_ 1/3 = ; 5M™ 1/2 = 1 , 

±7r/L — c — band — momentum — shift . (172) 

The ground-state transitions to such excited states involve a collective momentum shift ±ir/L in the c band. Since 
the initial ground state is a half- filled state, the lattice translational symmetry is explicit. Thus the momentum 
shifts +7r/L and —n/L lead in this case to the same c band occupancy configuration. The macroscopic momentum 
±7T generated by such a momentum shift is exactly cancelled by the momentum ±n of the created —1/2 unbound 
ri-spinon. 

The corresponding energy gapped spectrum shape is determined by the emergence of one hole in the c momentum 
band, one hole in the si momentum band, and the creation of one unbound 77-spinon of 77-spin projection —1/2. It 
reads, 

u(k) = e v _ 1/2 - e c (q) - e sl (q') ; k=[-q-q']; <?e[-7r,7r]; q' € [-tt/2, tt/2] . (173) 

Since the initial ground state refers to vanishing chemical potential /1 = 0, the unbound 77-spinon energy is £,,,-1/2 = A' , 
as given in Eq. (|88l) . The creation of the —1/2 unbound 77-spinon and corresponding rotated-electron doubly occupied 
site is behind the spectrum being gapped. 

The N + 1 = N a + 1 ground states are among the excited states whose number deviations are given in Eq. (|1 T3[) . 
They are generated by processes under which the c fermion hole and si fermion hole are created at a c Fermi points 
q = =p7r and si band limiting momentum values q' = ±7r/2, respectively. In that case the general spectrum given in 
Eq. (I173P simplifies to w(k) = e 7)j _i/2 = /-*° and k = ±ir/2 = ±kp. The (k,u>) plane domain of the UHB one-electron 
addition spectrum provided in Eq. (|173p is represented in Fig. [7] The figure refers to a reduced-zone scheme for spin 
density m — 0, electronic density n=l, and (a) U/t — 4.9 and (b) U/t = 13.6. 

G. Confirmation of the anti-binding 77-spinon (a = 77) binding spinon (a = s) character of the composite olv 

fermions from analysis of their energy dispersion 

The unbound 77-spinons and unbound spinons whose energies are given in Eqs. (|87|) . (l88l) . and (|89f have no energy 
dispersion. They neither interact among themselves nor with the /3 fermions. In the second paper it is confirmed that 
unbound 77-spinons and unbound spinons are neither scatterers nor scattering centers. The energies of such unbound 
objects play an important role as reference energy scales for the energies of the 2^-77-spinon composite 77^ fermions 
and 2/y-spinon composite sv fermions. 

The unbound 77-spinon and unbound spinon energies of Eqs. (|8"T|) . (|55| . and (|59")l and the energy dispersions e v ^(q) 
and £ S v{q) of Eq (|162[) alike are measured relative to the zero-energy level of the Hamiltonian H of Eq. ([1]), which 
is that of the initial ground state. The zero-energy level of the energy dispersions £®„(q) and £°„((?) of Eq. (|161[) 

corresponds in turn to that of the Hamiltonian H symm of Eq. ([2]). Relative to it, the corresponding unbound spinon 
and unbound 77-spinon energies read, 



£ a ,±i/2 = °> a = r),s. 



(174) 
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FIG. 7: The (k,uj) plane domain of the UHB one-electron spectrum of Eq. (|173[l within a reduced-zone scheme with the energy 
in units of t for spin density m = 0, electronic density n = 1, and (a) U/t = 4.9 and (b) U/t = 13.6. The shape of this spectrum 
is determined by the combined one-hole occupancy configurations of the c and si fermion energy dispersions plotted in Figs. 
[T] and [2] (a), respectively, plus the finite energy associated with creation of one —1/2 unbound 77-spinon. The latter energy is 
behind the spectrum being gapped and refers to creation of one rotated-electron doubly occupied site. 



On the other hand, relative to the zero-energy level of the initial ground state the rjv and sv fermion energy 
dispersions are such that, 

mme w {q) = e nv (±q vi> ) = v[e 1h+1/2 + -1/2] = (1 - <W) 2 ^ M + <*n,i 2i V > 

ma,xe slJ (q) = e sv {±q sv ) = v[e s ^ +1/2 + £3,-1/2] = \H\ . (175) 

This confirms the anti-binding and binding character of the 2^-ry-spinon and 2i/-spinon occupancy configurations of 
a composite rjv fermion and composite sv fermion, respectively. The energy scales on the right-hand side of the two 
equations given here are exactly v = 1, 00 times larger than the energy of a unbound 77-spinon pair (and unbound 
spinon pair), which is provided in Eq. (|87|l (and Eq. (|91l0 . That the minimum (and maximum) magnitude of the 
energy dispersion e vv (q) (and s sv {q)) is reached at the limiting momentum values ±q V i, (and ±q sv ) is below found to 
be related to important r\v fermion (and sv fermion) symmetries. 

The inequalities given in Eq. (|175|) are consistent with the behavior of the energy dispersions plotted in Figs. [21 
[31 [31 andO For m = and H = 0, the two energy dispersions £ s i(q) and £°i(q) have the same zero-energy level, 
so that e s i(q) = £gi(q)- They are plotted in Fig. [2] (a). On the other hand, for m ^ such two energy dispersions 
are as given in Eq (|162[) related as £ s i(q) = £gi(<?) + 2/jb \H\. Hence the zero-energy level of the energy dispersion 
e si(<z) plotted in Fig. H](b) for m — > n = 0.59 and H = H c corresponds to the energy H c of the dispersion e s i(q). 
Zero and 2/is H c are the maximum two-spinon si fermion energies of such dispersions, respectively. Moreover, the 
zero-energy level of the s2 fermion energy dispersion plotted in Fig. [3jfor m — > n = 0.59 and H = H c corresponds to 
the energy H c of the dispersion e S 2 (q)- That is the maximum four-spinon s2 fermion energy of that dispersion. 
The zero-energy levels of the rjl and i]2 fermion energy dispersions plotted in Figs. 2] and [5] correspond to the energies 
2|/j| and 4|//| of the dispersions £ v i(q) and £^2(9), respectively. 

The binding energy of one sv fermion with v > 1 spinon pairs vanishes at the sv band limiting momentum values 
q = ±q sv - Its maximum magnitude is reached at q = and thus vanishing velocity v s „(q) — 0. At q = ±q sv the sv 
fermion velocity also vanishes, v sv (±q su ) — 0. Relative to its energy at such limiting momentum values there is at 
q = an energy gain of —W su . That energy gain equals the maximum binding energy. 

It is useful for the definition of the binding energy at arbitrary sv band momentum values q G (—q su , q su ) to consider 
a sv fermion dispersion e k sl/ {q) = [e2„(<z) + Its zero-energy level coincides with that of the momentum q = 0. 

For small q values that energy plays the role of sv fermion kinetic energy. At small q the sv fermion acquires a finite 
velocity, v su {q). The corresponding energy s k sv (q) is an increasing function of Its maximum magnitude is reached 
at ±q sv , exactly when the binding energy vanishes. Out of the two terms of the energy £ k u (q) — [^(q) + W s „], the 
energy dispersion £ all (q) plays the role of binding energy. As mentioned above, its absolute value is maximum at rest 
when q = 0, v av (q) — 0, and £g„(0) = —W sv . For small q values, part of the sv fermion binding energy is transferred 
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over to the kinetic energy e^(q) = [e® v (q) + W sv \, which becomes finite. The amount of transferred energy exactly 
equals the decrease in the energy scale (e^^)! absolute value. Such an energy transfer occurs as well for larger q 
values, although for them the energy e^(q) is not anymore a free-like kinetic energy. When the energy e* v {q) reaches 
its maximum magnitude, W sv , at q = ±q su all the binding energy has been used up. 

The binding energy of one si fermion created onto a ground state may also be identified with the energy scale 
e^^q). However, due to the finite si fermion ground-state occupancy, this refers now only to the si fermion-hole si 
band momentum ranges q e [— q s %, — qFsi] and q <E [QFsiiQai]- At the si band Fermi points q = ±9fs1 it reaches its 
maximum absolute value, W% v . For \q\ £ [qFsiiQsi] it is a decreasing function of \q\, vanishing at q = ±q s %. A similar 
analysis reveals that for the nv fermions the anti-binding energy can be identified with the energy scale (q). The 
arguments involve now energy losses rather than gains. That energy vanishes at q — ±q V i,. Its maximum value, W^ v , 
is reached at q = 0. 

Hence the energy dispersions £% v {q) and s° av {q) defined in Eq. (I16ip . whose zero-energy level refers to that of 

the Hamiltonian H S ymm of Eq. @, have a well defined physical meaning. Consistent with the above analysis, 
upon combining the unbound ?7-spinon and unbound spinon energy-scales expressions given in Eqs. (|90[) and (|91l) . 
respectively, with such energy dispersions expressions provided in Eq. (|162p . one finds that they can be written as, 

e °i/(<?) = £ »?r/(?i) ~ v K -1/2 + £77,-1-1/2] > , 

£L(<?) = - v[e a _ 1/2 +e s ,+i/ 2 ] < 0. (176) 

The maximum absolute value of the anti-binding energy (a = 77) or binding energy (a = s) of one av fermion 
created onto a ground state is the energy bandwidth of that state unoccupied av momentum-band sea, W^ v . This 
holds for creation of one si fermion. For the av ^ si branches this follows from the equality W av — W^ u . A related 
quantity of physical interest is the maximum absolute value of the r\v fermion (and sv fermion) anti-binding energy 
per 77-spinon pair (and spinon binding energy per spinon pair). It is given by w^ v = W^ v /v and can be expressed as 
follows, 

l° av £[l/v,l/v 2 }, 

-- l/v 2 , u — ^ 00 . (177) 

Here W al = 2v/j,b \H\ is the si band hole energy bandwidth provided in Eq. (I166P and W\/ x = W v \. The latter energy 
scale reads = 2\fi\ for u — >• and vanishes in the limit u — > 00, as given in Eq. (|C23|) of Appendix [C] 

The dependence on the number of pairs v = 1, 00 of the quantity of Eq. (|177p . which is the energy w^ v in 
units of W£i , reveals that the maximum absolute value of the nv fermion (and sv fermion) anti-binding energy per 
?7-spinon pair (and spinon binding energy per spinon pair) decreases for increasing pair number. Such an effect is 
stronger for larger values of u, with the v dependence smoothly changing from l/v for u -C 1 to l/v 2 for u 1. 



1° = l/v, 7i^0; C 



H. Transformation laws of the av fermions under the electron - rotated-electron unitary transformation 

Except for u — > 00 , the numbers of electron doubly occupied sites, electron unoccupied sites, spin-down electron 
singly occupied sites, and spin-up electron singly occupied sites are not good quantum numbers. As confirmed in 
Ref. 22], in general the expectation values of such numbers are u dependent. Only in the u — > 00 limit do they 
equal the corresponding rotated-electron w-independent eigenvalues. The studies of that reference used objects in 
one-to-one correspondence to those of our description. In that reference such objects were directly associated with the 
BA quantum numbers and the spin and 7/-spin tower of states. (The relation of their operators to those of the original 
electrons was lacking and no relation to the representations of the model global .SO (3) <Ei SO (3) ® U(l) symmetry was 
established.) 

Combination of the results of Ref. [22[ with those of the present paper provides important information about the 
transformation laws of the av fermions under the electron - rotated-electron unitary transformation. If one accounts 
for the different notations used here and in that reference, its results confirm that the deviations in the numbers of 
rotated-electron doubly occupied sites, rotated-electron unoccupied sites, spin-down rotated-electron singly occupied 
sites, and spin-up rotated-electron singly occupied sites due to creation or annihilation of c fermions, av fermions, 
unbound ?7-spinons, and unbound spinons are u independent. They exactly equal the expected rotated-electron 
numbers eigenvalue deviations. 

Furthermore, the results of Ref. [22| confirm that except for u — > 00 the corresponding deviations in the expectations 
values of the numbers of electron doubly occupied sites, electron unoccupied sites, spin-up electron singly occupied 
sites, and spin-down electron singly occupied sites are a function of u. They are different from the deviations in the 
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corresponding rotated-electron numbers eigenvalues. That is consistent with for finite u values such electron numbers 
not being good quantum numbers. 

There are two exceptions, though. As discussed in Section llVB| the first results from the ±1/2 unbound 77-spinons 
and ±1/2 unbound spinons being invariant under the electron - rotated-electron unitary transformation. For the whole 
u > range, creation onto a ground state of one —1/2 unbound ry-spinon, one ±1/2 unbound ?7-spinon, one —1/2 
unbound spinon, and one ±1/2 unbound spinon then exactly leads to a u-independent deviation in the expectation 
values of the corresponding electron occupancy numbers. Specifically, the deviations in the numbers of electron doubly 
occupied sites, electron unoccupied sites, spin-up electron singly occupied sites, and spin-down electron singly occupied 
sites, respectively, exactly equal one. This results from the unbound 77-spinons and unbound spinons corresponding to 
the same 77-spin and spin degrees of freedom occupancy configurations, respectively, in terms of electrons and rotated 
electrons. 

There is a second exception. It refers to the av fermions whose anti-binding energy (a = rj) or binding energy 
(a = s), and velocity, v au (q) — de^ Jtl/ (q)/dq = de av (q)/dq, vanish. The requirement of av fermion vanishing 

anti-binding or binding energy is fulfilled at limiting momentum values q — ±q av . At them the av fermion dispersions 
are given by, 

£nu(±q V v) = v K ,-1/2 ± ^,+1/2] = (1 - S n ,i) 2v\fi\ + S nt i 2vfi° ; e° v (±q TI „) = Q, 

E sv (±q su ) = v[e s _ 1/2 +e Si+1/2 }=2v^ B \H\; £°„(±q sv ) = . (178) 

The results of Ref. [22} confirm that the av fermions with limiting momentum values q — ±gw are invariant under 
the electron - rotated-electron unitary transformation provided that their velocity vanishes, 

. , _ de av {g) , del v {q) , 

^ - Q q \l=±9c* = Q q \q=±q a „ = ! 

±q av = ±y (N aav - 1) = ±— (1 - l/N aav ) for N Uav odd , 

±q av = ±1 JV 0o „ = ±— for N aaiJ even. (179) 

L &OLV 

In general we use units of lattice constant a, so that N a = L. On the other hand, here we have expressed the limiting 
momentum value of Eq. (|117p in terms of the av effective lattice spacing a av , Eq. (|1 36[) . This reveals that provided 
that N aav /N a < 1 is finite and thus N aau ^> 1 for iV a 3> 1, the limiting momentum values ±g Q „ play the role of the 
Brillouin zone limits ±7r/a Qi , of the av momentum band. 

The results of Ref. [22| confirm that for the whole u > range creation onto a n < 1 (and n > 1) ground state 
of one T]v fermion fulfilling the above two requirements leads to a deviation in the expectation value of the number 
of electron doubly occupied (and unoccupied) sites that exactly equals the number of 77-spinon pairs v. Similarly, 
creation onto a m > (and m < 0) ground state of one sv fermion meeting the two above requirements leads for 
u > to a deviation in the expectation value of the number of spin-down (and spin-up) electron singly occupied sites 
that again exactly equals the number of spinon pairs v. In contrast, creation of one av fermion that does not meet 
both requirements leads to it-dependent deviations in such electron-number expectation values. 

The av fermion vanishing velocity requirement is fulfilled at the limiting momentum values, q — ±q av , provided 
that the initial ground state has a vanishing av fermion density, N av /N a —¥ 0, as N a 00. This requirement is met 
for u > and arbitrary electronic density and spin density values by all av branches other than the si branch. For 
the si branch it is met for u > in the spin-density limits m — > n for n € [0, 1] and m — > (2 — n) for n £ [1, 2]. This 
follows from the si band ground-state occupancies being otherwise finite. 

When the initial ground state has a vanishing av fermion density and thus for finite u the av fermion velocity 
vanishes at the limiting momentum values q = ±g Q1 ,, its u dependence has a singular behavior at u = 0, becoming 
finite. Such av fermion velocity behavior is confirmed by inspection of the slopes at the av bands limiting momentum 
values of the si, s2, 77I , and r]2 fermion energy dispersions plotted in Figs. [21 El SI and[SJ 

An important consequence of the invariances considered here refers to particular case of one sv =/= si fermion (and 
r\v fermion) created onto a S s = (and 5 r; = 0) ground state. For instance, for an excited state with number values 
N au — 1 and N av * — for v' > v branches one confirms from the use of Eqs. (I101[) and (|102p that N aav — 1. Thus the 
av band consists of a single momentum, q = ±q au — 0. Furthermore, such an object energy reads e sv — (a — s) and 
e. qv = 2vfi° (a — 77), consistently with Eq. (|1T8|) . In that single-momentum and single-energy state the av fermion 
has vanishing velocity. Hence it is invariant under the electron - rotated-electron unitary transformation. The same 
applies to a finite number of sv 7^ si fermions (and rjv fermions) created onto an initial S s = (and S n — 0) ground 
state. 
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I. The simplest excitations and the distinct types of elementary objects 

The operator formulation considered in this paper accounts for the representations of the model global SO(3) C3> 
SO(3) <g> U(l) — [SO (4) (g> U(l)]/Z 2 symmetry algebra. Here we discuss the relation of the simplest excitations to the 
representations of that symmetry algebra and the properties of the distinct types of elementary objects. 

Consistent with the results of Appendix 151 for a ground state with densities n € [0, 1] and m £ [0,n] the numbers 
of our formulation objects read, 

M a = M h s ° + 1/2 = N = n N a ; = 2N X ; M™ 1/2 = JV t - N x = m N a , 

M v = M™ = M™ 1/2 = N a - N = (1 - n) N a , 

N c — N ; N al = Ny, N^ = N t -N±=mN a . (180) 

Since the equality Nj? = [N a — N c ] holds, here and below no values of are given. For the present general ground 
state, the values of the following numbers vanish: M^ n _ 1 ^ 21 M^°, -^"™x/2j anc ^ se ^ °f au av fermion numbers {N av } 
except that of the si fermions. For it the si band momenta are occupied for q E [— qp s i, qFsi] = [— fcF4.jfc.F4.] and 
unoccupied for |g| e [qfb1,Qsi] = [fcF4,fcFt]- The c band has the same occupancies as the m = ground state whose 
object numbers are provided in Eq. (|168[) . 



1. The elementary objects within spin excitations 

Within the present representation the excitations associated with the si band holes refer to the spin degrees 
of freedom. However, alike the corresponding composite si fermions, the si fermion holes are spin-neutral objects. 

Consistent, in each unbound 77-spinon and unbound spinon subspace, the C N SV configurations of the N£ v sv fermion 
holes are accounted for in the number J\f S i ng iet(S s , M s ) of spin-singlet configurations, which is given in Eq. (|125[) for 
a = s. Such spin-singlet configurations arc different and independent from the (2S S + 1) = (M" n + 1) spin-multiplet 
configurations of the M" n unbound spin-1/2 spinons. Consistent with the expression given in Eq. (1851) . the spinon 
transformation laws under the electron - rotated-electron unitary transformation performed by the BA solution reveal 
that the spin S s and spin projection S^ 3 values of an energy eigenstate are fully determined by the occupancies of 
such M" n unbound spinons. 

Indeed, an important physical consequence of such transformation laws is the independence of the occupancy 
configurations of unbound spinons and bound spinons, respectively. The former and the latter are invariant and 
are not invariant under that transformation, respectively. One has found that the bound spinons exist within spin- 
neutral 2^-spinon sv fermions. Here v = 1, ...,00 is the number of spinons pairs. The latter objects live on an sv 
effective lattice with a well-defined number, N su , of occupied sites (sv fermions) and, AT™,, of unoccupied sites (sv 
fermion holes). The conjugate coordinates of such sites spatial variables are the sv band N asu = [N s „ + N^] discrete 
momentum values. A number N s „ of those are occupied (sv fermions) whereas N^ v are unoccupied (sv fermion holes). 
In units of 2n/N a , such discrete momentum values are the BA solution quantum numbers, Eq. (|113|) . 

The counting of the model global symmetry algebra representations of Section IV Bl refers to independent represen- 
tations. Their total number equals the Hilbert-space dimension A Na , Eq. (|124[) . Therefore, our above analysis of 
the J\f(S s ,M s ) spin configurations of each unbound 77-spinon and unbound spinon subspace explicitly confirms that 
within the present representation the M™ unbound spin-1/2 spinons are different from the Ng U sv fermion holes. 
The unbound spin-1/2 spinons play the passive role of unoccupied sites of both the spin effective lattice and the 
corresponding sv effective lattices. Thus the objects that move around in the spin effective lattice ring of length 
L = N a ^ a s = N a a are the sv fermions. In that motion they use all the M" ra = 2^ unbound spinons and the 2(v' — v) 
bound spinons in the sets of 2(v' — v) sites of the spin effective lattice out of the 2v' sites occupied by each sv 1 fermion 
with 2v' > 2v bound spinons. Note however that the holes in the corresponding BA sv momentum bands exist in their 
own right and must not be treated as being the M" n = 2S S unbound spinons and Yl v '=v+i ~ v )N S v' unbound 
spinons that play the role of unoccupied sites. The M" n = 2S S unbound spinons have their independent occupancies 
and the sv' fermion with 2v' > 2v bound spinons their own independent sv' momentum bands. Consistent, the sv 
effective lattice object motion can alternatively be described in terms of sv fermions and sv fermion holes, which in 
general have a finite velocity. 

That the unbound spinons and sv fermion holes are independent objects that describe different degrees of freedom 
of the spin SU(2) symmetry algebra is confirmed by the studies of the second paper on the states spin currents. 
Consistent with the passive role played by the unbound spinons in the sv fermion motion, such spin currents are 
found in that paper to be carried by the sv fermion holes or sv fermions. On the other hand, the effects on such 
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currents of the spin transverse fluctuations associated with spin-flipping the unbound spinons are fully controlled by 
such spin-1/2 objects, which exist on their own right. 

Within the N a — > oo limit, a sv fermion moving around in the lattice ring of length L = N a a is exactly described 
by the same sv fermion moving around in its sv effective lattice ring of length L = N asv a sl , = N a a and spacing 
a sl/ = L/N aav . Inspection of Eq. (|I02I) reveals that such a sv fermion does not "see" the spin-effective lattice sites 
occupied by sv" fermions with a smaller number v" < v of spinon pairs. On the other hand, it uses the M"™ sites of 
that effective lattice occupied by unbound spinons as sv effective lattice unoccupied sites. Furthermore, it uses as well 
as sv effective lattice unoccupied sites the 2(v' — v) sites of the former effective lattice out of the 2v' sites occupied 
by each sv 1 fermion with v' > v spinon pairs. 

Within the present representation this does not imply however that M" n out of the N^ v > M"" unoccupied sites of 
the sv effective lattice should be identified with the M"™ unbound spin-1/2 spinons. The remaining \N% V — M^ n ] > 
unoccupied sites of the sv effective lattice neither should be identified with the sets of 2(v' — v) spin-1/2 spinons of 
each sv 1 fermion of v' > v branches with finite occupancy in the state. Those are also used by the sv fermion as 
unoccupied sites of its sv effective lattice. 

For instance, let us consider an excited state of the ground state whose numbers are provided in Eq. (|180[) . That 
excited has one four-spinon composite s2 fermion. Its spin degrees of freedom object numbers are, 

M s = nN a = [Af™ + 2N sl + 4N s2 ] ; M s un =mN a , 
N sl = [(n - m) N a /2 - 2] ; N s2 = 1 , 

Njji = [Air +2] = [m N a + 2] ; N% = M™ = m N a . (181) 

The excited-state s2 effective lattice has N ae2 = [A^+A^] = [mN a + l] sites. Hence the s2 fermion can move along 
its s2 effective lattice by interchanging position with each of the M" n = m N a sites occupied by unbound spinons in 
the underlying spin effective lattice. The s2 fermion does not "see" the 2N s i spin-effective lattice sites occupied by 
the si fermions. On the other hand, the AT s i = [(n — m) N a /2 — 2] si fermions move along their si effective lattice by 
interchanging position with both each of the sites occupied by the Af "™ = m N a unbound spinons in the underlying 
spin effective lattice and two out of the four sites occupied by the spinons of the composite s2 fermion. 

However, within the present representation neither M™" out of the = [M"™ + 2] si fermion holes nor the 
Ng 2 = M™ n s2 fermion holes should be identified with the M" n unbound spin-1/2 spinons. Those exist in their own 
right, their occupancies generating in this case (2S S + 1) = (M"" + 1) = (m N a + 1) spin-multiplet configurations. 
Consistent, the occupancy configurations of the = [M™+2] si fermion holes and N^ 2 = [M™+2] s2 fer mion holes 

generate different and independent C^ sl = C™^° +2 and C^ b2 = C™^ a +1 — m N a + 1 spin-singlet configurations, 
respectively. All such arguments apply as well to the spin-singlet excited states whose deviations are provided in Eq. 
(11691) . which are a particular case of the excited states under consideration. 

Moreover, for an excited state whose generation from the general ground state with numbers provided in Eq. (| 180[) 
breaks the spinon pair of one si fermion, giving rise to two unbound spinons, the spin degrees of freedom object 
numbers are, 

M s = nN a = [M" n + 2iV sl ] ; M"" = [m N a + 2} , 

JV s i = [(n - m) N a /2 - 1] ; AT* = M"" = [m N a + 2} . (182) 

Here the values of the densities n and m are those of the initial ground state. The A^\ = M"™ si fermion holes of such 
an excited state should not be identified with the M" n unbound spin-1/2 spinons, yet their number exactly equals 
that of such objects. Those exist in their own right, their occupancies generating in this case (M™ 1 + 1) = (m N a + 3) 
spin-multiplet configurations. On the other hand, the occupancy configurations of the N' s \ = M™ 1 si fermion holes 

generate different and independent C^ sl = C™^^ 2 spin-singlet configurations. The same applies to the spin-triplet 

states whose deviations are given in Eq. (|170|) . which are a particular case of the excited states considered here. In 
that case there are (M"" + 1) = 3 spin-triplet configurations. Each of these three configurations involve M"" = 2 spin 
1/2 unbound spinons. Within our representation they are described by two unbound spinons with spin projection 
+1/2, one unbound spinon with spin projection +1/2 and one unbound spinon with spin projection —1/2, and two 
unbound spinons with spin projection —1/2, respectively. 

2. The elementary objects within excitations involving the three global symmetries 

Within the present representation 2S 1 ' 1 = [N a — 25 c ] is the number of rotated-electron doubly occupied plus unoccu- 
pied sites. The independent occupancy configurations of the corresponding M v = [N a — 2S C ] 77-spin 1/2 ?7-spinons and 
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N' c l = [N a — 2S C ] c fermion holes are associated with two degrees of freedom of the rotated-electron doubly-occupied 
sites and unoccupied sites occupancy configurations. Those refer to representations of the 77-spin SU{2) symmetry 
and c hidden U{\) symmetry algebras, respectively. 

As discussed in Section IV Bl a unbound 77-spinon and unbound spinon subspace is spanned by all states with 
fixed values for the numbers S v , M n = [N a — 25c], and = [N a — 2S C ]. For each such a subspace the number 
of configurations of the M v = [N a — 2S C ] ?/-spin-l/2 77-spinons is AT(S V , M v ), Eq. (II 2 1 [) . On the other hand, the 
number of configurations of the N£ = [N a — 2S C ] of c fermion holes equals the dimension d c , Eq. ()103j) . which may be 

rewritten as d c = C N " = C N a . The dimensions Af(S v , M v ) and d c , which appear in the overall subspace dimension 
of Eq. ()120j) . are found in Section IV Bl to correspond to completely different and independent configurations. They 
refer to the number of subspace representations of the 77-spin SU(2) symmetry algebra and c hidden U(l) symmetry 
algebra, respectively. The model global symmetry imposes that both the number M n = [N a — 2S C ] of 77-spin-l/2 
77-spinons and N% = [N a — 2S C ] of c fermion holes equal that [N a — 2S C ] of rotated-electron doubly plus unoccupied 
sites. Nonetheless within the present representation they are different and independent objects in their own right. 

N N h 

Furthermore, consistent with the form of the number of subspace configurations d c — C N " — C N " , the c fermions and 
c fermion holes alike are 77-spin-less and spin-less objects. 

We illustrate the issue discussed here in the case of the simplest one-electron addition excited states of the \i — 0, 
7i = 1, and m = absolute ground state. The generation from it of such one-electron addition excited states involves 
changes in the degrees of freedom associated with the c hidden U(l) symmetry and two SU(2) symmetries. Under 
such an excitation, one si fermion spinon pair is broken. The created electron combines with one c fermion and that 
of the two spinons of the broken si fermion whose spin projection is opposite to its own. Such a process gives rise to 
one rotated-electron doubly occupied site. Its 77-spin SU{2) symmetry degrees of freedom are described by the —1/2 
77-spinon created under the one-electron addition excitation. The c hidden U(l) symmetry degrees of freedom account 
for the new rotated-electron doubly occupied site through the creation of one c fermion hole. Its available momentum 
values q G [—if, 7r] correspond to c band occupancy configurations that generate representations of the c hidden U(l) 
symmetry algebra. They are associated with the dimension d c in the excited subspace dimension given in Eq. (I120p . 
Such finite- velocity c fermion hole should not be identified with the —1/2 77-spinon also created under the transition 
to the excited state. There is also a spinon left behind by the si fermion spinon-pair breaking. It has the same spin 
projection as that of the created electron. Such a unpaired spinon becomes a unbound spinon in the excited state, 
which does not contribute to the state momentum. It should not be identified with the finite- velocity si fermion hole 
of momentum q' G [— 7r/2, 7r/2], which is also created as a result of the si fermion annihilation. Such unbound spinon 
and si fermion hole describe different spin degrees of freedom. 

VI. QUANTUM PROBLEMS DESCRIBED BY THE ID HUBBARD MODEL IN SUBSPACES OF 

PHYSICAL INTEREST 

A. Symmetries and physically important subspaces 

The combination of the present operator formulation with the PDT of Refs. [8l-[l0j reveals that for the metallic 
phase of the ID Hubbard model nearly the whole spectral weight associated with one- and two-electron excitations is 
generated by transitions from the ground state to both LHB excited states and first UHB excited states, as defined 
in Section TV Fl For electronic densities n < 1 and n > 1, the LHB (and first UHB) subspace is spanned by energy 
eigenstates with no (and one) rotated-electron doubly occupied site and unoccupied site, respectively. 

The metallic-phase LHB subspace is a large subspace, which is spanned by 3^" energy eigenstates. Within it one 
has that = so that the eigenvalue 2S C of the c hidden ?7(1) symmetry generator and the 77-spin S v are trivially 
related as S v — [N a — 2S c ]/2. For the ID Hubbard model in that subspace the 77-spin SU(2) degrees of freedom 
play a minor role in charge and spin excitations that preserve the c fermion number N c — 2S C and thus the 77-spin 
Srj- Indeed, those of the energy eigenstates in that subspace that have the same 77-spin value S v have as well exactly 
the same 77-spinon occupancy configuration. It refers to M n = M^ n +1 p = 2S V (and M v — M'^ n _ 1 p — 2S V ) unbound 
77-spinons, all with the same 77-spin projection. 

Note however that excitations within that subspace as for instance one-electron addition and removal excitations 
under which the deviation 5M V = SM^ 1 ^ 2 = 2SS V is finite involve a change in such 77-spin SU{2) symmetry degrees 
of freedom occupancies. On the other hand, the charge and spin excitations within it conserve the c fermion and 
77-spinon numbers. Thus for them the quantum problem associated with the ID Hubbard model in that subspace has 
effectively a charge U(l) symmetry and a spin SU(2) symmetry. Our results reveal that that charge U{1) symmetry 
is the c hidden U(l) symmetry, which is independent of the 77-spin SU(2) symmetry and 77-spin U(l) symmetry 
contained within it. 
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B. The quantum problem described by the ID Hubbard model in the subspace with no rotated-electron 

doubly and unoccupied sites 

Another important subspace is that spanned by half-filling energy eigenstates with both no rotated-electron doubly 
occupied sites and no rotated-electron unoccupied sites. The 1 D Hubbard model in that subspace is the issue mainly 
addressed in this section. In it, that model becomes a spin-only problem. In the u 3> I limit it is the only finite- 
energy half-filling subspace. The relation within the m> 1 limit of the present representation spin- 1/2 spinons to the 
excitations of the related ID spin- 1/2 isotropic Heisenberg antiferromagnetic model [26[ is a problem briefly discussed 
in this section. 

The subspace with no rotated-electron doubly-occupied and unoccupied sites contains the absolute ground state 
of vanishing chemical potential, fj, = 0, electronic density n = 1, and vanishing spin density m = 0. We call it NDU 
subspace, which stands for no (N) doubly-occupied (D) and no unoccupied (U) sites subspace. An energy scale that 
plays an important role in the energy spectrum of such a subspace is the minimum magnitude of the excitation energy 
A D rot for creation of a number D rot = M v of rotated-electron doubly occupied sites and rotated-electron unoccupied 
sites onto a n = 1 and m = ground state relative to its energy level. Such a ground-state energy level coincides with 
the zero-energy level of the Hamiltonian H of Eq. (JXJ) . Consistent with Eq. ([55)1 and Eq. (II 75[) for n = 1, it reads, 

minAfl„, = (J?M V + (2a v ) pM™ v , at n = 1 and n £ [-/jP, fi°] , 

<r„=±l/2 

= n° M v , at n = 1 and fi = . (183) 

It follows from Eq. (|183|) that for one-electron and two-electron excited states the NDU subspace is for u > the 
only Hilbert-space subspace for an excitation-energy window uj 6 [0,/i°] and u> € [0, 2/jP], respectively. Indeed within 
that excitation energy range there is neither rotated-electron doubly occupied sites nor rotated-electron unoccupied 
sites for the Hubbard model given in Eq. (|2|) at chemical potential fi = 0. The energy eigenstates that span such 
a subspace are superpositions of a subset of local-rotated-electron states whose N a sites are all singly occupied by 
rotated electrons. Moreover, the c fermion momentum band of such energy eigenstates is full. 

The energy scale fiP, which defines the excitation-energy range of the present quantum problem, is one half the 
Mott-Hubbard gap 2fi°. Its exact u dependence is given in Eq. (|C3|) of Appendix [CJ whose u 1 corresponding 
limiting behavior is 2/i° ss [U — it], as given in Eq. (|C4[) of that Appendix. Hence for u —> oo one has that 2/jP — > oo, 
so that the NDU subspace is the only Hilbert-space subspace for all finite values of the excitation energy. 

For such a subspace there is no 77-spin effective lattice, whereas the spin effective lattice has N a sites and is identical 
to the original lattice. The spin degrees of freedom of the N = N a rotated electrons that singly occupy the N a sites 
are described by M s+1 / 2 — spin-up and M s _ 1 / 2 = N± = [N a — N^} spin-down spinons. The total number of 
spinons is thus fixed and given by M s = [M S;+1 / 2 +M S x/2] = -^o = N. The numbers M s and M Si ±i/2 can be written 
as the sum of two terms, which refer to unbound and bound spinons, respectively, 

00 

M s = N = M^ n + M b s ° = 2S S + 2Y^ vN sv , 

v=\ 

00 

M Si±1/2 = M™ ±1/2 + M h s °/2 = [S s T S^} + vN sv . (184) 

Here M«" = 2S S , M^ n ±1/2 = [S s T S**}, and M b s ° = 2^ = i » N s»- 

Since there is no ?7-spin effective lattice and the c effective lattice is full, the important dimension is that given in 
Eq. flUT} for a = s, 

M(S S ,M S ) = (2S s + l)M smg i et (S s ,M s ) 

= (N - M b s ° + 1) N slnglet (N/2 - M b s °/2, N) , (185) 

where 

A/7 9 M\ , - r M °/ 2 ~ S ° r M./2-S.-l 
■'V \o s , ivi s j 8ln gi e t — <^jvf s M s 

= cf>° /2 -cf°° /2 -\ (186) 

and we used that 25*8 = [N — M bo ]. The number of bound-spinon pairs M h s ° j2 may have the values M h s ° /2 = 
0, 1, 2, N/2, so that the Hilbert-space dimension of the quantum problem corresponding to the half-filled ID Hub- 
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bard model in the NDU subspace is, 

(N~M b s + l){cf°° /2 -cf°° /2 - 1 }=2 N = 2 N *. (187) 



N/2 

idu ^ ^ 

M|°/ 2 = 



Since in that subspace there is only rotated-electron single occupancy, upon acting onto it, the rjs local quasi-spin 
operators of Eq. (|66|) simplify to, 

Qj = «j> = 4,t ; ^ + = S J> = ( S J>) f ; ?T = S i,s = _ ( ft J.t ~ ™3\i)/ 2 ■ ( 188 ) 

Furthermore, again accounting for the lack of both rotated-electron doubly occupied sites and unoccupied sites, 
upon writing the Hamiltonian of Eq. (fT5|) in the NDU subspace, one finds that all its terms of odd order vanish and 
the terms of even order simplify to, 



= UV C 



H {2) = -— T_iT_ 



t 2 



u 

H (4) = ^ [f_if+if_if+i - l x f % x - f-iT 2 f +1 ] . (189) 

Within the NDU subspace, the contributing Hamiltonian terms of Eq. (| 1 89[) belonging to the general term T^ 1 Tl 1 
are of the form T—i-jjiT-^i-jj/T—i-ji/jniT^i-jujin where both nearest-neighboring sites of indices j and j' are different 
from the two nearest-neighboring sites of indices j" and j'". Moreover, in the Hamiltonian term containing TLiTgT+i, 
the contributing terms are of the forms T-i ; j t j'To-j t j»To ; j t j»T+i-j» j j///T+i- l j t ji or T-i-j^'Tq^' ,j'"Tb ; j/ t f»T+i-j» j/r/T+i-jj/ 
where j" ^ j 1 and j'" ^ j, respectively. In addition, after some straightforward algebra one finds that within the 
NDU subspace all products of operators F T j j> of Eq. (f7Tj) in the corresponding products of rotated kinetic operators 
T 7 = —[t/U]^2/j ji\ Fy>j,f Q\j\;jJ' °f -^l- GO]) appearing in the Hamiltonian terms provided in Eq. (|189l) can be 
expressed only in terms of local operators and q c y . In the NDU subspace one can then replace such operators by 
their eigenvalue, which is 1. Thus the Hamiltonian terms of Eq. (| 1891) can be expressed in terms of only spinon 
operators. Moreover, the zeroth-order term becomes a mere constant, flw = [t//4] N a , and may be ignored. 

In addition, by manipulations of the commutators of Eq. (|72l) . we find that for the Hamiltonian in the NDU 
subspace the term of sixth order can also be expressed only in terms of spinon operators. In the present ID case 
the problem simplifies, since many contributions vanish, and one expresses the commutators of Eq. (|72| directly in 
terms such operators. Fortunately, at n — 1 there is no need of deriving the Hamiltonian sixth order term expression 
explicitly in terms of the three rotated kinetic operators, in contrast to for densities n =/= I. That is a more involved 
problem that we do not address here. 

In order to make our notation compact, we use in the following Hamiltonian expression the spinon operator vector 
§j iS whose operator components are sj 1 ,,, sj 2 s , and sj 3 s . Upon expressing the Hamiltonian operators and of 

Eq. (|189[) along with the Hamiltonian operator in terms of spinon operators we find after some algebra, 

Hsymm = + + + ... , 



k{2) = -jf E^-^V^+m), 



u 
u 3 



3=1 



[-4(1 -4 4s • S 3+1 , s ) + (1 - 4 4, • s j+2 ,s) 



3=1 



t 6 r 

H (6) = -jjs 2 [ 15 ( : ~ 4 4« ' - 6 (! - 4 4« ■ *i+2,«) + (1 - 4 4* • *3+3,.) 



3=1 



+ 7/5 E P ~ 4 % s ' 5 3+3,«)( 1 - 4 4« ' S 3+2, s ) ~ (1 -45^ ■ S i+ 2, s )(l ~ 4 4s ■ Sj+3, S )J ■ (190) 
3=1 

For the corresponding spin excited states, such expressions are valid at fi = and excitation energy w < 2/x°. Note 
that in the NDU subspace the expression of the Hamiltonian operators and do not include terms multiplying 
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the non-universaly valued real-number parameter 6. Therefore, their expressiona are fully known and given in Eq. 

po|i . 

The fi = 0, n = 1, m = absolute ground-state energy is to t 6 /U 5 order readily obtained from the general ground- 
state energy expansion derived in Ref. (55j | for n G [0, 1] and m = 0. By choosing n = 1 in that general expansion we 
find, 

E° GS /N a = ~ 4 In 2 + ^ 9 C(3) - ^ 75 C(5) + ... , (191) 

where £(3) and £(5) are Riemann Zeta functions. This ground-state energy refers to the Hamiltonian expression to 
t 6 /U 5 order given in Eq. (|190l) . 



VII. THE ID HUBBARD MODEL IN THE NDU SUBSPACE FOR u > 1 

For !t> 1 the leading-order term of the Hamiltonian expression provided in Eq. (I190p becomes dominant, so 
that in the NDU subspace that Hamiltonian approximately reads, 




(192) 



where N = N a . Within the u ^> 1 limit the rotated electrons become electrons, so that the local spinon operators §j :S 
defined in Eqs. and (|57|) that appear in Eq. (|190[) were here replaced by the corresponding unrotated local spin 
operators, §j tS , provided in Eq. (|42l) . Consistent with the known u ^> 1 relation of the half- filled ID Hubbard model 
to the ID isotropic Heisenberg antiferromagnetic model of spins 1/2, for exchange constant J = 4t 2 /U the expression 
given in Eq. (|192p is exactly that of the latter model. 

Since in the u 1 limit the N rotated electrons become the N electrons, the M s = N spin-1/2 spinons of the 
present formulation become the spins of each of the N electrons. At N = N a the NDU subspace becomes the only 
finite-energy subspace as u — > oo. Within it the spin and charge degrees of freedom are decoupled, so that the spins 
of each of the half-filled ID Hubbard model N electrons become the N spins of the ID spin-1/2 isotropic Heisenberg 
antiferromagnetic model. Thus at half filling and for u — > oo the M s = N spin-1/2 spinons of the present formulation 
become such ./V spin-1/2 spins. 

This reveals that in such a limit the spin-1/2 spinons of the present electron and rotated-electron vacuum normal- 
ordered formulation and the spin-1/2 spin waves of the ground-state normal-ordered description of Ref. [26[ are 
completely distinct quantum objects. Indeed, the S v = 0; S s = 0; 2S C = N a absolute ground state of the ID Hubbard 
model is populated by N — N a spin-1/2 spinons, as defined in the present formulation, whereas it contains no spin- 
1/2 spinons, as defined in Refs. [29l. l30jj . Alike the spin-1/2 spin waves of Ref. (2|| are excitations of the zero-spin 
ground state of the ID isotropic Heisenberg antiferromagnetic model, the spin-1/2 spinons of Refs. [2!j|3(| have been 
constructed inherently to be excitations of the = 0; S s = 0; 2S C = N a absolute ground state of the ID Hubbard 
model. 

The relation for u ^> 1 of the TV spin-1/2 spinons of the present formulation to the spin-1/2 spin waves of Ref. [26[ 
is in terms of electrons ex actly the same as that for u > of the spin-1/2 spinons of the present formulation to the 
spin-1/2 spinons of Refs. [29l . |30| in terms of rotated electrons. Thus we forward the reader to the second paper, 
Ref. [18j |. for further information on the relation between the present spin-1/2 spinons and those of the ground-state 
normal-ordered traditional spinon descriptions. 

Importantly, in the second paper it is confirmed that there is no contradiction whatsoever between the spin-1/2 
spinons of the present formulation and those of Refs. 28-3(3] ■ The former and the latter spinon definitions are 
normal ordered relative to the electron and rotated-electron vacuum and an initial S s = ground state, respectively. 
The two types of spin-1/2 spinons are confirmed to be distinct elementary objects associated with uniquely related 
scattering-state basis choices. The possibility of such different choices stems from the degeneracy of the excited energy 
eigenstates of the S s — ground state with the same number of spin-1/2 spinons as defined in Refs. [28-30]. 



VIII. CONCLUDING REMARKS 



Since the BA solution of the ID Hubbard model is exact, most results of studies by means of that solution are 
correct independently of the interpretation of which elementary objects are involved in its excitations. However, the 
relation of the operators of the elementary objects to those of the original electrons is an issue of highest importance 
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for the physical interpretation of such results. Unfortunately, the unique definition of the elementary objects operators 
and their relation to the electron creation and annihilated operators remains for most object representations an open 
problem. That unsolved problem is one of the main motivations of this paper. Indeed, the charge-spin separation 
and corresponding interpretations in terms of exotic objects such as holons and spinons play a major role in the 
interpretation and description of the low-dimensional correlated systems physics. 

The operator formulation introduced in this paper goes beyond the BA solution subspace. It has been constructed 
inherently to spinon, ?y-spinon, and c fermion occupancy configurations generating all 4 Na energy eigenstates. Con- 
cerning the general problem of the ID Hubbard model in its full Hilbert space, our results confirm that there occurs 
a three degrees of freedom separation associated with the two previously known spin SU (2) and fy-spin SU (2) sym- 
metries plus the c hidden U(l) symmetry found in Ref. [19j. From such three symmetries three types of elementary 
objects follow: The spin-1/2 spinons, the ?7-spin-l/2 7/-spinons, and the spin-less and 77-spin-less c fermions, respec- 
tively. The corresponding c fermion operators and rjs quasi-spin operators whose expressions are provided in Eqs. (|61[) 
and (|67p , respectively, are mapped from the rotated-electron operators by an exact local transformation that does not 
introduce constraints. The same applies to the related spinon and ?7-spinon operators defined in terms of the local c 
fermion number operators and rjs quasi-spin operators in Eq. (|66p. On the other hand, the electron - rotated-electron 
unitary operator is defined in Section Til Bl it terms of its A Na x 4^° matrix elements given in Eqs. (f5T|) and (|52l) . 
Profiting from the interplay of the exact BA solution with the model global symmetry algebra representations, we 
have found which of these matrix elements vanish. Those that are finite have been expressed in terms of exact BA 
amplitudes, which we have extended to the non-LWSs generated from the Bethe states. 

The transformation laws under the electron - rotated-electron unitary transformation of the M s = = 2S C spinons 
(and M n = N% = [N a - 25 c ] ^-spinons) imply that M s bo = [25 c - 2S S ] of such spinons (and M%° = [N a - 2S C - 25„] of 
3nt 
Mt 

rbo 
1 1 



such 77-spinons) are bound (and anti-bound) within spin-neutral (and ^-spin-neutral) composite sv fermions (and rjv 
fermions). Although the physical meaning of the eigenvalue 2S C of the c hidden C/(l) symmetry and of the numbers 
Mg° and remained masked by the underlying mathematical formalism, we have shown that they appear in 



the exact BA thermodynamic equations introduced in Ref. [IJ. As confirmed in this paper, that solution explicitly 
accounts for the model global symmetry and corresponding transformation laws under the electron - rotated-electron 
unitary transformation. This is for instance explicit in the form that the rotated-electron amplitudes given in Eqs. 
(1551) and (|54"|) have for all 4^° energy eigenstates. 

The model global [50(4) <g> U(l)]/Z 2 = 50(3) ® 50(3) ® U(l) symmetry accounted for by the operator formulation 
introduced in this paper applies as well to the Hubbard model on any bipartite lattice. The studies of Ref. [56[ 
present an approximate operator formulation to the model on the square lattice. Its motivation and starting point is 
that introduced here for the ID Hubbard model, in terms of quantities extracted from the exact BA solution. 

Several properties predicted by the ID Hubbard model were observed in low-dimensional complex materials [5ll57j. 
The investigations presented in Refs. [TJ HI1 confirm that the PDT successfully describes the unusual finite-energy 
spectral features observed by angle-resolved photoelectron spectroscopy in quasi-lD organic metals. The results 
presented here and in the ensuing paper contribute to the further understanding of the microscopic non-perturbative 
mechanisms behind these properties. In the second paper it is shown that the scatterers and scattering centers that 
naturally emerge from the present elementary-objects operator formulation for the ID Hubbard model are the spin-less 
and ry-spin-less c pseudofermions, the spin-neutral composite sv pseudofcrmions, and the 77-spin-neutral composite 
r\v pseudofermions. Those are related to the c fermions, sv fermions, and r/v fermions, respectively, by a unitary 
transformation. For the excited states, it slightly shifts the discrete momentum values of such objects. Otherwise 
the c fermions, sv fermions, and rjv fermions that emerge from our rotated-electron operator formulation and the 
corresponding c pseudofermions, sv pseudofermions, and rjv pseudofermions, respectively, studied in the second paper 
have the same properties. 

The pseudofermion scattering theory and corresponding dressed 5 matrices and phase shifts are shown in the second 
paper to control the spectral- weight distributions of the PDT |8l4ldl| . Specifically, such dressed 5 matrices control the 
momentum, electronic density, spin density, and repulsive interaction dependence of the spectral-function exponents 
associated with the high-energy behavior beyond the linear Luttinger-theory predictions revealed by the PDT for 
the ID Hubbard model and more recently by the methods of Refs. [Tij - flil 17 | for a wider class of integrable and 
non-integrable ID correlated problems. The relation of the present formulation elementary objects to traditional 
spinon and holon descriptions as those of Refs. [28| - |30| is an issue that is also investigated in the second paper, Ref. 
[18| . As mentioned in previous sections, it is found in that paper that there is no contradiction whatsoever between 
such distinct elementary objects choices and corresponding scattering theories. The possibility of different choices of 
both elementary objects and scattering states is found to stem from the degeneracy of the excited energy eigenstates 
that span the subspaces of the corresponding scattering theories. 
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Appendix A: Useful quantities and transformations 



In this Appendix we introduce the known model L-matrix expression. Furthermore, it is confirmed that the 
operational expressions of the rotated-electron operators in terms of rjs quasi-spin operators and c fermion operators, 
Eq. (|69l) . give the correct results upon their application onto the one-site states \j,Q,h c ), \j, fj,, h c ), \j,t,Pc), and 
\j,i,p c ) appearing on the right-hand side of Eq. ([82]) . 

The known expression in terms of electron creation and annihilation operators of the BA inverse-scattering method 
on-site 4x4 L-matrix £y(A) of Eq. (j6U|) is (43j. 



C, 



where 
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ibn 
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[an jt(T - ibrij t0 -] , 
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and 



3,° 



(1 — nj. a ). It is shown in this paper that it can be expressed in terms of the seven generators of the local 



gauge SU(2) <g) SU (2) ® U(l) symmetry of the ID Hubbard model at t = and corresponding spin-less and 7/-spin-less 
fermion operators, as given in Eq. (|60[) . 

The seven local operators obtained by electron rotation of such generators are expressed in Eq. (|67p in terms of the 
•qs quasi-spin operators and c fermion operators defined in this paper. Application of the operational expressions of 
the rotated-electron operators in terms of such rjs quasi-spin operators and c fermion operators provided in Eq. (|69|) 
onto the one-site states \j,Q,h c ), \j, fi, /i c ), \j,t,p c ), and \j,i,p c ) leads to, 



fj,c ( 2 



1 

2 

(-1.1'//.. (~ 



q? 



b',t,Pc) = o, 

|jU,Pc) = fj,c\j,UPc) = \j,u,h c ), 

fl Q - <?f) l/> ©, h c ) = fl li, 0, h c ) = li, t,p c ) , 
~qf ) \j,U,hc)=0, 



V 



(A3) 



and 



4ji,0,M = fi c qf\j,Q,h c ) = (-iy fl c \j,U,h c ) = \j,UPc) 
c]jj, U,h c ) = fl c qf\j,U,h c }=0. 
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Hence for the rotated-electron annihilation operators the result is, 

C*,ttfi V, c ) = fj.e (~ - |j, t,Pc) = fj,c \j, t,Pc) = \j, 0, he) 



and 



gj-.tli.tiA) = (-i)7j, c {\ + ^A \j>n,h e ) = (-iy f]jj,n,h c ) = \j,±, Pc ) , (A5) 



culi.t.Pc) = f 3 ,c.qj \jA,Pc) = 0, 

i,Pc) = fj.c qj \j, i,Pc) = fj,c \j, t,Pc) = \j, 0, h c ) , 

c 3 ,i\j, 0, K) = (-IY /t c qj \j, 0, h c ) = 0, 

c u \j, U, fee) = /t c g- |j, tl, /ic) = /j, c |i, 0, h c ) = \j, t,p c > ■ (A6) 

In the above equations it was used that, as given in Eq. (|67| . the operator gj where Z = 2:3, ± reads gj = si B (and 
gj = ) when the site j is singly occupied by rotated electrons (and unoccupied or doubly occupied by rotated 
electrons.) 

Analysis of Eqs. (|A3[) - (IA6[) confirms that the correct final states are obtained. 



Appendix B: Unbound spinon and unbound 77-spinon energies and ground-state occupancies 

In this Appendix the energies of the ±1/2 unbound 77-spinons are explicitly derived. A similar analysis provides 
those of the ±1/2 unbound spinons. The BA solution subspace is smaller than and contained in the ID Hubbard 
model full Hilbert space. Thus analysis of the energy spectra of the Bethe states alone is not enough to clarify the 
problem of the ground-state occupancies. Here we combine the model global symmetry with the BA solution to also 
address that problem. 

For simplicity, to start with we consider m = ground states of the Hamiltonian of Eq. (TTJ at zero magnetic field, 

H = Hsymm - ix 25? = H symm -fi(N-N a ), (Bl) 

where H S ymm is given in Eq. ([2]). 

For u > we denote by Eq S and Eqs the ground-state energy eigenvalues of the Hamiltonians H S ymm and H given 
in Eqs. @ and (|B1|) . respectively, 

^gs(^v) = {GS\H symm \GS) , 
E GS (S V ,S^) = (GS\H\GS)=E GS (S ri )-(i(N-N a ) = E° GS (S, 1 )+fixN a . (B2) 

Here x = (1 — ri) denotes the hole concentration. 

The Hamiltonian provided in Eq. (IB1[) describes the same quantum problem as the Hamiltonian H S ymm, Eq. ([2J, 
but refers to the grand-canonical ensemble. Since the latter Hamiltonian commutes with — fi2S^ 3 , one has that the 
two Hamiltonians commute, [H symmi H] = 0. Thus they have the same energy eigenstates. Only the corresponding 
energy eigenvalues are in general different. Those of H arc shifted relative to the energy eigenvalues of H symm by 
—fj, (N — N a ) = fix N a . Hence the following relation holds for all 4 Na energy eigenstates \^i oi i & , u ), 

(*; ,; A ,«l^l*io,(A,«) = (^i ,iA,u\H symm \^i ,i^,u) - fJ-(N - N a ) 

= (^l ,l^,u\Hsymm\^l ,l & ,u)+fJ'XN a . (B3) 

We start by assuming that for electronic densities n < 1 (and n > 1) the ground state is a LWS (and HWS) of the 
77-spin symmetry algebra. We then use symmetry arguments to confirm that this is so. The zero-energy level of the 
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Hamiltonian H, Eq. (IB1|) . corresponds to the initial ground-state energy. For that Hamiltonian the set of 2S V tower 
energy eigenstates generated by application onto the n < 1 ground state of the off-diagonal generator of the 77-spin 
algebra SJ. a number of times n v — 1,...,25^ have energy eingenvalue (^i ,i^, u \Hsymrn\^>i ,i/s,u) + 2n v \fi\. Thus out 
of the corresponding tower of 2S V + 1 energy eigenstates, that corresponding to n n — has lowest energy. 

However, to confirm that for n < 1 the ground state is a 77-spin algebra LWS we must consider a n < 1 ground state 
other than that belonging to the 77-spin tower. Specifically, n < 1 ground states with the same electron numbers as 
the tower states of energy (\E , i ,i A , u |i? S i/mm| x I / ; ,i A ,ti) + 2n n \[i\. Hence, out of the such tower states, only those having 
electronic densities 77, < 1 and thus n n < S n arc of interest. The 77 < 1 ground states under consideration have r\ 
spin and 77-spin projection {S^ — SS^) and — (S^ — SS^), respectively, where SSr) = n v . Here S v is the 77-spin value 
of the LWSs of the 77-spin tower and n v < S v . Provided that the inequality < 5S n < S v holds, one has in the 
thermodynamic limit that the two ground-state energy eigenvalues under consideration of the Hamiltonian H of Eq. 
(|Bip are such that, 

E GS (S V , -S v ) = E GS (S V - 6S V , -S n + SS V ) . (B4) 

This relation follows from in the thermodynamic limit and for 77 7^ 1 the chemical potential dependence on n associated 
with the ranges provided in Eq. (|Cip of Appendix [Cl being a smooth continuous function. On the one hand, as given 
in Eq. (|B4|) . the ground state energies Eqs^S^ — SS V , —S„ + 5S V ) and EQs(S n , —S v ) have the same value. On the 
other hand, the non-LWS generated from the ground state of energy Ecsi^ri, —S v ) have higher energy than that 
state. Therefore, for the Hamiltonian H of Eq. (|B1[) the energy eigenvalue E(S, V —S n + n v ) with 77^ = SS^ of the 
energy eigenstate belonging to such a tower has the same number of electrons N = N a — 2S V + 2n v as the 77 < 1 
ground states of 77 spin and 77-spin projection (S v — SS V ) and — (S^ — 5S. q ), respectively. The former state has larger 
energy than those ground states. Specifically, the following relation holds, 

E(S n , -S v + n n ) = E GS {S ri - n n , -S v + SS^) + 2?7^|/j| , SS n = n ri . (B5) 

Note that for SS V — n v — the two energy eigenstates whose energies obey this relation are either the same state or 
two degenerate LWSs of the 77-spin algebra. Since for 77^ > one has that Eas{S v — SS V , —S V +5S V ) < E(S n , — S' ?? +77 ?? ), 
one finds full consistency with for n < 1 the ground state being a LWS of the 77-spin algebra. Similar arguments 
confirm that n > 1 ground states are HWSs of the 77-spin SU(2) algebra. 

Concerning the ground-state occupancies in terms of the objects or our operator representation, we note that the 
unitary transformation that relates electrons to rotated electrons preserves their number. Thus the number of rotated 
electrons equals that of electrons. Given that each one-site rotated-electron occupancy has two degrees of freedom, 
associated with the c hidden U(l) symmetry and the 77s quasi-spin SU(2) symmetry, respectively, the rotated-electron 
(and electron) numbers obey several exact relations. The following, 

N a - N = M™ 1/2 - M«"_ 1/2 = M rh+l/2 - M Vt _ 1/2 , 

AT t - JV4. = 1/2 - 1/2 = M S)+1/2 - M s ,_ 1/2 , (B6) 

refer to the 77-spin and spin degrees of freedom, respectively. Moreover, the relations, 

N = N c + M^° + 2M^ n _ 1/2l 

= M s + M*° + 2M™_ 1/2 , (B7) 

refer to the c hidden U(l) symmetry and 77-spin degrees of freedom and to the spin and 77-spin degrees of freedom, 
respectively. The same applies to the following relations, 

N h = N c + M*° + 2M™ 1/2 , 

= M s + M b v ° + 2M; n +1/2 , (B8) 

where N h — (2N a — N) is the number of holes. (Such holes are different from those associated with the hole 
concentration x = (N a — N)/N a . They are such that an unoccupied site has two holes, a singly occupied site has one 
hole, and one doubly occupied site has no holes.) 

Above we have found that the 77-spin SU{2) algebra implies that n ^ 1 ground states are LWSs and HWSs of 
the 77-spin algebra for 77 < 1 and 77 > 1, respectively. A 77 = 1 ground state is then both a LWS and HWS of that 
algebra. Furthermore, in Section IV Gl it is confirmed that the 77-spinons that are not invariant under the electron - 
rotated-electron unitary transformation have an anti-binding character. In the N a — > 00 limit the minimum energy 
for creation onto a n ^ 1 ground state of a 77-spin singlet pair of anti-bound 77-spinons is found in that section to 
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exactly equal that for creation onto it of a pair of unbound 77-spinons of opposite 77-spin projection. Relative to the 
zero-energy level of the Hamiltonian H, Eq. (|B1|) . the latter energy reads (e r) ,+i/2 + Ejj.-i/a) = for n ^ 1, as 
given in Eq. (|90j) . (The validity of the latter equation is confirmed below.) Such a zero-energy level coincides with 
that of the initial 1 ground state. Hence relative to it, creation of anti-bound 77-spinons is a finite-energy process. 

One then concludes that a n < 1 (and n > 1) ground state has no unbound —1/2 77-spinons (and no unbound 
+1/2 77-spinons) and no anti-bound ?7-spinons. The number of anti-bound 77-spinons = [N a — 2S C — 2S V ] given 
in Eq. (|86p is, due to symmetry, always an integer even number. This is consistent with the corresponding ?7-spm- 
singlet configuration involving an equal number, Mr°/2 = [N a /2 — S c — S v ], of 77-spinons of 77-spin projections —1/2 
and +1/2, respectively. Hence the above ground-state vanishing occupancies may be rephrased in terms of general 
77-spinon occupancies (unbound and anti-bound) as follows: A n < 1 (and n > 1) ground state has no —1/2 77-spinons 
(and no +1/2 77-spinons). Indeed, both the energy for creation of one pair of unbound ?7-spinons of opposite 77-spin 
projection and the minimum energy for creation of one 77-spin-singlet pair of anti-bound ?7-spinons onto a n < 1 (and 
n > 1) ground state stem from the corresponding —1/2 77-spinon (and +1/2 77-spinon). This is consistent with a —1/2 
?7-spinon (and a +1/2 77-spinon) describing the 77-spin degrees of freedom of a rotated-electron doubly occupied site 
(and rotated-electron unoccupied site), independently of it being a unbound or anti-bound 77-spinon. Hence one may 
again rephrase the above ground-state vanishing occupancies as follows: A n < 1 (and n > 1) ground state has no 
rotated-electron doubly occupied sites (and no rotated-electron unoccupied sites). 

One then finds by combining the above ground-state vanishing occupancies with the general exact relations of Eqs. 
(|B6[) - (|B8|1 that a m = and n < 1 (and n > 1 ) ground state of the u > ID Hubbard model has N c = 2S C = N (and 
N c = 2S C = N h = [2N a — N]) c fermions, an equal number M s — N c of spinons, Nf? = \N a — N\ c fermion holes, and 
an equal number M^ )+1 / 2 = M n = 2S V — N a — N (and M rh _i/ 2 — M v = 2S V = N — N a ) of unbound +1/2 77-spinons 
(and unbound —1/2 77-spinons). This is consistent with a m = and n = 1 ground state having no 77-spinons and 
no c fermion holes and thus no rotated-electron doubly occupied and unoccupied sites. For such a ground state the c 
fermion momentum band is full. 

Similar results are reached for finite spin densities, concerning the ground-state spin degrees of freedom. In the 
corresponding derivations, the energy scale 2/ib\H\ plays the same role as 2|/x| for the 77-spin degrees of freedom. For 
instance, the spin SU (2) symmetry algebra implies that ra/0 ground states are LWSs and HWSs of that algebra for 
to > and m < 0, respectively. A m — ground state is then both a LWS and HWS of the spin algebra. Hence for 
to > (and 777 < 0) the ground state has no unbound —1/2 spinons (and no unbound +1/2 spinons) and at m = it 
has no unbound spinons. The sv fermions were shown in Section [V Gl to involve binding spinon configurations, rather 
than anti-binding. From the interplay of this property with the spectra extracted from the exact BA solution, one 
then finds that ground states have no sv fermions with v > 1 spinon pairs. On the other hand, except in the case 
of fully polarized states, they have a finite population of si fermions. For instance, for the electronic density range 
77 G [0, 1], the spin degrees of freedom of a m > (and 777 < ) u > ground state are described by M s = N spinons, 
M b s ° = 2N i (and M h ° = 2A t ) bound spinons within N sl = N± (and N sl = A t ) si fermions, and M Si+1/2 = JV t - JVj, 
(and M s _\j2 — [N± — iV-j-]) unbound +1/2 spinons (and unbound —1/2 spinons). 

Our above analysis confirms that for densities in the range n £ [0, 1] and m £ [0, 77] the ground states are inside the 
BA solution subspace. After this has been confirmed, analysis of the BA solution spectra straightforwardly leads to 
the ground-state c and si band occupancies given in Section |VT1 

Next we derive the unbound 77-spinon energies for 77 7^ 1. As a result of the 77-spin SU{2) symmetry of the 
Hamiltonian H sym m, its ground-state energy Eq S = Eq S (S v ) is the same for the LWS and HWS ground states for 
which S n — —S* 3 and S v = S* 3 , respectively. Indeed for that Hamiltonian all 2S V non-LWSs generated from a S v > 

ground state have the same energy. In contrast, the ground-state energy of the Hamiltonian H of Eq. (IB1[) depends 
both on the S v and S* 3 values, as given in Eq. ()B2|) . 

The continuous character for N a — > 00 and 77 7^ 1 of the chemical-potential curve /1 = ^(n) of Eq. (|C1[) of Appendix 
[C] implies that, except for terms of order 1 /N a , the following relation between ground state energies of the Hamiltonian 
H provided in Eq. (JBTJ holds, 

E GS (S r , + l/2,S* 3 - sgn{x}l/2) = E GS (S v ,S* 3 ), x = (1 - 77) , S v > . (B9) 

On combining Eqs. (|B2[) and (|B9|) we confirm that, 

M = E GS (S n + 1/2) - E° GS (S V ) > , x £ [-1, 0] , 

= E° GS (S n ) - E° GS (S V + 1/2) < , x £ [0, 1] . (B10) 

Indeed, when expressed in terms of the electron numbers rather than of the 77-spin values, such relations are equivalent 
to the usual chemical-potential definition, 

A 4 = e gs\n+i - E GS \ N . (Bll) 
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Furthermore, it follows from Eq. (|C1[) of Appendix [Cl that the sign relation sgn{/i} = -sgn{x} where x = (1 — n) holds 
for n ^ 1, so that, 

\ f i\=E° GS (S v + l/2)-E° GS (S v ), 5„>0, n^l. (B12) 

The relation Egs(S v , —S v ) = Egs{S v — SS V ,—S V + 6S V ) given in Eq. (|B4|) for SS n < S v between the energies 
of two n < 1 ground states reveals that creation or annihilation onto a n < 1 ground state of a finite number 
= 2<5 S 1 !) of unbound 77-spinons of 77-spin projection +1/2 is a vanishing-energy process. A similar relation 
Egs(S v , S v ) — Egs{St) — SS n , S n — SS n ) between the energies of two n > 1 ground states where SS V < S n confirms 
that creation or annihilation onto a n > 1 ground state of a finite number -M^^ , 2 — 25S n of unbound 77-spinons of 
77-spin projection —1/2 is also a vanishing energy process. 

On the other hand, a n < 1 (and n > 1) ground state state is a LWS (and HWS) of the 77-spin algebra. Hence 
it has no 77-spinons of ?7-spin projection —1/2 (and projection +1/2). A n < 1 ground state \GS) = \^i .io A . u ) °f 
ry-spin Sn = —S* 3 — (x/2)N a has M™ — M^ n +l ^ 2 — 2S V unbound 77-spinons of 77-spin projection +1/2. Those 
describe the 77-spin 5/7(2) symmetry degrees of freedom of a number 2S V of rotated-electron unoccupied sites. The 
corresponding degrees of freedom of the c hidden U(l) symmetry are described by an equal number Nf? = 2S V of c 
fermion unoccupied sites. 

We apply onto that n < 1 ground state the off-diagonal generator of the 77-spin algebra given in Eq. (jUJ). Such a 
process flips the 77-spin projection of one unbound 77-spinon. It transforms the ground state \GS) = \&i a jQ onto an 
energy eigenstate | x E'/ ,/ A , u ) with one —1/2 unbound 77-spinon, n n — M^ n _ x ^ 2 = 1. Within the occupancy configurations 
of that state, one rotated-electron unoccupied site has been replaced by a rotated-electron doubly occupied site. Such 
a process affects the rotated-electron site occupancy 77-spin SU(2) symmetry degrees of freedom. The corresponding 
degrees of freedom of the c hidden U(l) symmetry are not affected. The number of c fermion unoccupied sites and thus 
of c fermion holes are preserved. There are neither changes in the spinon occupancy configurations associated with the 
spin SU (2) symmetry degrees of freedom. Analysis of the corresponding BA solution energy spectra, confirms that 
no energy contributions arise from the c fermions and spinons, whose occupancy configurations remain unchanged. 

Furthermore, by the well known properties of the 77-spin SU(2) symmetry algebra, for the Hamiltonian H symm 
of Eq. @ such a rotation in the 77-spin space does not change the energy. Hence (^i j« & , u \Hsymm\'&i 0l ia &iU ) — 
(*/ ,; A ,u|-ff S ymm|*/ ,/A,«) where the energy (^^^JH^^^J is that of the initial ground state, \GS) = 
\^i a ,i° & ,u}- The operator commutes with the electron - rotated-electron unitary operator. Thus such a process 
enhances the number of both electrons and rotated-electrons by two. From the use of Eq. (]F33|) one then finds, 

(*J ,Ja,«I-H1*Jo,1a,«) = (^lo,iA,u\H svmm Wia^,u) -2fi= (^i .i & .u\H symm \^i o .i & . u ) + 2|/j| , n<l. (B13) 

In the above process, one unbound +1/2 77-spinon is annihilated and one unbound —1/2 ?7-spinon created. We have 
shown above that annihilation and creation of one unbound +1/2 77-spinon is a zero-energy process for an initial n < 1 
ground state. Hence it follows from Eq. (|F3 13|) that relative to that ground state, creation onto it of the unbound 
77-spinon of 77-spin projection —1/2 and corresponding rotated-electron doubly occupied site is a process of energy 

21m!- 

If instead one applies the 77-spin algebra generator S v of Eq. (|4ip onto a n > 1 ground state \GS) — |^ r ; ,jo , u ) of 

77-spinon S n = S^ 3 = (— x/2) N a and thus M~ n = M^.-i^ = 2S n unbound 77-spinons of 77-spin projection —1/2, one 
finds similar results. Specifically, we find that relative to such a ground state, creation onto it of the unbound +1/2 
77-spinon and corresponding rotated-electron unoccupied site is a process of energy 2\^l\. One then concludes that the 
energy £^,±1/2 for creation onto a n + 1 1 ground state of one unbound 77-spinon of 77-spin projection ±1/2 is given by 
Eq. (|87l) . Our analysis focused on the 77-spin SU(2) symmetry algebra. A similar analysis involving the spin SU(2) 
symmetry algebra leads to the unbound spinon energy provided in Eq. I|89p. 

Finally, we derive the n = 1 unbound 77-spinon energies from analysis of the energy spectrum associated with 
creation and annihilation of one electron onto and from the n = 1 and S n = ground state, respectively. The 77-spin 
rotations considered above generate changes in S* 3 that are restricted to integer values. In turn, the transitions 
between ground states whose electron numbers differ by one considered here are associated with half integer ±1/2 
deviations in S^ 3 . The well-defined discontinuity at n = 1 in the chemical-potential curve fi = /j(ti), Eq. (|C1|) of 
Appendix [C] reveals that the energy eigenvalues of the Hamiltonian of Eq. (|B1[) two ground states with numbers 
Sj) = 0, Sf t 3 — and 1/2, — sgn{x}l/2, respectively, are related as follows, 

£ G s(1/2,+1/2) = £ gs (0,0) + (m ±m), S v = 0, M e [-/j ,/]. (B14) 

Here /i° is the one-electron gap and 2/i° the corresponding Mott-Hubbard gap given in Eqs. (|C3[) and (|C4[) of Appendix 
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The transitions between ground states whose electron numbers differ by one affect all three degrees of freedom, 
associated with c hidden U(l) symmetry and two SU(2) symmetries. Creation of one electron onto such a ground 
state involves annihilation of one c fermion, annihilation of one spinon, and creation of one —1/2 ?7-spinon. On the 
other hand, annihilation of one electron from a n — 1 ground state involves annihilation of one c fermion, annihilation 
of one spinon, and creation of one +1/2 ?7-spinon. 

For simplicity, we consider initial ground states of zero spin density, m = 0. The final unbound ?7-spin energy 
expressions are the same at m = and for finite m. Only the energy scales fi and /j° in such expressions have m 
dependent values. The unbound ry-spin energies considered here arc relative to the zero-energy level of the Hamiltonian 
H, Eq. (jBip . (Due to symmetry, relative to that of the Hamiltonian H sym m of Eq. ([2]) such energies vanish.) The 
use of the exact BA solution reveals that the above annihilation of one c fermion and one spinon involved in the 
ground-state - ground-state one-electron excitations considered here are in the thermodynamic limit and relative to 
the zero-energy level of the Hamiltonian H of Eq. (|B1[) zero-energy processes. It then follows that that Hamiltonian 
energy eigenvalues deviations under such one-electron addition and removal excitations correspond to the energy for 
creation of the —1/2 7/-spinon and +1/2 r/-spinon, respectively. 

We consider two types of excitations, (i) Addition of removal or one electron. And (ii) addition or removal of one 
electron plus one 77-spin rotation. The latter rotation has as initial state a 77-spin 1/2 ground state. It is populated 
by a single unbound 77-spinon of 77-spin projection +1/2 or —1/2. The ?7-spin rotation flips it to —1/2 or +1/2, 
respectively. In the derivation of the energy spectra of the two types of excitations relative to the zero-energy level 
of the Hamiltonian H of Eq. (|B1[) the ground-state-energy relation of Eq. (|B14[) is used. One then readily finds that 
relative to the zero-energy level of that Hamiltonian the energy for creation onto a n = 1 ground state of one unbound 
±1/2 77-spinon is that given in Eq. (|88|) . Finally, the relations given in Eqs. (|90|) and (19T1) follow directly from the 
expressions provided in Eqs. (|57|) - (|M| . 



Appendix C: Useful energy scales 



In this Appendix we introduce several energy scales that are extracted from the BA solution and play an important 
role in the studies of this paper. 

The chemical potential curve \i = (J-(n) and magnetic energy scale 2fis H — 2fis H{m) associated with the magne- 
tization curve are defined in Eq. (|163l) . The corresponding chemical-potential dependence on the hole concentration 
x = (1 — n) is such that, 

n(x) = -fj,(-x); n(x) € [/A /J 1 ] , EG [-1,0]; [i G [—(J?, (J?] , x = 0, 

+/1 1 = lim fi(x) ; +/_t° = lim fi(x) ; /a 1 > fi° . (CI) 

The finiteness for u > of the one-electron gap, fx , implies that the chemical potential curve \i = /i(x) has a well- 
defined discontinuity at n = 1. For u 3> 1 and m — the chemical potential curve is such that the energy scale 
2/i = 2/i(x) is of the form, 

2[i{x) = — sgn{x}[[/ — 4i cos(7ra;)] , 16 [—1, 0] and x G [0, 1] , 

€ [-{U -4t),({7 -4t)], x = 0. (C2) 

For u > and n = 1 the Mott-Hubbard gap remains finite for all spin densities m 6 (—1, 1). It is an even function 
of m. For instance, at to = and m = —1,1 it is given by Q, 



^ J w(l + e" 2 ") U J 1 sinh(2^) 

= v/(4i) 2 TIP -At, m = -l,l, (C3) 
respectively, where Ji(w) is a Bessel function. For u <§; 1 and u>l, this energy scale behaves as, 



El 

8t ' 

2/i° rj [U-U], me [-1,1], u>l. (C4) 



o 

2/j° « -Vlf/e-Mij) j TO = 0; 2/j°w^— , to = -1,1, it < 1 , 



The energy scale 2/i 1 associated with the minimum — fi 1 and maximum chemical-potencial values is for all m 
magnitudes given by, 

2/j 1 = [U + At] . (C5) 
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On the other hand, the magnetic energy scale 2fi B H in Eq. (|163[) dependence on the spin density to is such that, 



2n B H(m) = -2fx B H(-m); 2fi B H (to) G (0, 2/i B H c ) , m G (-[1 - 0) , 
2fi B H(0) = ; T^fJ-B H c = lim 2fi B H(m) , x = (l-ri). 

x^±\l — \x\] 



(C6) 



A closed-form expression for the dependence on U, t, and density n of the energy scale 2a B H c associated with the 
critical field H c can be derived from the general 2[i B H expression provided in Eq. (| 163[) 59]. Since 2fi B H c is a even 
function of the hole concentration x = (1 — n), we expressed it in terms of it , 



2fi B H c 



u 2 



1 H — arccot 
7T \ 



c/ 2 



4* 



t/ (1 — \x\) H cos(7rx) arctan 



4* sin(7r|x|) 



tan(7r|a;|) 



(1 



(C7) 



At it = and for u 3> 1 its limiting behaviors are, 



2fi B H c = 4tsin^— 
8(l-|a:|)f 



2 r 



U 



u = , a; = (1 — n) , 
sin(27r|a;|) 



2tt(1 - |.t|) 



u > 1 



(1 - n) . 



(C8) 



As a function of n, it has for instance the following values, 

2fi B H c = 0, 7i = 0,2, 
1 
2 



C/ 2 - {/ 



1 3 

71 = — , — 

2' 2 



= yj{tt) 2 + U 2 ~U, n = 1 



(C9) 



Other energy scales involved in our studies are the c and av fermion energy dispersions given in Eqs. (I161[) and 
(|162p . Both the momentum widths 27r(l — n) and 27tto of the r\v and sv ^ si momentum bands, respectively, and their 
energy-dispersion bandwidths, [£ av {q a u) — e a^(0)L where av = r/v and av = sv ^ si, vanish in the n — > 1 ; m — > 
limit. Hence for the n = 1 and m = ground state the corresponding av ^ sl fermion energy dispersions defined 
in Eqs. (|161l) and (|162j) do not exist. On the other hand, the c fermion energy dispersions e®(q) and s c (q) and si 
fermion energy dispersions £®i{q) and £ s i(q) have closed-form expressions, which read, 



e°M = -| -2* cos fc°(<?)-4* 

Ec(Q') = e°(<?)+jU°, 9 e [— 7T, 7r] , 

cos(w A^g)) 



cos(a; sinfcg(q)) 

- -T7777T^ Jl(U) , <? € [— 7T, 7rJ 



-2* 



w cosh(w M 



JiH, ge [-7r/2,7r/2] 



(CIO) 



In this equation, Ji(w) is a Bessel function and the n = 1 and to = ground-state rapidity functions k®(q) (such that 
A®(q) = sin k®(q)) and A° 1 (g) are defined in terms of their inverse functions as follows, 



q = 
q = 



sm(u sinfcg(q)) 
aw - ^— - Jq[oj) , g G [-7T, 7rJ , 



') 



, sin(a; A", (q j) T . . r 
— v , J Q (w) , q G [-7T/2, tt/2] . 



co cosh(w u) 



(Cll) 



Here Jo(w) is again a Bessel function. 

On the other hand, for m — > and n € [0, 1] all sf ^ sl bands momentum and energy bandwidths vanish. For a 
to = ground state the corresponding sv sl fermion energy dispersions defined by Eqs. (I16ip and (|162l) do not 
exist. For densities n € [0, 1] and to = the c fermion, sl fermion, and rjv fermion energy dispersions £ c (q), £si(q), 
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and e° l/ (g), respectively, defined in such equations have the following limiting behaviors for u — > and u 1, 

£c(q) = -4i cos (I)- cos (Jy) , M < 2fc F = 7rn, u , 

, 2fc F = 7r n < \q\ < 7r , u — > , 



= -2* 



to\ /to 
cos ( \q\ ~ -yj ~ cos i^ — 



= —2t [cos(g) — cos(-to)] 
8nt 2 



U 



ln(2) [ sin 2 (q) -sin 2 (7 



\q\ <tt, it> 1 



(C12) 



e,i( 9 ) = e2 1 ( 9 ) = -2t 
27rnt 2 

1 



?7 



/■ro\] 

(g)-cos(^ — JJ , |g| < fc F = irn/2, u^O, 
s (~) ) | 9 | < fci? =7rn/2, u>l, 



sin(27r n) 



2tt n 



(C13) 



£ °„(g) - At cos fM±I?^ , | 9 |<( 7 r-2fe i r)= 7 r(l-n), u^O. 



8t 2 (l - n) 
vU 



1 - 



sin(27r(l - n)) 
2tt (1 - n) 



cos 



2(1 -n) 



, \q\ < (tt - 2fc F ) = tt (1 - n) , u»l. (C14) 



Also in the limit of an initial fully polarized ground state corresponding for instance to m — > n for electronic 
densities n G [0, 1], the energy dispersions e°(q) and e c (q), e° sv {q) and £ si ,(g), and £^(g) and £ V v(q) defined by Eqs. 
(I16ip and (|162l) have closed-form expressions, 



e°(q) = -^-2tcosg, g G [— w,n], 

£ c {q) = — 2t[cos q — cos(to)] , q € [— ir, 7r] , 

, x 2f /™ „ , , /[sinfc- A° (g)l 
= / dfc sin fc arctan 

7T V-™ \ VU 



q G (— 2fc F , 2k F ) = (-to, to) , 



£ S v{q) = + VK sly , 5 G (-2fc F ,2fc F ) = (-to, to) , 

e°„(g) = -^CZ + ^Rcfyi - 



2£ r" / [sinfc -A°„(g)] 

— / dfc sin fc arctan 1 

7T \ vu 



gG (-[7r-2fc F ],[7r-2fcir]) = (-7r[l-n],7r[l-n]), 



euv(g) = + gG (-[tt- 2fc F ],[7r-2fc F ]) = (-7r[l-n],7r[l -n]) 



(C15) 



In this limit one has that k®(q) — g, whereas the rapidity functions A°„(g) appearing in the above expressions are 
defined by their inversion functions, 



q = 2Re[arcsin(A° (qj) + ivu)\ 

1 r" , / [sinfc -A^(g)] 
afc arctan : 
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g G (-[tt - 2fc F ], [tt - 2k F \) = (-7r[l -n],7r[l - n]) , 



1 f™ „ /[sinfc- A°(g)l\ 
g = / dk arctan siA^i ^ g (_ 2 fc F , 2fc F ) = (-to, to) . 



7T 



(C16) 



Moreover, for the electronic density range n G [0, 1] and spin density m — > n the sv and r\v energy dispersion 
bandwidths of Eq. (|164[) appearing in some of the expressions provided in Eq. (|C15[) are given by, 



tt y (yU) / 



4t ( 4£sm(7m) , 

vU n cos(7rn) arctan — , n G [0, 1J , m — > n . 

7T \ vU 



(C17) 
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and 



2 



H — arccot . 

7T \ {VU) 



TT , 4t / 4f sm(7rn) 

(1 — n) H cos(7rn) arctan 

7T \ i/(7 



tan(7rn) 
ti 6 [0, 1] , m — > n , 



(C18) 



respectively. The energy bandwidths W s „ and given here are an increasing and decreasing function of n, respec- 
tively. For instance, for n — > 0, n = 1/2, and n = 1 they read, 



and 



, n — » , m — > n . 



n = - , m — > ra , 
2 



(C19) 



n — > , m — > n , 



V(4i) 2 + (i/(7) 2 - (7 



, n = 1 , to — > n . 



n = — , m — > n . 
2 ' 



(C20) 



respectively. 

For spin density to = and electronic density n = 1 all energy bandwidths vanish except that of the si fermion 
energy dispersion, which reads, 



W sl = 2t 



duj- 



uj cosh(w u) 



n — 1 , m = . 



(C21) 



Finally, for u — > and electronic densities n € [0,1], and spin densities to = and to — > n the sv and 771/ 

dispersion energy bandwidths of Eq. (|164[) have the following limiting behaviors, 



2t 



1 — cos 



TO = , U . 



9 /7TO\ 

4t sin I — ) = 77 c , to — > n , w -> , 



= Vi 



irn' 



U 



1 - 



(£»)]} 

B H c , 
sin(27rn) 



27T7T, 



= 5 v ,i - (2/is i? c ) , m = , it > 1 , 



i/l7 



1 - 



sin(27rn) 



and 



= 4t cos 1 —11 

= At 



2nn 



— 2/ibH c , to — > n , !i> 1, 



(C22) 



1 • 2 /7T"\ 

1 - sur — 



2|/u| , m = , u ->• 0, 

= 2|ju| , to — > n , it — > , 



2 7 



1 



8(1 -n)f 2 
i/l7 

8nt 2 sin(27rn) 
vU 2irn 



sin(27r(l - n)) 
2tt(1 - n) 

, m — > n , m> 1 



, m = , u ^> 1 , 



(C23) 



respectively. 
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